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Outline  of  the  Final  Report 


The  research  work,  which  has  been  accomplished  at  the  Regis 
College  Research  Center  for  AFGL  contract  F19628-80-C-0016 ,  includes 

(1)  adaptive  procedures  to  correct  ionospheric  variability  in  time 
delay  and  range  error  measurements  of  radio  waves  (2)  seasonal  and 
diurnal  studies  of  total  electron  contents  (TEC)  at  Ascension  Island 
(3)  effects  of  ionospheric  modification  caused  by  high  power,  HF  radio 
waves  (4)  Faraday  polarization  fluctuations  (FPF)  of  satellite  signals 
in  transionospheric  propagation. 

Multiplicative  and  additive  prediction  functions  are  used 
as  adaptive  procedures  to  correct  ionospheric  variability  in  time  delay 
and  range  error  measurements  of  radio  waves.  The  adaptive  algorithms 
employ  real  time  measurements  to  reduce  rms  errors  in  ionospheric  pre¬ 
dictions.  It  is  found  that  significant  reductions  in  rms  errors  are 
achieved  in  daytime,  at  solar  maximum,  when  time  delay  and  range  errors 
are  greatest,  using  a  multiplicative  prediction  function. 

TEC  data  recorded  at  Ascension  Island  from  September  22,1978 
through  December  21,1980  have  been  analyzed,  displaying  two  character¬ 
istic  features  of  the  night-time  ionosphere.  They  are  (1)  a  post¬ 
sunset  enhancement  which  lasts  for  several  hours,  and  follows  a  rapid 
decrease  unattributable  to  the  usual  ionospheric  decay  processes,  and 

(2)  the  sudden  sharp  changes  in  the  TEC  along  the  ray  path  of  satellite 
signals,  that  are  often  superimposed  on  those  enhancements  and,  to  a 
lesser  degree,  on  the  subsequent  normal  night-time  ionosphere.  F  layer 
dynamo  theory  and  generation  of  ionospheric  bubbles  are  invoked  to  in¬ 
terpret  these  observations. 

Many  unexpected,  nonlinear  plasma  processes  can  be  triggered 
in  the  ionosphere  by  a  powerful,  HF  radio  wave  near  the  F  peak.  Selec¬ 
ted  for  studies  are  the  following:  ion  line  enhancement,  the  excitation 
of  upper  hybrid  modes  or  field-aligned  plasma  lines,  and  the  induced 
field-aligned  ionospheric  irregularities  with  a  broad  scale-length 
range.  Various  parametric  instabilities  are  responsible  for  these 
phenomena.  Beating  current  at  zero  frequency  is  also  an  important 
nonlinearity  in  producing  field-aligned  high-frequency  modes  (upper 
hybrid  waves)  and  large-scale  non-oscillatory  modes  (irregularities) 
in  addition  to  the  differential  Ohmic  heating  force  and  the  nonlinear 
Lorentz  force.  In  the  analysis  of  non-oscillatory  modes,  it  is  more 
reasonable  to  employ  four-wave  than  to  use  three-wave  interaction 
process.  A  criterion  has  been  established,  allowing  the  excitation 
of  upper  hybrid  waves  by  radio  waves  with  ordinary  rather  than  ex¬ 
traordinary  polarization.  Ionospheric  inhomogeneity  can  be  shown  to 
impede  generally  the  excitation  of  plasma  instabilities  and  consequent¬ 
ly  suppress  the  growth  of  the  induced  ionospheric  irregularities.  The 
thresholds  of  various  parametric  instabilities  concerned  are  rather 
low  as  compared  to  the  power  typically  used  in  the  ionospheric  heating 
experiments. 

Faraday  polarization  fluctuations  (FPF)  of  beacon  satellite 
signals  origin  in  the  differential  phase  shift  and  the  differential 
changes  in  logarithmic  amplitudes  of  ordinary  and  extraordinary  modes 
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when  a  linearly  polarized  radio  wave  propagates  in  a  birefringent 
medium  (i.e.,  the  ionosphere).  FPF  has  been  observed  for  signals 
transmitted  from  either  low-flying  orbitting  satellites  or  geosta¬ 
tionary  satellites  in  the  HF  and  VHF  bands.  Ionospheric  irregularities 
are  identified  to  be  the  scatterers  of  satellite  signals  in  trans- 
ionospheric  propagation,  that  lead  to  the  FPF.  Ionospheric  irregu¬ 
larities  of  power-law  type  and  quasi-Gaussian  type  correspond  res¬ 
pectively  to  the  early  phase  and  the  late  phase  of  ionospheric  distur¬ 
bances.  It  is  shown  that  power-law  type  irregularities  can  cause  the 
FPF  of  satellite  signals  with  frequencies  up  to,  approximately, 

250  MHz  under  the  typical  ionospheric  conditions  at  Ascension  Island. 

In  contrast,  Gaussian  type  irregularities  can  only  effect  the  HF 
satellite  signals,  e.g.,  20  MHz.  Although  the  saturation  of  amplitude 
scintillation  of  radio  waves  is  due  to  the  multi-scattering  process, 
it  is  accurate  enough  to  use  the  single-scattering  process  to  describe 
the  physics  of  FPF  even  in  the  VHF  band. 


I*'-. 


4% 


2 


Radio  Seine e,  Votaac  16,  Number  2,  pa*et  261-269.  March-April  198! 


Tiae  cefis  for  adaptive  prediction  of  total  electron  content 

D.  E.  Donatelli 

Regis  Colt* ge  Research  Center,  Weston,  Massachusetts  02193 

R.  S.  Alien 

Air  Fore*  Geophysics  Laboratory,  Ha  ns  com  Air  Force  Base,  Massachusetts  01731 
(Received  March  20.  1900,  reviled  October  2.  1980;  accepted  October  6,  1980.) 


Adaptive  procedures,  uata|  multiplicative  aad  additive  prediction  f unctions  of  the  standard  deviation 
of  the  observations,  arc  aaalyred  to  estimate  and  compare  their  effectiveness  in  correcting  for 
ionospheric  variability  in  time  delay  and  range  measurements  The  adaptive  algorithms  use  real 
time  measurements  to  reduce  rms  errors  in  ionospheric  predictions.  Total  electron  content  data 
from  Hamilton.  Massachusetts,  for  1969-1976  are  used  to  evaluate  diurnal,  seasonal,  and  solar 
cycle  impact  on  radar,  navigation,  and  communication  systems.  Time  cells  are  defined  from  the 
growth  rate  of  rms  residual  error.  Significant  reductions  in  rms  error  are  achieved  in  daytime, 
at  solar  maximum,  when  time  delay  and  range  errors  are  greatest,  using  a  multiplicative  prediction 
function. 


INTRODUCTION 

The  potential  precision  of  radar,  navigation,  aad 
communication  systems  may  be  increased  through 
advances  in  device  technology,  but  daily  variability 
of  the  ionosphere  continues  to  be  a  constraint  on 
achievable  system  performance.  Numerical  maps 
that  provide  monthly  median  correction  for  iono¬ 
spheric  effects  have  been  derived  from  a  worldwide 
climatology  of  ionospheric  parameters;  their  use 
alone  reduces  the  monthly  rms  residual  error  to 
about  20-23%  of  the  median  correction  in  daytime 
and  30-35%  at  night  [ Du  Long ,  1977] .  Many  preci¬ 
sion  systems  cannot  tolerate  the  individual  depar¬ 
tures  of  30-200%  that  are  averaged  away  in  monthly 
statistics.  For  these  systems,  adaptive  prediction 
algorithms  are  presented  using  field  measurements 
to  provide  real  time  correction  of  ionospheric  ef¬ 
fects. 

The  local  ionosphere  in  a  system  coverage  area 
may  be  monitored  through  measurements  of  iono¬ 
spheric  parameters,  such  as  the  critical  frequency 
of  the  F  region  (fOF2),  total  electron  contect  (TEC), 
the  gradient  in  the  lower  F  region,  or  through 
measurements  derived  from  the  system  itself. 
Since  some  system  users  may  not  have  the  capability 
of  implementing  the  computer  models  based  on 
worldwide  climatology,  they  may  use  locally  moni- 
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tored  TEC  or  f0F2  to  estimate  monthly  median 
behavior  and  short-term  persistence  effects  by  pro¬ 
ducing  means  of  several  days  to  predict  the  subse¬ 
quent  hourly,  daily,  and  weekly  response  of  their 
system  to  the  current  state  of  the  ionosphere. 

Adaptive  prediction  algorithms  have  been  exam¬ 
ined  using /0F2  data  [Wilkinson,  1979]  and  TEC 
data  [Donatelli  and  Allen,  1978;  Leitinger  et  al, 
1978;  Donatelli  and  Allen,  1979] .  These  include  TEC 
as  derived  from  measurements  of  Doppler  shift 
using  passes  of  the  Navy  Navigation  Satellites 
(NNSS)  and  from  Faraday  rotation  measurements 
using  VHF  beacons  from  geostationary  satellites. 
The  procedures  apply  to  any  type  of  synoptic  data, 
and  rely  on  a  reasonable  prediction  of  the  monthly 
mean  of  the  observations,  the  standard  deviation, 
and  an  assumption  on  the  correlation  of  the  ob¬ 
servations  in  the  time  cell  being  considered.  Two 
basic  approaches  will  be  comparatively  analyzed: 
the  first  is  a  multiplicative  prediction  function,  and 
the  second,  an  additive  prediction  function.  The 
second  approach  has  been  analyzed  in  detail  by 
Gautier  and  Zacharisen  (1965)  with  respect  to/0F2 
data  and  by  Da  Rosa  (1974]  with  respect  to  TEC. 
An  evaluation  of  an  adaptive  procedure  using  a 
multiplicative  function  to  update  a  monthly  mean 
refraction  correction  with  real  time  measurements 
of  TEC  is  presented  here.  The  temporal  growth 
of  the  rms  residual  error,  within  a  time  cell  extending 
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from  a  reference  (ime  with  residual  error  defined 
as  zero,  to  an  upper  limit  defined  to  lie  within  the 
error  budget  of  a  particular  system,  is  examined. 

The  ionospheric  time  delay  of  a  radio  wave  is 
directly  proportional  to  the  electron  content  along 
the  ray  path,  such  that  the  results  presented  here 
in  TEC  may  be  translated  by  the  user  to  system 
parameters  as  follows: 

Time  delay  ( Davies ,  1966] : 

40.3 

AT  =  — rxTEC  (sec) 

cf 

Range  error: 


of  the  applied  correction.  The  predicted  value  is 
defined  by  the  function  being  considered,  i.e., 

Case  1,  multiplicative  function: 

A  =  a.  —  M-,  (2  a) 

►*■'0 

Case  2.  additive  function: 

A,  =  U,  +  a2(0,#  -  p,,„)  (26) 

where  p.  is  the  observed  mean,  r0  is  the  reference 
time,  and  a ,  and  a2  are  constants  to  be  determined. 
Substituting  for  P,  in  (1),  the  residual  error  is: 

Case  I: 


40.3 

AA  *  cAT  =  — —  x  TEC 

r 


(m) 


where  /  is  the  frequency  of  the  radio  signal  in 
megahertz,  c  is  the  speed  of  light  in  a  vacuum 
in  meters  per  second,  and  TEC  is  in  units  of 
electrons  per  square  meter  times  10  “. 


-  2e,  — 

•S 

Case  2: 


(3a) 


ANALYSIS 


The  basis  for  the  procedure  is  the  assumption 
that  major  components  of  daily  variability  in  the 
ionosphere  can  be  described  as  large-scale,  large- 
magnitude,  slowly  varying  fluctuations  and  are, 
therefore,  reasonably  correlated  in  time.  These 
fluctuations  result  from  daily  variables  such  as  solar 
activity,  magnetic  activity,  lunar  tides,  neutral  winds 
that  could  produce  multiplicative  or  additive  effects 
on  the  ionosphere.  The  rms  error  using  the  adaptive 
procedure  is 


by  expanding 


xl  *  a*  +  -  2flJa,<T,()p,  (36) 

where  p,  is  the  correlation  coefficient: 

a  ~  H.XO,.  ~  »*.„) 

*"(X(O,-pl),I(O,0-p,.),)'/1 

*0.0.,  A 

0,0, a  XNh,?.,'  ) 

Differentiating  x,  with  respect  to  a,  and  a2 ,  setting 
equal  to  zero,  and  solving  for  a,  and  at 

Case  1: 

o]o,  p,  +  n,1 

®i  *  5  J"  *7" 

0,0,  +  I L, 

Case  2: 


X,1  -  *?  +  +  —  (A? -20, A.) 

N 

where  O,  is  the  observation  at  the  time  when  the 
correction  is  applied,  P,  is  the  prediction  at  the 
time  of  correction  using  the  observation  from  the 
reference  time,  o,  is  the  standard  deviation  of  the 
observations: 

o]  -  <; i/NHO ,  -  n,)1  -  -  o?  - 
N 

S  is  the  number  of  observations,  and  t  is  the  time 


p, 


where  a,  and  3,  are  defined  as 

o,  ■ - 

o, 

By  substituting  for  a,  and  a,  in  (3a)  and  (36), 
respectively,  the  minima  for  x,  are  obtained: 
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Case  1: 


P?  + 


P?(°,  ~  P,)}  ] 

°2,  «*?  +  P?  J 


Case  2: 


(4a) 


X!«a?(l-p?)  (46) 

Equation  (46)  is  the  result  obtained  by  Gautier  and 
Zackarisen  [1965].  Equation  (4a)  reduces  to  (4b) 
if  a,  =  p,. 

To  reduce  the  average  monthly  residual  error  such 
that  x,  s  co,,  where  c  is  the  degree  of  reduction 
desired,  the  implied  condition  for  case  .  is: 

1  -  pj  +  (a,  -  p,)J(-j — ^jst1  (5a) 

which,  for  a,  =  p(,  reduces  to  the  condition  for 
case  2: 


l-plst1  (56) 

If  a,  p,,  the  quality  of  the  prediction  degrades 
proportional  to  the  difference  between  a,  and  p,. 

To  minimize  the  error,  the  correlation  coefficient 
must  be  known.  If  p,  is  not  known,  a,  and  a2  cannot 
be  defined.  Therefore,  to  provide  an  operational 
definition  for  P,,  a,  and  a2  in  (2a )  and  (26)  are 
arbitrarily  set  equal  to  one.  Then  (3a)  and  (36) 
can  be  written  as 


Case  1: 

X 2  “  [1  -p?  +  (<*,-p,)2] 

(6a) 

Case  2: 

X,”  -  o?  |l -P,* +(».-*,)’) 

(6b) 

Comparing  (6a)  with  (4a)  and  (66)  with  (46)  it  is 
seen  that  in  case  1,  x)  increases  by  elimination 
of  the  factor  (p.2)/(<r?#  +  p.,20)  in  the  third  term 
of  (4a),  and  case  2  increases  by  addition  of  a  third 
term,  O,  -  p,)2.  The  limiting  condition  for  case 
1  is 

1 -p?  +  (a,-p,)JscJ  (7) 

which  defines  limits  for  a,  and  p,: 

a,sP,±  Ipf  -  (1  -c3)],/J  (8) 

Since  o,  is  always  real  and  positive,  the  condition 
p*  2  1  -  c2  &  0  must  always  be  satisfied.  By 
definition,  c  s  1;  therefore  0  <  a,  <  2p, .  (The 
analysis  for  case  2  is  identical  to  case  1  with  (3, 
substituted  for  a,.)  For  example,  if  c  ■  0.5,  which 


implied  a  50%  reduction  in  residual  error,  p,  2s 
0.866  is  required.  However,  if  p,  =  0.5,  a  reduction 
in  residual  eiror  of  more  than  10%  is  possible  if 
a.O,)  -  P, 

The  TEC  data  used  to  evaluate  this  procedure 
were  reduced  from  measurements  of  the  Faraday 
rotation  of  signals  from  the  beacon  of  the  ATS 
3  geostationary  satellite  observed  at  Hamilton, 
Massachusetts,  during  the  1969-1976  period  cover¬ 
ing  the  maximum  and  minimum  of  the  last  solar 
cycle.  The  location  of  the  420-km  subionosphere 
point  [ Titheridge ,  1972]  along  the  ray  path  from 
Hamilton  is  39°N,  289°E.  These  data  provide  a 
basis  for  evaluation  of  average  solar  cycle  behavior, 
since  this  solar  cycle  closely  approximates  the 
average  of  the  past  20  solar  cycles  with  Rt  =  100, 
£  =  155  at  maximum,  and  R,  =  10,  S  =  71  at 
minimum,  where  Rt  and  5  are  the  12-month  running 
mean  sunspot  number  and  solar  flux  at  2800  Mhz, 
respectively.  To  estimate  TEC  for  other  values  of 
Rt  (or  S),  it  is  possible  to  interpolate  to  a  reasonable 
approximation  by  assuming  a  linear  relationship 
between  Rt  (or  S)  and  TEC  [ DuLong ,  1977] . 

By  using  the  definition  of  a,  and  setting  p,  = 
1,  the  potential  of  an  adaptive  procedure  using  a 


Fi*  ••  Diurnal  curves  of  percentage  standard  deviation  in 
TEC  for  the  months  of  January.  April.  July,  and  October  1969 
at  Hamilton,  Massachusetts. 
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Fig.  2.  Sues  u  for  Figure  1  but  for  197S. 


multiplicative  function  can  be  estimated  with 
knowledge  of  the  percentage  standard  deviation  for 
the  reference  time  and  the  time  of  applied  correction 
(6a).  This  percentage  was  calculated  for  each  month 
of  data  and  is  shown  in  Figures  1  and  2  for  the 
months  of  January,  April,  July,  and  October  1969 
and  1973,  respectively,  representing  diurnal,  sea¬ 
sonal,  and  solar  cycle  effects. 

The  growth  of  the  percentage  error  in  time  at 
any  location  may  be  estimated  by  taking  progressive 
ratios  along  the  curve  from  some  choice  of  reference 
time  on  that  curve.  When  a,  *  c,  the  limit  of  a 
time  cell  is  defined,  since  the  residual  error  will 
then  equal  the  monthly  mean  error. 

One  of  the  limitations  on  time  cell  extent,  apparent 
from  these  curves,  is  that  large  errors  may  be 
generated  when  a  scaling  factor  from  nighttime  is 
used  across  the  sunrise  terminator  to  adapt  daytime 
predictions.  This  has  a  greater  effect  at  solar  mini¬ 
mum,  since  the  percentage  variability  fluctuates 
more  radically  at  this  time,  a  reflection  of  the  lower 
mean  values  at  solar  minimum  where  a  small  abso¬ 
lute  variation  in  TEC  will  produce  a  large  percentage 
variation,  particularly  at  nighttime.  In  general,  ex¬ 
tended  time  cells  are  possible  at  solar  mar  num 
but  must  be  greatly  reduced  *  e^l«r  r  .mum. 
However,  the  remaining  variabu»rv  in  >  .C  units 


((electrons/m2)  x  10>6),  x,<  is  still  less  at  solar 
minimum. 

When  an  additive  function  is  considered,  the 
potential  of  the  adaptive  technique  may  be  estimated 
using  Figures  3  and  4.  These  are  the  curves  for 
the  standard  deviation  in  TEC  units  corresponding 
to  the  percentage  curves  of  Figures  1  and  2.  Again, 
the  process  is  considered  for  p  *  1,  where  the 
error  x,  is  defined  by  (66),  and  is  estimated  in  TEC 
units  by  taking  ratios  between  points,  analogous 
to  Figures  1  and  2.  From  these  figures  it  appears 
that  additive  functions  may  well  describe  nighttime 
and  solar  minimum  conditions  but  would  be  less 
descriptive  of  daytime  conditions  at  solar  maximum. 

RESULTS 

Since  the  ionospheric  impact  on  systems  is  most 
severe  in  daytime,  at  solar  maximum,  the  adaptive 
procedure  applying  a  multiplicative  function  is  ex¬ 
amined  using  the  assumption  that  within  the  data 
base  the  correlation  is  expected  to  monotonically 
decrease  to  zero  with  time,  and  negative  correlation 
is  not  considered.  The  procedure  is  initiated  by 
obtaining  a  measurement  of  TEC  at  reference  time 
0,o,  determining  a  scaling  factor  from  the  ratio 


Fig.  3.  Diurnal  curves  of  standard  deviation  in  TEC  for  the 
mooths  of  January,  April,  July,  and  October  1969  at  Hamilton. 
Massachusetts 
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Fig.  4.  Same  u  for  Figure  3  but  for  1973. 


of  the  observation  to  the  baseline  prediction  for 
the  reference  point  M,o  and  scaling  the  baseline 
map  for  a  succeeding  time  period.  The  adapted 
prediction  P,  at  a  later  time  is  then  P,  = 
yo,o/Mlo)M,,  where  M,  is  the  baseline  prediction 
at  the  later  time.  As  noted  previously,  the  actual 
monthly  mean  is  used  for  the  baseline  prediction 
to  remove  any  model  dependent  error  from  the 
ensuing  results.  Thus  from  case  1,  P,  *  I#  )|t, . 

The  TEC  data  are  reduced  at  15-min  intervals. 
The  adapted  prediction  is  applied  at  these  intervals 
over  several  hours.  Comparisons  are  made  with 
the  corresponding  observations  and  the  rms  error 
X,  is  computed  each  month  at  each  1 5-min  increment 
from  the  same  reference  times  each  day.  The  results 
for  January  1969  are  presented  in  Figure  5,  including 
a  comparison  of  results  using  the  actual  monthly 
mean.  Figure  (5a),  and  a  model  predicted  mean 
based  on  ionospheric  climatology  [Barghausen  et 
at.,  1969] ,  Figure  5b.  The  outer  scale  on  the  figures 
is  the  mean  and  is  4  times  greater  than  the  inner 
scale,  which  is  the  standard  deviation  with  respect 
to  the  mean  and  the  rms  error  of  the  adapted 
prediction  with  respect  to  the  observations.  Al¬ 
though  there  is  up  to  20%  difference  between  the 
predicted  and  the  observed  mean,  the  standard 
deviation  with  respect  to  each  mean  is  comparable, 


indicating  that  without  adaption,  the  monthly  cli¬ 
matology  prediction  is  as  useful  as  the  mean  itself. 

The  growth  of  the  rms  error  of  the  adapted 
prediction  x,  is  described  by  the  thin  tines  originat¬ 
ing  at  2-hr  intervals,  with  x,  =  0.  In  (5a),  x,  extends 
in  time  from  t  =  r0  to  t  =  t0  +  3  hr;  in  (56), 
from  t  =  r0  to  r  «  <0  +  12  hrs.  The  scaling  factor 
has  been  determined  by  the  observation  at  t  -  t0 
and  held  constant  over  the  3-  and  12-hr  time  cell; 
only  the  scaling  of  the  mean  u,  and  the  predicted 
mean  M,  is  progressed  with  time  (Figures  5a  and 
56,  respectively). 

The  dashed  lines  define  the  residual  error  for 
a  3-hr  time  cell  (x,0*j*)  and  indicate  that  a  scaling 
factor  that  is  3  hr  old  is  useful  throughout  the 
daytime  hours,  unless  it  was  determined  prior  to 
sunrise  and  used  thereafter.  The  difference  between 
X>0+»  in  (^°)  and  (56)  is  the  model  dependent  error, 
most  evident  in  the  sunrise  area.  Comparing  these 
curves  with  the  standard  deviation  of  the  observed 
o,  and  predicted  mean  £,,  it  is  apparent  that  the 
time  cell  in  which  the  adapted  prediction  is  useful 
is  bound  by  the  reference  time  on  one  end  and 
the  succeeding  solar  terminator  on  the  other. 

An  analysis  of  x,  for  case  1,  showing  the  effect 
of  the  model  bias,  can  be  found  in  appendix  A; 
Gautier  and  Zacharisen  [1965]  have  presented  an 
equivalent  analysis  for  case  2.  It  has  been  assumed 
that  the  data  follows  a  normal  distribution .  In  reality , 
however,  it  is  somewhat  skewed  toward  higher 
values  in  daytime,  such  that  a  model  that  predicts 
a  slightly  greater  mean  than  that  observed  will  often 
succeed  in  reducing  the  monthly  variance  at  that 
time. 

The  performance  of  this  adaptive  technique  is 
summarized  in  Table  1  using  the  archive  data  of 
TEC  from  Hamilton,  Massachusetts,  for  1969 
through  1976.  The  daytime  and  nighttime  maxima 
of  the  mean,  standard  deviation,  and  three  time 
cells,  averaged  over  season,  at  both  solar  maximum 
and  minimum,  in  TEC  units,  (electrons /m2)  x  10‘\ 
are  compared.  The  table  shows  that  an  adapted 
prediction  can  reduce  the  residual  error,  on  average, 
by  60%  in  a  1 -hr  time  cell  and  still  as  much  as 
30%  in  a  3-hr  time  cell  for  daytime  solar  maximum. 
Detailed  results  may  be  found  in  DuLong  [1977] 

The  time  cell  is  theoretically  determined  by  the 
correlation  time  of  fluctuations  in  TEC.  Since  the 
basic  data  for  this  study  were  reduced  at  15 -min 
intervals,  fluctuations  with  periods  less  than  30  min 
are  not  observable.  Their  ampbtude  is  assumed  to 
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Fig.  5.  (a)  Comparison  in  TEC  for  January  1969  at  Hamilton.  Massachusetts,  of  the  observed  mean,  standard 
deviation,  rms  error  of  a  J-hr  time  cell,  and  error  growth  in  time  using  an  adapted  prediction  up  to  3 
hrs  after  a  reference  time.  (6)  Similar  to  5a  but  using  a  model  prediction  of  mean  TEC.  The  error  growth 
in  time  using  the  adapted  prediction  is  shown  for  12  hr  after  the  reference  time. 


be  much  smaller  than  the  longer-period  variations 
that  are  being  tracked  here.  The  present  simple 
adaptive  technique,  predicting  a  persistence  of  IS 
min  to  3  hr,  tends  to  mitigate  against  just  those 
long-period  TEC  variations,  which  seem  to  be  a 
major  component  at  solar  maximum,  with  a  residual 
of  the  same  absolute  magnitude  as  the  residual  at 
solar  minimum.  This  suggests  that  short-period 
variations  of  TEC  predominate  at  solar  minimum 
and  are  still  present  with  comparable  amplitude  at 


solar  maximum,  masked  by  the  characteristic 
larger-scale  variations. 

During  sunrise  and  sunset  periods,  the  effects 
of  production  and  loss  create  gradients  that  domi¬ 
nate  all  other  fluctuations.  Similar  gradients  occur 
during  magnetic  disturbance,  particularly  in  the 
region  of  the  aurora]  zone  and  the  trough  in  high 
and  mid-latitudes.  The  rate  at  which  these  changes 
occur  should  determine  the  time  cells  for  these 
periods.  This  i$  illustrated  in  Figure  6.  In  Figure 
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TABLE  1.  Summary  of  effectiveness  of  various  time  cells  in  reducing  prediction  error 


Winter 

Vernal  Equinox 

Summer 

Autumnal  Equinox 

Day 

Night 

Day 

Night 

Day 

Night 

Day 

Night 

34.4 

5.7 

Solar  Maximum  (A, 
44.8  8.2 

*  110;  S  -  ISS) 
25.4 

8.2 

41.1 

6  7 

«.o 

1.8 

7.8 

2.5 

4.0 

2.2 

7.8 

2.2 

)k 

4.0 

1.5 

4.5 

15 

2.7 

1.5 

4.0 

18 

XlQ*  III 

2.2 

1.0 

2.5 

1.0 

1.6 

1.2 

1.6 

0.7 

Xgg-e  Mm 

1.5 

0.6 

1.5 

0.6 

1.2 

06 

13 

0.6 

10.5 

1.5 

Solar  Minimum  (Rt 
10.5  1.5 

-=  10.  5  -  71) 
9.0 

2.1 

11.2 

2.1 

1.8 

0.6 

2.5 

0.6 

1.5 

0.5 

1.5 

0.6 

X<o«3h 

1.8 

0.9 

2.1 

0.5 

1.5 

0.5 

2  1 

06 

Xiq*  lb 

1.2 

0.3 

1.2 

0.1 

0.6 

0.3 

1.2 

0.3 

XgQ-v-Mhn 

0.9 

0.1 

0.6 

0.1 

04 

0.1 

0.6 

0  1 

Symbols  represent  the  following,  in 

TEC  units:  p,. 

the  mean  value  of  TEC;  a,, 

the  rms  deviation  of  TEC 

from  the  mean, 

representing  day-to-day  variability  of  the  ionosphere;  the  residual  error  using  an  adapted  prediction  in  a  3-hr  time  cell: 

the  residual  error  using  an  adapted  prediction  in  a  l-hr  time  cell;  the  residual  error  using  an  adapted  prediction 

in  a  30-min  time  cell. 


(6a)  the  observed  TEC  for  the  day  O,  is  compared 
with  the  prediction  for  the  month  p.,  and  the  adapt¬ 
ed  prediction  for  a  30-min  time  cell,  P,a+i0n,  and 
a  3-hr  time  cell,  P,  +Jh  ■  The  difference  between 


Fig.  6.  (a)  TEC  on  a  magnetically  disturbed  day  compared 
to  the  observed  TEC,  the  mean,  and  the  adapted  predictions 
for  the  day  using  30-min  and  3-hr  time  cells  (b)  Residual  error 
m  TEC  using  the  predictions  of  6a  for  the  day. 


the  observed  TEC  and  the  predictions  in  each  case 
are  shown  in  Figure  (66).  The  error  using  the  3-hr 
time  cell  |0Jo+3k  -  P,a.ih  |  is  equivalent  in  magnitude 
to  that  using  the  monthly  mean  | O,  -  p,|,  but 
large  excursions  have  been  introduced  at  time 
periods  when  they  would  otherwise  not  exist.  The 
adapted  prediction  using  a  30-min  time  cell,  P, (  „  30ni . 
provides  a  60%  reduction  in  the  maximum  error 
despite  the  steep  gradients  in  TEC  on  this  disturbed 
day. 

CONCLUSION 

When  using  adapted  predictions,  it  must  first  be 
determined  whether  the  dominant  features  are  de¬ 
scribed  by  additive  or  multiplicative  functions.  For 
instance,  for  a  system  propagating  through  the 
auroral  region,  large  changes  of  TEC  related  to 
precipitating  particles  may  be  dominant  sources  of 
error  that  are  best  treated  as  additions  to  a  back¬ 
ground  ionosphere,  whereas  for  a  system  at  mid¬ 
latitude,  large  changes  of  TEC  related  to  changes 
in  neutral  composition  might  best  be  treated  as  a 
multiple  of  the  expected  ionosphere.  From  this 
analysis  it  was  shown  that  the  rms  error  using  an 
additive  function  is  dependent  on  the  absolute 
difference  in  the  standard  deviation  of  the  two  points 
under  consideration.  For  the  multiplicative  function 
the  absolute  difference  in  percentage  standard  de¬ 
viation  with  respect  to  the  mean  is  the  relevant 
parameter.  Examination  of  these  parameters  indi¬ 
cated  that  for  the  solar  maximum  period,  when 
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the  ionospheric  impact  on  systems  is  most  severe, 
multiplicative  functions  best  describe  the  dominant 
features.  An  evaluation  for  an  adaptive  procedure 
that  corrects  for  these  features  has  been  presented. 
The  procedure  is  successful  if  the  condition  1  2 
p,  >  (a,/2)  2  0  exists,  where  p,  is  the  correlative 
coefficient  at  the  time  of  applied  correction  and 
a,  is  the  ratio  of  the  percentage  rms  deviation  at 
the  reference  time  to  that  at  the  time  of  applied 
correction. 

In  daytime  an  improvement  of  50 %  is  readily 
achievable,  at  solar  maximum  within  a  3-hr  time 
cell  and  at  solar  minimum  within  a  1-hr  time  cell. 
At  nighttime,  the  percentage  rms  deviation  is  usually 
greater  than  daytime,  but  the  absolute  deviation 
is  less;  therefore  the  need  for  correction  is  not 
as  critical.  However,  if  adapted  predictions  are 
used,  a  scaling  factor  obtained  at  nighttime  should 
not  be  projected  across  the  sunrise  terminator,  as 
the  rapid  change  in  percentage  rms  deviation  during 
the  sunrise  growth  period  introduces  errors.  The 
sunset  decay  period  is  less  critical,  since  the  per¬ 
centage  deviation  usually  increases  at  nighttime. 
This  implies  that  a,  would  then  decrease,  allowing 
a  greater  range  for  p, ,  thereby  increasing  the  proba¬ 
bility  that  the  requirement  p,  >  (a,/ 2)  will  be  met. 

It  is  seen  that  time  cells  with  an  extent  of  several 
hours  are  operationally  useful  at  solar  maximum. 
At  solar  minimum  the  absolute  value  of  the  rms 
deviation  is  low;  therefore  the  need  for  correction 
is  again  less  critical.  For  example,  the  residual  error, 
using  a  1  -hr  time  cell  at  solar  maximum,  is  equivalent 
to  the  average  monthly  variance  at  solar  minimum, 
indicating  the  adaptive  procedure  succeeds  in  re¬ 
ducing  the  error  at  solar  maximum  to  the  lower 
impact  solar  minimum  levels. 

APPENDIX  a 

For  actual  system  use,  the  monthly  mean  iono¬ 
spheric  behavior  would  be  predicted  by  a  numerical 
model  derived  from  worldwide  climatology.  Any 
model  has  some  inherent  bias,  which  may  have 
little  consequence  on  rms  error  when  the  primary 
concern  is  predicting  the  mean  behavior;  but,  as 
shown  in  Figure  5(b),  this  bias  can  affect  the 
adaptive  procedure.  The  following  derivation  will 
provide  some  perspective  on  the  import  of  this  bias. 

The  rms  error  is  defined  as 

I 

x!~-  (o.-py  (ad 

N 


where  P,  is  the  prediction  of  the  adaptive  procedure. 
For  the  case  of  a  multiplicative  function 

0„ 

Fm~zrM'  (A2) 

where  M  is  the  mean  as  predicted  by  the  model 
and  the  subscripts  f0  and  t  denote  the  reference 
time  and  the  time  of  applied  correction,  respec¬ 
tively.  Substituted  in  (Al)  and  expanded  in  terms 
of  o,,  cr,o  and  p,,  then 

X,j  -  of  +  p,1  +  (o*  +  pf0) 

M, 

~  2  —  («,<*, fi,  +  P,0lO  (A3) 

Now  defining  the  model  bias  B  such  that 


B 


*.« 


W, 

P. 


and  substituting  B,o  and  Bt,  the  results  analogous 
to  (6)  are 


(A4) 


It  can  be  seen  immediately  that  if  B,  =  B,  that 
is,  if  the  model  bias  does  not  vary  between  the 
times  being  considered,  use  of  the  model  will  be 
as  effective  as  the  use  of  the  observed  mean  in 
the  adaptive  procedure.  However,  empirical  models 
are  biased  by  their  data  base,  implying  preferred 
times  when  the  bias  is  less  variable,  and  therefore 
of  less  importance. 
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WTOXUCTiaN 

Over  the  past  decade  there  has  beoi  a  synthesis  in  the  understanding  of  nighttime  processes  in  the 
equatorial  icnsphere  using  contributions  cron  such  varied  obeervational  techniques  aa  incoherent  radar 
backacatter,  radio  scintillations,  optical  missions  and  in-aibu  measurements  (Baeu  and  Kelly,  1979; 
Fejer  and  Kelly,  1980  and  refexenaertherein)  •  Cons  ids  rations  from  aaronaty,  plasma  physics  and 
numerical  simulations  (Qssakow  et  al,  1979)  illustrate,  to  a  larger  extant,  the  unity  of  the  diverse 
obeecvations.  utilizing  the  fact  that  moat  of  the  contribution  to  TBC  ocmaa  from  the  vicinity  of  the 
F-regicnpeak,  TBC  observations  are  able  to  contribute  to  f1"— H  understanding  of  the  equatorial 
ionosphere. 

The  eyatnrial  icnoaphere  has  been  continuously  monitored  far  over  two  years  from  Ascension  Island 
(Ion.  15  W. ,  Lat.  8°  s.. )  through  measurements  of  total  electron  content  (TBC)  using  the  Faraday 
rotation  of  136  (tiz  signals  from  the  geostationary  satellite  SIRIO.  THo  characteristic  features  of 
the  nighttime  ionosphere  will  be  dl  aniseed  hare:  a  post-sunset  enhancement  which  lasts  for  several 
hours,  and  follows  a  rapid  sunset  decrease  ^attributable  to  the  usual  ionospheric  decay  processes; 
and  the  sudden  sharp  changes  in  the  electron  content  along  the  ray  path  that  are  often  superimposed 
on  these  enhancements  and,  to  a  lesser  degree,  on  the  subeequant  "normal"  nighttime  ionosphere. 

Similar  measurements  have  been  used  previously  in  the  equatorial  region.  Hunter  (1969)  in  Nairobi 
and  Koster  (1972)  in  Ghana  have  reported:  rapid  sunset  decreases  in  Faraday  rotation,  primarily 
emphasized  in  the  Ousia  data;  and  post-sunset  anhanoaiants  which  nay  equal  or  exceed  daytime  levels, 
greatly  emphasized  in  the  Nairobi  observations. 

The  sudden,  sharp  variations,  to  be  defined  here  aa  superimposed  "structures",  are  most  often  seen 
as  depletions  in  TEC.  similar  features  are  shown  by  Tsunoda  and  Ttwle  (1979)  to  be  associated  with 
radar  backacatter  plumes.  Kauahika  and  de  Mandonca  (1974)  and  Aarons  and  ttiitney  (1980)  have  shown 
that  sharp  variations  in  TBC  are  coincident  with  scintillation  patches.  Yeh  et  al.,  (1979)  tie 
these  structures  to  F-region  deletions  measured  by  in-si tu  probes  and  the  plaena  "bubbles"  of 
theoretical  studies.  The  variation  in  oocumoe  of  thse  structures,  both  seasonally  and  diumally, 
within  the  27  months  of  observation*  are  axmainad  here. 


OBSERVATIONS 

The  data  base  far  this  study  consists  of  aver  two  years  of  oontinous  measurements  of  TEC,  from 
September  22,  1978  through  Dsowbei  21,  1980.  The  ray  path  to  the  SHOO  satellite  from  Ascension  Island 
is  directly  north  at  an  elevation  angle  of  80°.  A  schematic  representation  of  the  geographic  and 
magnetic  relationships  is  ahovn  in  figure  1.  It  is  seen  that  field-aligned,  plam  density  structures 
in  tte  F-region  that  cross  this  ray  path  map  to  altitudes  of  600  to  1000  km.  above  the  magnetic 
equator.  The  observation  period  is  near  the  maximum  of  the  present  solar  cycle;  topside  density 
measurement*  in  the  evening  sector  over  the  Atlantic  Ooaan  near  the  maxinun  of  the  last  solar  cycle  at 
altitudes  greater  than  600  ton.  (Burks  et  al,  1979) ,  suggest  that  TBC  measurements  should  be  highly 
structured  at  Ascension  Island. 

Figures  2  3  contain  examples  of  nighttime  TEC  aa  a  function  of  T JT,  illustrating  daily  and 

seasonal  variability.  Data  from  figure  2  (3)'  acme  from  fivs  (six)  consecutive  days  near  the  June 
(Daomtber)  solstice.  local  tins  lags  nr  by  ana  hour.  Grand  sunset  occurs  at  about  1845TTT  (1915CT) . 
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Figure  1.  Schematic  of  Ascension  Island-  SHOO  ray  path  relationship  to  flAld  lines  over  the 
magnetic  equator. 


A  quiet,  well-behaved  ionosphere  is  ««p*>**-*a  during  Jins  at  Asoenaicn  Island  and  fowl  on  the  10,  11, 
and  13  of  June  1980.  Oh  acme  nights  (June  10  and  13)  the  decay  that  begins  during  the  sunset  period 
persists  as  a  gradual  decrease  over  the  entire  night  with  only  small  hour-bo-hour  variations.  Ch 
other  nights  (June  9,  11,  and  12)  then  is  a  post-awaat  increase  in  TIE  peaking  about  3  to  4  hours 
after  sunset.  Often  this  local  mnxlmMi  exceeds  daytime  TEC  values.  After  the  imxisun,  TEC  usually 
decays  to  its  expected  pre-dawn  minimum.  Evan  during  this  period,  however,  there  are  exanples 
of  sudden,  sharp  variations  in  TEC.  Depletions  of  '15%,  aocoRpanied  by  intense  scintillation,  appear 
on  nights  (JUne  9  and  12)  with  a  particularly  large  post-mnaet  enhancement.  These  structures  have 
distinct  boundaries  which  suggest  east-west  dlmenstions  of  300  to  500  km,  as  discussed  below. 

The  examples  of  figure  3  are  from  Dsoenfear  1980.  UhllJoe  June,  the  atypical  case  is  the  night 
(December  18)  of  no  super inposed  structure  or  scintillation  activity,  an  all  other  evenings  there 
are  post-sunset  enhancements  of  large  amplitude  and  2  to  3  hours  duration,  the  largest  occurring  on 
the  17th.  This  day  had  several  structures  in  TEC  which  appear  as  depletions,  and  a  large  finger  like 
enhancement.  If  this  enhancement  is  pert  of  the  background  the  depletions  are  £  30%.  These 
structures  appear  as  early  as  2100  UT  (December  13  and  17)  or  as  late  as  0100  UT  (December  16) .  The 
first  structure  tends  to  have  sharp  bmmdarlas  in  tine.  Structures  appearing  thereafter  are  less 
distinctly  bounded,  suggesting  either  a  decay  process  or  the  existence  of  two  or  more  structures 
long  the  ray  path. 

In  the  examples  presented  hart, TEC  could  be  measured  continuously  throughout  the  day  including  the 
rapidly  changing  evening  periods.  This  is  not  always  the  case.  On  many  evenings  in  September 
throuoh  March  period,  fluctuations  in  Faraday  rotation  occur  so  rapidly  that  the  polarization  angle, 
and  theurtVMC,  are  indeterminate.  This  occurs  moat  often  with  the  first  appearance  of  structure, 
infrequently  in  the  late  night  period. 

The  nocturnal  and  seasonal  variability  are  mamerizsd  in  the  histograms  of  figure  4.  Each 
histogram  is  a  plot  of  the  paromvtage  of  nights  in  which  structures  were  observed  within  each 
12 -minute  time  cell  during  the  monthly  period  beginning  with  the  22nd  day  of  each  month.  There  is 
a  sharp  increase  (decrease)  in  the  fxequ»cy  of  occurrence  near  the  September  (March)  equinox  with 
a  broad  maximal  (minimal)  between  equinoxes.  There  is  a  hint  of  solar  cycle  dependences  in  the 
increased  frequency  of  occurrence  far  the  September  1979  to  March  1980  period.  This  behavior  is 
similar  to  that  for  scintillation  activity  disc  used  by  Aarons  et  al  (1980)  for  the  0°  to  70° 
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Iheee  structures  fT*“*-  earliest  *  in  ths  "f  "*»■  to  Daaater  parlod  and  latest  in  Am.  B» 
serlUst  structures  an  also  those  of  shortest  duration.  Typically ,  thay  last  oa  has  tut  vary 
trm  Less  than  ana  ainute  to  aawaral  hours  in  duration. 

Tte  "nr-Himai  behavior  shows  a  distinct  saaaonal  oontrol.  A  pronounced  peaking  in  oocurzance 
during  aarly  evening  hours  oocura  at  tha  equinoxes,  but  la  obacarad  at  tha  aolatioaa.  This  is 
particularly  obvious  naar  tha  Dacantaar  aolatica  whidi  has  an  oocursnos  fraquan^  oayarahla  to  that 
imr  tha  — for  the  aarly  evening  hours.  Tha  aarly  mmning  behavior  bands  to  persist 
throughout  tha  ni£it,  however. 


DISCUSSION 

Than  are  basically  two  effects  being  considered?  first,  tha  large  post-sweet  enhancements; 
second,  the  sharp  structures,  usually  depletions,  associated  With  thaeu  The  first  has  been  observed 
by  Hunter (1969)  and  foster  (1972)  and  discussed  by  foster  and  Beer (1972)  and  foetar(1972) .  They 
™nr.ii«w.  that  the  F-region  dynamo  thoeryteiahbeth,  1971)  offers  tha  most  plausible  explanation  for  this 
effect.  Since  then  vertical  drift  measuremBite  at  Jicanarca  showing  Urge  enhancements  in  ugMard 
drift  near  aunset(Fajer  et  al  1979)  give  dredanoa  bo  Rishbsth'a  theory.  The  use  of  this  theory  by 
feelis  et  al(1974)  in  modelling  these  drifts  produced  the  effect  observed  at  Jicanarca. 

RLahbeth(1971,  1977)  suggests  that  the  neutral  winda  in  the  P-xegicn  produce  polarization  fields 
that  short-circuit  through  the  E- region  in  daytime,  but  build  at  sweet  through  the  mxJdert'drcp 
in  E-region  acnJuctivity  This  produoes  large  v^ward  drifts,  resulting  in  a  rapid  lifting  of  the 
F- layer  and  massive  transport  of  plane  along 'fiAId  lines.  Near  the  magnetic  equator  this  o^eare  as  a 
rapid  swsst  drf-Hn,.  in  TEC  as  noted  by  foster  (1972)  with  oBeervatians  at  dip  latitude  4<>.  This 
effect  is  lass  pronounoed  in  Hunter' a (196 9)  observations  at  13P  dip  latitude,  and  even  less  in  the 
Ascension  T«i«nd  observations  at  15°  dip  latitude.  The  subsequent  enhancement,  however,  is  more  pro¬ 
nounced  in  Hunter's  observations.  The  plasm  oust  be  transported  along  the  field  lines  corresponding 
to  the  height  above  the  magnetic  equator  to  which  tha  F-layer  is  lifted,  indicating  tha  layer  is 
to  altitudes  i  600  tan.  The  topside  density  measurements  reported  by  Young  et  al (1981) 
conplaneit  time*  observations.  DMSP  measurements  in  the  post-swaet(19.5  LT)  sector  at  840  km. 
showed  mnooth  density  depletions  in  the  imnadiata  vicinity  of  the  magnetic  equator  with  enhancement 
in  the  +  leP  to  20°  magnetic  latitude  range. 

The  second  feature  considered  hare  is  tha  sharp  structure,  often  depletions, in  TEC.  Observational 
studies  of  the  equatorial  P-region  suggest  that  these  structures  are  plana  bubbles  along  the  ray 
•path.  These  bubbles,  have  been  observed  in  the  bottcxneide  (folly  et  el,  1976) ,  topside  (Burke  et  el, 
1979)  end  near  the  peak  (McClure  et  al  1977)  of  the  F  layer.  Presently  accepted  theory  suggests  that 
these  begin  as  an  E  x  B  or  oollisional  Rayieigh^feylor  instability  in  the  bottomeide  of  the  F  layer. 
The  growth  of  the  initially  nail  perturbations  depends  on  (1)  the  altitude  of  the  initial 
perturbations,  (2)  the  strength  of  the  bottomeide  density  gradient,  (3)  the  altitude  of  the  F  peak,  and 
(4)  the  strength  of  the  poet- sunset,  eastward  elactrit  field  (Oeeakow  et  al,  1979  ;  Anderson  and 
Haerendel,  1979) .  Thus,  large  plasms  depletions  at  and  above  the  peak  of  the  F  layer  most  likely 
exist  than  the  F  layer  is  moving  vertically  toward. 

It  taw  noted  above  that  most  of  the  contribution*  to  TEC  at  Aaoansian  Island  oomee  from  flux  tubes 
which  cross  the  magnetic  equator-  at  altitudes  between  600  and  1000  km.  Yowg  et  al  (1981)  have  found 
that  in  the  Atlsntic-African  longitude  sector,  plane  bubbles  appear  near  the  magnetic  equator  at  840 
km.  on  more  than  80%  of  DMSP  passes  between  2000  and  2200  LT  and  never  before  2000LT.  This  oanicides 
with  the  local  time  at  which  the  first  structure  of  the  evening  appears  in  TEC  at  Ascension  Island 
during  tha  seasons  of  highest  occurrence  frequency. 

The  first  structure  of  the  evening  is  usually  well  defined,  and  as  noted  abovev  lasts  for  about 
one  hour.  This  can  be  viewed  as  a  structure  of  finite  longitudinal  extant  drifting  across  the 
Ascension  Island  ray  path.  It  has  been  established  (Rishbrth,  1971)  that  rwt-sunaet  polarization 
electric  fields  cause  the  ionosphere  to  drift  eastward  *t  a  faster  speed  than  oorotative.  A  typical 
value  for  this  drift  is  100  m/sec.  Thus  the  structures  appear  to  hove  longitudinal  widths  of  about 
300-400  tan. 

Structures  that  appear  later  in  the  evening  usually  have  less  well-defined  boundaries.  Ttoo  possible 
explanations  are  presented.  First:  Cceta  and  Kelley  (1978)  suggest  that  plasma  bubbles  grow  as 
large  amplitude  Rayleigh-  Taylor  waves.  Sharp  density  gradients  develop -at  the  boundaries  of  the 
bubbles  which  give  rise  to  mall  wavelength  drift  wave  ihztabilitiee.  The  dfift  waves  grow  at  the 
expense  of  tha  density  gradient  causing  cross-  field  line  diffusion,  leading  to  decaying  boundaries.' 
Second:  the  ray  path  passes  through  more  than  one  structure.  This  occurs  if  the  wind  speed  varies 
with  height.  This  would  produce  a  shearing  effect  on  the  drifting  structures,  allcwirv}  portions  of  two 
or  more  structures  to  cross  the  ray  path  simultaneously. 


seasonal' ‘variation  in  the  occurrence  of  these  structures  hoars  a  similarity  to  the  behavior  of 
the  vertical  drift  measurements  of  Fejer  et  al  (1979) .  The  moat  structure  is  observed  in  norths  for 
which  the  (as  reversal  ehhanoenant  of  vertical  drift  is  greatest.  This  —  ses  the  relationship 
of  these  structures  to  the  evening  enhancement  in  7 BC,  i.e.  the  enhancement  ae  a  pre  condition  for 

the  existence  of  structure-  A  clue  to  the  difference  in  nocturnal  behavior  with  season  may  be 
provided  by  McClure  et  al  (1977) .  They  observed  that  acrne  "bubbles’  drift  more  slowly  than  others, 
and  in  sane  cases,  they  move  with  the  velocity  of  the  background  plasma,  if  the  faster  moving  ones 
were  to  oocur  at  the  equinoxes,  the  slower,  at  the  solstice,  those  occurring  at  the  would 

be  longer  lived. 
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Abstract.  Ion  line  enhancement  detected  by  the  Arecibo  430  MHz  radar  is  interpreted 
to  be  a  result  of  nonlinear  interaction  of  weak  Langmuir  waves  via  ponderomotive 
force.  These  weak  Langmuir  waves  originate  in  the  scattering  of  the  O-modc  heater 
wave  and  of  the  excited  Langmuir  waves  from  thermal  ion  acoustic  waves.  Theoretical 
analysis  shows  that  those  enhanced  density  irregularities  associated  with  ion  lines  are 
forced  ion  acoustic  modes;  but  their  frequencies,  determined  by  the  beat  frequencies  of 
Langmuir  waves,  are  not  necessarily  their  characteristic  frequencies.  It  is  expected  that 
much  stronger  ion  line  enhancement  would  be  observed  if  back-scatter  radar  of  lower 
frequencies  (e.g.  SO  MHz)  was  employed. 


1.  Introduction 

It  has  been  observed  in  ionospheric  heating  experiments  (see  Radio  Science,  November 
1974)  that  the  density  fluctuations  of  the  ionosphere  can  be  greatly  enhanced  by  powerful 
radio  waves.  Figure  1  shows  a  typical  power  spectrum  detected  by  the  Arecibo  incoherent 
back-scatter  radar  from  the  ionospheric  F  region,  which  is  illuminated  by  a  powerful  hf 
O-mode  radio  wave  known  as  a  heater  wave.  This  power  spectrum  contains  two  com¬ 
ponents,  namely,  plasma  lines  and  ion  lines.  Plasma  lines  are  the  signals  with  upshiflted 
and  downshifted  frequency  slightly  less  than  the  heater  wave  frequency  (/h).  They  are  the 
pair  of  narrow  lines  located  at  430  MHz  ±/h  in  figure  1  (where  430  MHz  is  the  Arecibo 
radar  frequency).  Ion  lines  correspond  to  the  double-humped  frequency  spectrum 
straddling  the  radar  frequency.  The  pair  of  ion  lines  are  broadened  by  the  Landau 
damping  with  a  Doppler  spread  of  the  order  of  the  ion  acoustic  frequency. 

Existing  theories  (Fejer  and  Kuo  1973,  Perkins  et  al  1974)  show  that  the  O-mode 
heater  wave  can  generate  a  spectrum  of  unstable  Langmuir  waves  near  its  reflexion 
layer  in  the  ionosphere  via  parametric  decay  instability.  These  unstable  Langmuir  waves 
propagate  within  a  narrow  cone  centred  on  geomagnetic  field  lines.  Their  angles  of 
propagation  arc  estimated  to  be  less  than  20°  for  the  power  of  the  heater  wave  used  in 
the  experiments.  However,  the  angle  between  the  Arecibo  radar  beam  and  the  geomag¬ 
netic  field  line  in  the  modified  ionospheric  F  region  is  about  45°.  Hence,  the  weak  plasma 
lines  at  45°  observed  at  Arecibo  are  interpreted  to  result  from  the  scattering  of  the 
heater  wave  and  of  the  unstable  Langmuir  waves  from  thermal  ion  acoustic  waves  (Fejer 
and  Kuo  1973,  Perkins  et  al  1974).  But  much  stronger  plasma  lines  at  90°  observed  at 
Boulder,  Colorado,  have  been  suggested  to  be  caused  by  a  similar  process  from  intense 
field-aligned  density  irregularities,  rather  than  from  thermal  ion  acoustic  waves  (Fejer 
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Figure  1.  A  typical  power  spectrum  detected  by  the  Arecibo  430  MHz  monostatic 
radar  during  ionospheric  heating  experiments  (adapted  from  Showen  and  Kim  1978) 
showing  HF-induced  enhancement.  Plasma  lines  are  the  pair  of  narrow  lines  located  at 
430  MHz±/i  where /h,  the  heater  wave  frequency,  is  of  the  order  of  several  MHz.  Ion 
lines  at  thermal  and  enhanced  levels  correspond  to  the  broad  double-humped  portion 
straddling  the  radar  frequency  (Hagfors  and  Zanilutti  1973). 


1975).  According  to  the  theory  of  parametric  decay  instability,  a  spectrum  of  unstable 
ion  acoustic  waves  propagating  nearly  along  geomagnetic  field  lines  is  also  concurrently 
generated.  It  is  clear  that  the  enhanced  ion  lines  at  45°  detected  at  Arecibo  cannot  be 
these  ion  acoustic  waves.  We  suggest  that  this  ion  line  enhancement  arises  from  the  non¬ 
linear  interaction  of  Langmuir  waves  via  ponderomotive  force.  Thermal  interaction  of 
Langmuir  waves  is  discussed  by  Perkins  (1974)  and  by  Lee  and  Fejer  (1978),  who  interpret 
the  generation  of  intense  density  striations. 

This  paper  presents  a  qualitative  investigation  of  the  Arecibo  ion  line  enhancement. 
Section  2  describes  the  power  spectrum  of  induced  low-frequency  number  density 
fluctuations  associated  with  the  enhanced  ion  lines.  As  shown,  these  induced  fluctuations 
are  forced  ion  acoustic  modes  whose  frequencies  (determined  by  the  beat  frequencies  of 
Langmuir  waves)  are  not  necessarily  their  characteristic  frequencies.  Interpretation  and 
discussion  of  the  Arecibo  ion  line  enhancement  are  finally  presented  in  §3. 


2.  Induced  low-frequency  number  density  fluctuations 

2. 1.  Ponderomotive  force 

In  the  presence  of  two  high-frequency  waves  (e.g,  two  Langmuir  waves,  or  one  Langmuir 
wave  and  one  electromagnetic  wave),  electrons  experience  a  wave-like  nonlinear  Lorcntz 
force  whose  frequency  is  the  difference  between  these  two  high-frequency  waves.  If  both 
thermal  motion  of  electrons  and  geomagnetic  effect  are  neglected,  this  nonlinear  force 
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acting  on  an  electron  has  the  following  expression 

xexp  [i(o>i  —  <02)  /  —  i  (*i  —  k°)  •  r — i  (91  —  92)] j  (1) 

where  E,  to,  k,  <p,  e  and  m  are  wave  electric  field  amplitude,  wave  frequency,  wave  vector, 
wave  phase  angle,  electric  charge,  and  electron  mass  respectively.  This  nonlinear  Lorentz 
force  is  generally  called  the  pondcromotive  force. 

Equation  (1)  is  appropriate  for  the  case  concerned,  since  the  phase  velocities  of  those 
high-frequency  waves  are  much  greater  than  the  thermal  velocities  of  electrons  in  the 
ionosphere,  and  the  wave  frequencies  are  about  five  times  the  electron  cyclotron  fre¬ 
quency.  Furthermore,  since  the  wave  frequencies  of  high-frequency  waves  (the  heater 
wave  and  the  excited  Langmuir  waves)  arc  just  slightly  larger  than  the  electron  plasma 
frequency  tup,  equation  (1)  can  be  approximated  by  (Kuo  and  Fejer  1972) 

F=  Re  [“’2^/2*  (£  •*;)  (*2 -ki)  exp  [i(<«,- W2)t-i(kl-kt)*r-i(<pi-tp2)lJ  (2) 

The  ponderomotive  force  experienced  by  a  singly-charged  ion  is  F(m/A/)~  10  4  F, 
where  (m/Af)  is  the  mass  ratio  between  an  electron  and  an  ion.  Ions  are,  therefore,  con¬ 
sidered  to  be  immobile  in  the  electric  fields  of  high-frequency  waves.  However,  electron 
bunching  caused  by  F  results  in  a  charge-separation  electric  field  Ecs ■  Hence,  the  net 
force  acting  on  an  electron  is  (F-eEc»),  while  that  acting  on  an  ion  is  approximately 


2 . 2.  Dispersion  relation  of  density  perturbations 

The  corresponding  number  density  perturbations  in  electrons  and  in  ions  can  be  obtained 
from  the  linearisation  of  their  equations  of  continuity 


^(ene)  +  V»ne<Je  ^-£C8j=0 

,  (em)+  V*(m<7i£n)~0 
ot 


(3) 

(4) 


by  assuming  that  all  perturbations  have  a  space-time  dependence  of  the  same  form  as  F, 
i.e.  of  the  exp  [i(ou-4U2)f-i(iri-ifc2)«r]  type.  Thus,  equations  (3)  and  (4)  lead  to 


and 


(tui  -  (Oi) 


m=no 


o\ 

e 


(ki-ki) 
(un- utt) 


(5) 

(6) 


where  «n  is  the  unperturbed  plasma  density,  and  oK  and  cr4  are  the  electric  conductivities 
for  electrons  and  ions  respectively.  The  density  perturbations  zic  and  zji  are  approximately 
equal,  since  the  plasma  remains  neutral  for  low-frequency  perturbations  on  scales  much 
larger  than  the  Debye  length.  In  the  case  of  interest,  the  wavelength  of  the  low-frequency 
perturbations  picked  up  by  Arecibo  430  MHz  monostatic  radar  is  about  35  cm,  which  is 
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quite  large  as  compared  to  the  Debye  length  (~0.3  cm).  Therefore,  cither  equation  (5) 
or  equation  (6)  may  be  written  as 


„  F»(ki-k2)ae0l 

«12  =  M0  (7) 

e-  (aii  —  u)«)((Te  4-  <7i) 

which  represents  the  low-frcqucncy  number  density  perturbations  induced  by  the  non¬ 
linear  interaction  of  the  two  high-frequency  waves  via  ponderomoiivc  force. 

In  the  absence  of  external  disturbance,  i.e.  F=0,  the  dispersion  relation  of  these  low- 
frequency  modes  can  be  obtained  by  setting  the  denominator  of  (7)  to  be  zero,  namely 

Oe  +  oi=0.  (8) 


2.3.  Determination  of  oc  and  oi 

Because  the  wavelength  (~35cni)  of  induced  fluctuations  of  interest  is  much  shorter 
than  electron  collision  mean  free  path  (several  hundreds  of  meters  in  the  F  region),  the 
electric  conductivities  ae  and  m  are  determined  by  Vlasov  plasma  dynamics.  In  addition, 
the  wave  frequency  of  these  induced  fluctuations  ai(  =  aii  —  oj>)  satisfies  the  following 
inequalities: 

til  <^w~kCg  <  aii  Dp  <  aip 

where  Dj,  Cs,  an,  Qc,  and  a>p  are  the  ion  cyclotron  frequency,  ion  acoustic  velocity,  ion 
plasma  frequency,  electron  cyclotron  frequency,  and  electron  plasma  frequency  res¬ 
pectively;  fc=  |Ai  —  Aii|.  Ions  can,  therefore,  be  treated  as  unmagnetised,  and  electrons  as 
strongly  magnetised.  Nevertheless,  since  the  electron  cyclotron  radii  (a  few  centimetres 
in  the  F  region)  are  smaller  than  the  low-frequency  wavelength  (~35  cm),  electrons  can 
move  freely  along  the  geomagnetic  field  to  neutralise  ionic  perturbations  as  long  as  these 
low-frequency  perturbations  do  not  propagate  perpendicularly  across  geomagnetic  field 
lines.  As  mentioned  before,  ionic  fluctuations  at  45''  arc  detected  by  the  Arccibo  mono¬ 
static  radar.  Thus,  neglecting  the  geomagnetic  effect  on  electrons  as  well  is  still  a  goou 
approximation. 

Mathematically,  ion  dynamics  are  described  by 


9 

dt 


f\  +  v 


9 

dr 


E  • 


9 

dv 


/i=0 


where  fi  is  the  ion  distribution  function  in  the  presence  of  electrostatic  disturbance 
(E  =  Ecs)  due  to  electron  bunching  caused  by  pondcromotive  force:  f  is  represented  by  a 
linear  combination  of  perturbed  and  unperturbed  components,  namely,  f\  =fw  +/io-  The 
unperturbed  distribution  function /io  is  assumed  to  be  of  the  isotropic  Maxwellian  type 


It  can  be  shown  that  the  ion  conductivity  is  given  by 


,Ti 


oifi.  /  a/\'/2  f  ucxpi-Mu-nro 

k2hr  \2Ti )  J  ,  V-(co/k) 


(*>) 


where  L  indicates  that  integration  follows  the  Landau  prescription;  h\  is  the  ion  Debye 
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length.  If  oj  is  real,  equation  (9)  may  be  written  as 

<ji  =  (iaifo/Ac®Ai2)  (1  —  2  an  exp  (  — ar)  J*‘  dr  exp  ai  exp  ( — oti2)] 

(10) 

or 

at  =  (iwrolkW)  [  1  -  «,Z(  -  a,)]  (11) 

where  Z(-*i)  is  the  so-called  plasma  dispersion  function  (Fried  and  Conte  1961); 
a,=(o./*)(A//2r,)i/2. 

Electron  dynamics  are  represented  by  the  following  Vlasov  equation 

f  fe  +  f  *  fe-e  (E+vxBo)^  A=0 
at  or  m  ov 

where /« is  the  electron  distribution  function  under  electrostatic  disturbance  £=£CB-F/e 
associated  with  ponderomotive  force;  Bt>  is  uniform  geomagnetic  field.  fe  is  similarly 
composed  of  two  components,  namely,  /e=/ei+/eo.  The  unperturbed  distribution 
function  feo  is  also  assumed  to  be  Maxwellian,  as  given  by  equation  (12) 

/,o = (m/2wreF2  exp  (  - mi’2/27'e).  (1 2) 

The  electron  conductivity  can  then  be  shown  to  have  the  form  (Bernstein  1958) 

»-!z [' -* u. />-»(-■ 2st »>)] 

where  8  is  the  angle  between  the  wave  propagation  vector  k  and  the  geomagnetic  field  B0; 
re  is  the  electron  cyclotron  radius  defined  by 

re=(27V/Qc2m)1/2. 

For  the  parameters  of  present  interest  the  term  with  sin2  f)  sin2  0  in  equation  (13)  can  be 
neglected.  Then,  after  some  straightforward  mathematical  manipulation,  equation  (13) 
can  be  written  as 

-^(•^--P(— drexp(-r2)j  (14) 

or,  in  terms  of  the  plasma  dispersion  function,  equation  (14)  takes  the  form 

C7„  =  (i  eowlkVlv2)  {1  +  OeZ(«e)]  (1 5) 

where  ott  =  (wik)(mj2T,.)'i2  (cos  0)~l.  Since  «t.<l,  equation  (15)  can  be  well  approxi¬ 
mated  by 

ov~eF{)Uijk2het.  (16) 

This  final  form  of  o,.  shows  that  the  geomagnetic  effect  on  electrons  turns  out  to  be  neg¬ 
ligible  for  the  low-frequency  fluctuations  of  interest.  Recall  that  the  wavelength  (~35 
cm)  of  these  induced  fluctuations  is  larger  than  the  electron  cyclotron  radius  ( ~2  cm) 
by  one  order  of  magnitude  and  that  their  wave  propagation  angle  (~45°)  is  far  from 
being  perpendicular  to  geomagnetic  field  lines. 

Substituting  equations  (9)  and  (16)  into  equation  (8)  and  separating  the  real  and 
imaginary  parts  leads  to  the  following  two  equations 

Re  (at)  —  Cj/c 

Im  («,)  =  tr»/»  (C/lvft  k  exp  (-  re/r,). 


(17) 

(18) 
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Equation  (17)  indicates  that  these  low-frequency  plasma  fluctuations  are  essentially  ion 
acoustic  modes,  whose  Landau  damping  rate  is  given  by  equation  (18). 

However,  in  the  presence  of  external  disturbance,  i.e.  F ^  0,  w  is  real  and  determined 
as  the  difference  frequency  of  two  high-frequency  waves.  Substituting  equations  (2), 
(11)  and  (16)  into  equation  (7)  yields  the  specific  form  of  number  density  perturbations 
caused  by  nonlinear  interaction  of  two  high-frequency  waves  via  the  ponderomotive  force 

»u«,  r>-  -  (**.£•)  •,l:;<w7’)l  (19) 

27i  \ki  k-i)  {(re/ri)+[l-aiZ(-a,)]  !} 

where  k =ki-k»  and  oj(  =  un-  cot)  is  the  driving  frequency  of  the  ponderomotive  force. 
These  low-frequency  density  perturbations  are,  consequently,  forced  ion  acoustic  modes, 
which  are  not  necessarily  excited  at  their  characteristic  frequencies. 

In  terms  of  the  Fourier  transform  nu,(w,  k )  of  nn(t,  r),  i.e. 

nn(w,  k)  =  j  y  t^  m2(t,  r)  exp  [  -  i(wl  -  k  •  r)]  d t  dr 

equation  (19)  becomes 

»»(«.*)-  ~~Eklfpf  °  n°  (*)  *£)  yL'i-  (^  +  n  •  Zt-c,))-1^1  exp  (-i(9i-«P2)] 

(20) 

where  Kl  is  a  very  large  volume  and  t\.  a  very  long  time  interval. 

2.4.  Power  spectrum  and  radar  cross-sections  of  low-frequency  perturbations 

The  power  spectrum  <|/i(Jfc,  w>|2>  of  low-frequency  number  density  perturbations, 
caused  by  the  heater  wave  and  a  spectrum  of  randomly  phased  Langmuir  waves,  is 
obtained  from  the  generalisation  of  equation  (20)  as  follows: 

<  I »(»,  *>  I  *>  -  (z  «  |  -  v  (•W  (I  £  +  [I  -  «Z(  -  I  | 


(1^ 


i ki-*)' 


exp  [— i(?>i 


-9,)]|2> 


where  <  |  1 2>  indicates  that  the  absolute  value  of  physical  quantities  concerned  is  squared 
and  phased-averaged.  It  should  be  noted  that  k2 s  in  equation  (20)  have  been  replaced 
by  (ki— k)'s  in  equation  (21).  The  summation  in  equation  (21)  is  taken  over  all  pairs  of 
high-frequency  waves  whose  beat  products  satisfy  the  prescribed  value  (at,  k). 

If  those  randomly  phased  Langmuir  waves  are  characterised  by  a  spectral  density 
W(ki) 

«o Ehmn  =  2  W(ki,  km,  kn)(&k)3 

where  ki,  km  and  kn  are  the  three  perpendicular  components  of  k\  in  a  rectangular 
coordinate  system,  ki  =  /(Ak),  where  /=...,  -  2,  - 1 , 0, 1 , 2 . . . ,  and  AAr  is  an  infinitesi¬ 
mal  defined  by  27tFl_,/3.  It  follows  from  Lee  and  Fejer  (1978)  that  the  summation  in 
equation  (21)  is  replaced  by  integration  and  thus 

<|n(cu,  *)|2>=  16*4  VjJuiozTr2  {\(TelTi)+  [1  —  aiZ(  — ai)]_1|}"2 

x/J  W(ki)  W(ki-K)(dw!dkixY' d kiv  d*u  (22) 
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where  the  z  axis  is  taken  along  the  direction  of  geomagnetic  field,  and  the  double  integral 
is  taken  over  the  surface  w= constant. 

The  power  spectrum  A)|'->  is  related  to  the  average  differential  scattering 

cross-section  <p(<o0  +  to)  by 

o(u)o  +  ui)  dui=(2nVLti)-1  o'r(\n(w,  A)|2>  do>  (23) 

where  is  the  incoherent  radar  frequency  (430  MHz  in  the  present  case),  and  or  the 
Thomson  differential  cross-section.  Substituting  equation  (22)  into  (23)  leads  to  the 
final  result 

<t(w0  +  to)  do, r=  8  Tr>rto27V  2<7t{ \(Te/Ti)  +  [I  -  <x(Z{ - ai)]"1 1}"2 

*  [||  WikiWiki-k^duildkiz)--'  6kiy  dJtuj  dto  (24) 

which  specifically  expresses  the  radar  measurement  (per  unit  incident  power  per  unit 
solid  angle  per  unit  volume  and  per  unit  frequency  interval)  of  average  power  scattered 
incoherently  from  the  enhanced  plasma  density  fluctuations  due  to  nonlinear  wave-wave 
interaction  via  the  ponderomotivc  force. 

Computation  of  equation  (24)  seems  to  be  complicated,  but  a  crude  estimate  of 
ct(<o0+<o)  can  be  made  as  follows.  Remember  that  w  in  equation  (24)  is  the  frequency 
difference  between  two  Langmuir  waves:  one  with  the  wave  number  kt  and  the  other  with 
(ki-k),  namely 

to  =  <oi(*i)  -  m(ki  -  k). 

If  the  Bohm-Gross  relation  is  used  as  the  dispersion  relation  of  Langmuir  waves,  then 
<jti=\hez<o0(2ki  •k-ke) 

and  therefore 


du>ldkiT  =  'ihe'i<t>pk  sin  9  (25) 

where  6  is  the  angle  between  k  and  the  earth’s  magnetic  field  which  coincides  with  the  z 
axis,  k  lies  in  the  x-z  plane. 

It  follows  from  equation  (25)  that  (24)  can  be  written  as 

</(too  +  to)  =  8;rVto27'i  *3-1Ae"2ft>P  'A  '1  sin  }  0A(rxi)or 

*JJ  W{kx)W{kl-k)^kyyAkU  (26) 

where  i4(<xi)«{{(7V7i)+[l  -atiZ(-Qa)]-1  |}~z  determining  the  shape  of  ion  lines.  A(a,) 
versus  «i  for  Te/T\  =  1 .0,  1 .5  and  2.0  arc  plotted  in  figure  2.  For  a  larger  ratio  of  Te  to  7), 
/4(«i)  has  a  sharper  hump  located  at  a  higher  frequency.  This  indicates  that  Landau 
damping  broadens  the  ion  lines. 

It  is  interesting  to  note  that  o(ioo  +  u>)  in  equation  (26)  shows  a  simple  k~1  behaviour 
or  an  /  '  dependence  on  radar  frequency.  In  other  words,  the  theory  predicts  that 
radars  of  lower  frequencies  would  detect  stronger  ion  lines.  These  ion  lines  have  less 
Doppler  shifts  as  explained  by  an  example  given  in  §3.  Much  simplification  of  equation 
(26)  can  be  achieved  if  wc  assume  that  the  energy  density  »F(Ai)  of  Langmuir  waves  is 
constant  inside  a  cylinder  of  radius  kr  and  height  Ah  in  k  space.  This  means  that  the 
double  integral  in  equation  (26)  is  expressed  approximately  as 

|{  W{kx)W(ki-k)&kiytikuzz2W*ktk* 
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Figure  2.  Graphs  of  A  (ai)={  |(re/7i)+ [1  —  oti.Z(— ai)]~' 1[  !  as  a  function  of  ai  = 
(<o/*)(mi/27j)'«  for  r«/r,  =  ( A)  1.0,  (B)  1.5,  and  (C)  2.0. 


or 

JJ  W{k{)  W(k\-k)  dkiy  dku~2-*e0*E*kr3n^kn-' 

in  terms  of  the  amplitude  of  Langmuir  wave  field  E.  Then,  equation  (26)  can  be  written  as 


where 


o(u»o  +  oj)  =  tjE* 

Tj=4wMo2rr2eozfcr~3fctT,3~1/»e~2‘Op~1&  1  sin-1  0OT.d((Xi). 


Equation  (27)  is  a  convenient  form  of  the  theory,  from  which  one  may  estimate 
roughly  what  the  Langmuir  wave  amplitude  would  have  to  be  to  produce  the  observed 
radar  cross-sections. 


3.  Interpretation  and  discussion  of  Arecibo  ion  line  enchancement 

It  is  clear  from  figure  ]  that  the  ion  lines  observed  at  Arecibo  (Hagfors  and  Zamlutti  1973) 
are  modified  appreciably  after  ionospheric  heating.  However,  the  wave  number  of  these 
ion  density  fluctuations  is  1 8  m  ’,  which  is  larger  than  the  linear  dimension  of  the  A-space 
region  occupied  by  the  spectrum  of  unstable  Langmuir  waves  (see  figure  3).  This  spectrum 
of  unstable  Langmuir  waves  is  adapted  from  the  theoretical  calculation  of  Chen  and  Fcjer 
(1976).  According  to  the  theory,  these  unstable  Langmuir  waves  are  excited  by  an 
O-mode  heater  wave  near  its  reflection  layer  in  the  ionosphere  via  parametric  decay 
instability. 

Obviously,  the  ion  lines  detected  by  the  Arecibo  radar  at  45"  cannot  be  generated  as 
beat  products  of  these  unstable  Langmuir  waves.  They  may,  however,  be  so  generated 
by  a  spectrum  of  stable  Langmuir  waves  which  originate  in  the  scattering  of  the  heater 
wave  and  of  the  unstable  Langmuir  waves  from  thermal  ion  acoustic  waves.  These  stable 
Langmuir  waves  are  much  weaker  than  those  unstable  Langmuir  waves,  but  they  have  a 
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Figm  3.  Illustration  of  the  relative  magnitudes  of  the  ion  acoustic  mode  detected  by  the 
Arecibo  430  MHz  radar  and  of  the  saturated  spectrum  of  Langmuir  waves  in  the  k  space. 
The  saucer-shaped  regions  represent  the  latter,  which  is  adapted  from  the  theoretical 
calculation  of  Chen  and  Fejcr  (1976).  The  diagonal  arrow  represents  the  wave  vector 
of  the  ion  fluctuations.  Bo  is  the  local  geomagnetic  field  direction. 

broader  distribution  in  Ar-spacc.  Calculation  of  stable  Langmuir  waves  at  45°  has  quite 
successfully  reproduced  the  Arecibo  weak  plasma  line  echoes  (Fejer  and  Kuo  1973, 
Perkins  et  a!  1974).  We,  therefore,  interpret  the  weak  ion  line  enhancement  observed  at 
Arecibo  as  the  consequence  of  nonlinear  interaction  of  stable  Langmuir  waves  via 
ponderomotive  force.  Such  a  physical  process  is  formulated  in  detail  in  §2.  This  inter¬ 
pretation  leads  to  the  expectation  that  ion  density  fluctuations  of  longer  wavelengths 
(i.e.  smaller  wave  numbers)  will  be  enhanced  strongly  by  the  nonlinear  interaction  of 
intense  unstable  Langmuir  waves.  In  other  words,  it  is  predicted  that  back-scatter  radar 
of  lower  frequencies,  say,  50  MHz,  would  detect  at  Arecibo  ion  lines  which  are  strongly 
enhanced  by  intense  sources  (i.e.  unstable  Langmuir  waves)  in  addition  to  their  lc~l 
behaviour  shown  in  equation  (26). 

We  note,  from  figure  2,  that  the  ion  lines  are  located  approximately  at  |af|  =  ]  .45  in 
the  case  of  TcITt=*  1.5,  where  at  is  the  ratio  of  the  phase  velocity  (oj jk)  of  ion  acoustic 
waves  to  their  thermal  velocity  (27’t/Afi)’/2.  If  7)  =  1000  K,  the  frequency  of  forced  ion 
acoustic  modes  detected  by  the  Arecibo  430  MHz  radar  is  about  4.2  kHz,  which  agrees 
with  the  observation  (Hagfors  and  Zamlutti  1973).  Under  the  same  ionospheric  condi¬ 
tion,  ion  lines  detected  by  a  50  MHz  radar  are  estimated  to  have  a  frequency  of  0.5  kHz. 
That  is,  the  theory  expects  that  much  narrower  ion  lines  would  be  observed  by  a  50  MHz 
radar.  Moreover,  since  the  ion  lines  detected  by  a  50  MHz  radar  are  primarily  generated 
by  the  saturated  spectrum  of  parametrically  excited  Langmuir  waves,  the  Langmuir  wave 
field  £  in  equation  (27)  is  therefore  given  by  E = REo,  where  £o  is  the  threshold  field 
amplitude  of  the  decay  instability,  namely 

£o2  =  3.42  eo_I »o(7e  +  7i)  vel Wp 

and  R  is  the  ratio  of  £t2  to  £02,  where  Et  is  the  largest  field  amplitude  of  the  heater  wave  at 
the  standing  wave  maxima,  neglecting  absorption  and  the  magnetic  field,  that  is. 
Ex1 -(BPfj cqc)  (///Z),/2  (Fejer  and  Graham  1974). 

Currents  flowing  in  the  k  direction  of  ion  density  fluctuations  may  cause  the  two- 
humped  ionic  spectrum  to  become  asymmetrical  (Roscnbluth  and  Rostoker  1962).  It  is, 
however,  impossible  to  generate  such  strong  currents  through  ionospheric  heating.  Fejer 
(1977)  suggested  that  the  asymmetry  of  enhanced  ion  lines  might  be  caused  by  parametric 
interaction  of  the  incoherently  back-scattcrcd  waves  with  the  upgoing  probing  wave.  This 
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process  involves  a  parametric  instability  known  as  stimulated  Brillouin  scattering.  The 
back-scattered  waves  representing  the  lower  frequency  ion  line  are  amplified,  while  those 
representing  the  upper  frequency  ion  line  are  attenuated.  Fejer  (1977)  estimated  that  an 
asymmetry  of  ±  1  %  in  the  Arecibo  ion  lines  was  caused  by  stimulated  Brillouin  scattering. 

This  paper  has  presented  a  qualitative  investigation  of  ion  line  enhancement;  quanti¬ 
tative  analysis  of  it  requires  knowing  the  spectra  of  stable  and  unstable  Langmuir  waves. 
Experimental  spectra  of  Langmuir  waves  are  still  not  available,  in  spite  of  two  attempts  u 
obtain  them.  It  is  by  no  means  easy  to  calculate  the  theoretical  spectra  to  meet  the 
exact  experimental  conditions,  since  some  effects,  e.g.  the  self-focusing  of  the  heater  wave 
(Lee  1979),  are  now  qualitatively  clear,  but  it  is  difficult  to  evaluate  them  quantitatively. 
It  is  crucial  to  prove  the  existence  of  the  saturated  spectra  of  Langmuir  waves  now, 
because  interpretations  of  some  unexpected  phenomena  observed  in  ionospheric  heating 
experiments  are  based  on  their  existence.  Continuing  effort  to  measure  these  spectra  of 
Langmuir  waves  is  badly  needed. 
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which  leads  to  r So  criterion  of  exciting  tho  instability  by  an  i-mod#  heater 
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D  io*  #  o  os 

Tho  threshold  power  has  a  minimus  valuo  at  4  *4^  -  J*v^Jwq2| wn«r#  g  . 


Abstract  Jpfor  hybrid  aide-bands  and  a  ixjrely  (rowing  ar'-i*  an  bo  [>a r wsot n  - 
rally  excited  by  In  0-mod#  healer  wavs.  Thu  instability  can  .oopote  with  tha 
perimetric  dscay  instability  in  tho  loruispher ic  heating  sapor. ments.  It  may  bo 
responsible  for  tho  production  of  tho  plassM  lm*v  at  90°  at  Bouldor.  'o  lor  ado . 
I .  Introduction 

It  has  boon  generally  bolisvod  that  plasma  lino  enhanceeient  at  46°  dotsetod 


(aJ*<2kJv  *.1  t  \l  ^ )  l[  C  |2.  |C  I*  ■l1*11.  vhor#  a  ■  lo 
t  *  o  1  px  I  I  px I  th 

and  h  *  (m  n  1  •  (l/-*  Ik^v  *,()  *). 

•  i  t  # 


In  th*  cuss  of  Y -mods  I 


E  ■  -iE  Ths  threshold  powor  for 


tho  instability  can  bo  similarly  dorivod  as  IsE  'm  \ 

p*  0*1 

In  /»  *k^v  * ,  2tl  J)(v  ^  (w  2*{1  V*w  ^f^Hdw  'w  ~i7  l  »v  *(• 

si  t  •  *0*  O  ao*e 


J,v  ’.o  -a  4, 

os  Co* 


by  th#  Arecibo  410  *0ia  incoherent  radar  indicates  tho  excitation  of  parametric 
docay  instability  by  an  0-mod#  hoatsr  wav*  in  the  ionosphere  r#)#r  and  Kuo. 
19*11)  Pork  ms  ot  ol.,  1974)  Thoso  enhanced  plasma  linos  aro  interpreted  to  re- 
suit  from  tho  Acattennq  of  eh*  hooter  wave  and  of  th#  peramet r tool  ly  excited 
Langmcir  wives  off  thermal  ion  acoustic  wave*.  Hur  ouch  sort  in tans*  "field- 
alijned-  plasm*  lines  have  been  observed  *t  Boulder  tnd  era  characteristic  of 
extremely  narrow  spectral  widths.  radicating  .1  dj  f  ferent  mechanism  of  producing 
the#  (HinJtoff  and  Kreppel .  1976).  In  this  paper,  wo  present  a  theory  sugqest- 
inq  that  those  f leld-aligned  plasma  lines  bo  upper  hybni  modes  excited  para- 
metrical  ly  with  i  purely  qrowinq  ■rx'-e  by  an  O-mude  heater  wav*.  This  pora- 
metnr  instability  ...innot  bo  jxciteo  ,y  an  «stra<>rd mary  heater  wave,  howovar. 

1.  a  parametric  instability 

Near  th*  rwflotv.-tn  hoiqht  of  a  heater  wive,  the  pump  field  can  bo  approxi¬ 
mately  represented  ly  a  dipole  field  fp  •  5pe«pi  -  iw^t)  ,  whore  w0  it  th#  wave 
frequency,  and  ie  tho  *i#ld  atrenqth.  The  paramotrl*'  instability  under  in- 
vosii.iat ion  is  expected  to  produce  upper  hybrid  modes  and  purely  qrowinq  (i.o. , 
toro  frequency)  modes.  Since  tho  purely  qrowinq  eode  is  involved,  the  excita¬ 
tion  of  both  the  downshifted  and  upshlftsd  upf-er  hybrid  modes  hevo  to  bo  con- 
sulersd.  Tho  ;«rturbat ion*  in  physical  quantities  le.q.,  density)  duo  to  those 
modes  can  bo  represented  respectively  by  Sn^oxplyt  -ika),  Xn^espf -i (kx*wt)  I , 
and  dn^exp' Hhx-wti ) .  whore  w  •  it*iy,  end  v  is  tho  qrowth  rate.  A  roctanqular 
coordinate  system  has  boon  So  chosen  that  its  s  axis  coincides  with  tho  700-  ,• 
maqn#*lc  'le’d  1 B  I  and  it*  x  .1x1s  with  the  direction  ot  k. 

The  c"i*if  .eq  mode  -t-iuations  for  excitino  »  ..per  hybrid  wave-,  ny  scattorinq 
lh*  *  eater  v**'«i  .»ff  >h*  electrostatic  density  perturbation*  can  bo  derived  from 
*!ect.r''n  on». . nu.ty  equation  «nd  momentum  equoti  m .  Thoy  aro  given  by  '3^  • 


where  *■  the  electrostatic  (wren»*#l  of  eh*  downshifted  'upshlfted) 

upper  hybrid  waves,  Jn^  i»  th#  electron  density  «.erturt»ation  duo  to  th*  oxcitod 

aovo-fr*  ,uO«>vy  >m>dn.  ami  A  *  <w  •  1  •  »  l  -fi  C  1  /  <i« 

,  ,  o»*px#py  look 

where  >  •  w  *1  rv  *-w  ',  w  '  •  w  *  w  ,  0  ,  v  .  n  ,  and  v  aro  electron 

keeuk  o  tpeeo  t 

plasma  'r*  money,  -le:tron  ryr<.  frequency,  electron- ion  roliisiort  rmpioncy.  th# 
^ackqr  vund  ;  I  msj  J#^.>  *•/,  100  electron  velocity  respectively. 

The  >*h#r  ;on  relatinq  Jn  *0  0,  and  ‘i  .  'ibtamod  from  eVoctron  and 


rttuo  -tquati'.n*,  elect r'-n  heat  oqua* 


snd  riaet -  neutrality  approxima- 


*»ue,  can  oe  showr  •<>  ho  -  C'ay,  (2) 

wher#  •  •  •*. •ti4/l<‘)eiM)y||/f<J  and  C  •  12/1) 

i»V*  't?)  *  (w  *-  a  !i  )  (v  -iv  .«  ,’w  i  -iv  ik^k^)  .  whore  /  and  v  aro 
4  *  *  *  #  px  py  *  n  r»  3  pm  y-f 

the  porpenli.-  4j#f  -rngx-rent 4  of  the  induced  electron  v*ior-.ty  in  th*  presence 


of  the  tej’fr  wave  *ieid.  I  *1  end  k^  •  •*^1^  S  »ost itut  inq  ill 

into  fj)  leads  »3  tn*  .»  *por« ion  equation  n  •  CA  •  .  .c  I 


k“)  (w  •»*)  I  .w  *0  >/w  ■  <61 

Dos  040 

Note  that  5  is  defined  to  bo  w,  2*Q  2*v  ^ -w  (tear  th#  reflection  height  of 

k  a  *  o 

X -mode  heater  waves,  w  ■  (1/2)  (0  ♦((!  ‘*4w  "l*1).  Therefore.  4  ■  k^v  *v  "*11.2) 
^  o  a  *  p  t  « 

0  -  (3  .  2 \  (.7  ^*4w  *jV  In  th#  ionospheric  F  roqion:  T  -T  1000°k.  a  m  1.4* 

*  e  e  p  ei  oi 

10  s.  V  l  KHz ,  a  /2*  1  .4  mz .  W  /:v-6  HH*.  2*/k  -JxiO  V.  4  car.not  be  positive 

a  a  p  D 

•  «t  th*  criterion  of  the  instability  ((6))  can  raver  be  satisfied.  Hence,  the 
instability  of  concern  cannot  bo  excited  by  X-oodo  heater  wave*  ;n  th*  iono- 


sphoro.  In  contrast,  in  th*  caso  of  0-« 


4 .  Balov ant  roles  in  the  ionosphere 


1  heater  waves,  ic  can  ba  shown  that 


For  ‘(-2*A)«  10  m.  th*  minimum  threshold  power  ( f4> ‘  can  be  approximated 

*•  "’“‘O111!2  *"•>  th.  »ro«h  T-S<|  Ep,l/lep«Uln-l> 

The  threshold  field  etrenqth  for  excitinq  th*  modes  with  X  •  1  e  is  14  tnv'm. 

which  is  far  below  the  field  strength  'a  few  hundred*  of  mv/mj  of  th*  incident 

hooter  wav*  field,  and  the  corresponding  growth  rat*  (time)  is  •  ^0  sec*1 

(20  ms)  if  I E  | 'IE  1  is  taken  to  be  10.  It  it  interesting  tc  not#  that 
'  px’  *  pxmnn 

the  ihreihold  and  th*  growth  rates  of  these  modes  are  comparable  to  these  of 
th#  modes  excited  bv  th#  parametric  decay  instability.  In  other  words,  th# 
paraeettlc  instability  of  interest  can  competa  with  parametric  decay  »r stabi¬ 
lity  in  ionospheric  heating.  Thia  seem*  to  be  supported  by  observation* of  the 
"mini -over shoot”  in  the  measured  power  of  the  enhanced  plasma  lines  (Shower,  and 
KiNk  i 97 Aj  ,  which  occurs  in  a  few  tens  of  mili-seconds.  Much  lower  threshold 
field  strength  (-2  *v/m)  is  needed  to  excite  modes  with  >  *  5  »,  and  th#  growth 
rat#  (time)  is  2  sec  1  <o.S  sec).  The  excitation  of  modas  with  longer  wav#- 
langths  probably  contributa  to  tha  "main  ovarshoot"  reported  by  Showen  and  Kim 
(1979).  For  1  >  10  m.  JeE^/m^  t^m4n  which  is  independant  of  th# 

waveienqth  ind  r  )  «  «  10_i{  { l*-5. <»*102k (|  C  |3  |e  |J  - 1 ' ** - 1 ;  .  The  threshold 

px’  I  px  I  min 

field  for  exciting  ths  modes  of,  say,  50  s  11  -0.0  mv/m  and  the  growth  rate 
(time)  is  0.12  sec  ^  (7. 5  *#cl .  Our  proposed  instability  may  be  responsible 
for  th*  plasma  line  enhancement  at  90°  and  stay  site  contribute  to  the  produc¬ 
tion  of  f inld-aitqned  ionospheric  irregularities  with  scale  sites  ranging  from 
meters  to  hundreds  of  meters  in  ionospheric  heating  experiments. 
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the  coupling  mode  equation  for  the  purely  growing  node  end  the  coupled  aodr 
equation*  for  the  upshlfiad  and  dovnshitted  decay  arxlti.  Solving  theee  equa¬ 
tion*  will  Iced  to  the  dlepereion  relation  for  the  instability  of  interest 
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£#tract.  A  parametric  decey  Instability,  which  involvee  four  wave  interaction 
with  an  electroatatlc  pump,  la  propoead  ea  a  process  of  producing  both  the 
large-  and  the  ahort-acala  lrragularltlae  in  the  loooapherlc  heating  experi¬ 
ments.  It  can  explain  the  detection  of  discrete  rather  than  continuous  spectre 
of  large-ec*le  irregularities  In  the  experiments. 

1.  Introduction 

The  generation  of  f lald-allgned  Ionospheric  Irregular It  las  wee  one  of  the 
unexpected  modification  affects  due  to  Ionospheric  hasting.  The  field-aligned 
structure  with  scale  lengths  of  hundreds  of  meters  or  larger  la  associated  with 
the  occurrence  of  artificial  spread  F,  while  that  with  scale  lengths  of  meters 
or  shorter  is  responsible  for  the  creation  of  strong  HF  through  UHF  backscatter 


cions.  The  identification  of  possible  mechanisms  has  stimulated 


s  flurry  of  thaorstlcsl  Investigation.  Thermal  self -focusing  instability 
(Parkins  and  Valeo.  1974;  Cragln  and  Fsjcr.  1974)  has  been  generally  believed 
to  be  the  cause  of  those  large-scale  Irregularities.  However,  discrete  rather 
than  continuous  spectre  of  irregularities  as  predicted  by  the  theories  wars  de¬ 
tected  by  the  radar  beams  scanned  rapidly  across  tbs  modified  region  In  the 
ionosphere  (Duncan  and  Rehnk*,  1971).  The  causae  of  ahort-acala  irregularities 
seen  to  ^a  very  coalesced  and  many  possibilities  have  been  indicated  (Inhester 
st.  si.,  1901  and  references  therein).  In  this  psper.  we  lavestlgste  the  para¬ 
metric  decay  process  of  O-mode  heater  wsvee  In  the  ionosphere  end  show  that  It 
may  be  an  efficient  process  of  producing  flsld-sligosd  lrragularltlae  with  s 
bread  range  of  scales  from  motors  to  kllomotoro. 

2.  Source  of  pfp  weveo 

0-mode  heater  wove*  con  oxclto  parametrically  o  Langmuir  wove  and  an  lon- 
scoustlc  wavs  near  Its  reflection  height  In  the  ionosphere.  A  saturated  spec¬ 
trum  of  Langmuir  waves  Is  expected  theoretically  to  be  formed  vie  the  nonlinear 
len  Landau  damping.  The  occurrence  of  such  e  parametric  decay  process  has  been 
strongly  indicated  by  the  plsamo  line  enhancement  observed  at  Areclbo.  Thaos 
•netted  Len^ulr  wave#  have  already  been  considered  eo  the  potential  source  of 
••aerating  tho  short-scale  Irregularities  (Perkins,  1974;  Lae  and  Fejer ,  1970). 
Along  tho  same  lino,  we  Investigate  four-wave  other  then  the  three-wave  Inter¬ 
action  discussed  by  Perkins  (1974)  as  the  mechanism  of  causing  the  Induced 
Ionospheric  Irregularities  In  the  ionospheric  heating  experiments. 


Langmuir  waves  excited  by  O-mode  heater  waves  ere  used  as  the  pump  to 

initiate  another  per«otrlc  instability,  through  which  upehiftod  and  downshifted 

<*eay  bands  and  a  purely  growing  decoy  mode  ere  produced.  The  matching  wave- 

veetor  and  wave  frequency  relatione  for  this  Instability  ere  b^kj+k , 

w  where  the  subscript  "I"  a  tends  for  the  Langmuir  p«ap.  "O"  end 
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"2"  for  the  upshlfted  end  downshifted  decay  bende,  respectively,  and  "l"  for  th* 
puroly  growing  docey  mode.  The  basic  equation#  employed  in  our  analysis  of  the 
Frobla  include  electron  and  Ion  continuity  equations,  electron  end  ion  momenta* 
•Oust ions,  end  electron  hoot  equation,  for  simplicity,  o  rectangular  system  of 
ceerdlmetee  le  eo  chosen  that  all  the  wave  vectors  He  on  the  x-e  plena.  The 
Um^ktli  pump  1#  iieummd  to  be  elong  the  geomagnetic  field  line,  which  coincide* 
•1th  the  ■  axle.  It  is  a  tedious  but  at t eight forward  procedure  to  formulate 


pump  field,  grtwch  rate,  gyrof reqwenry,  electron- ion  collision  frequency,  inn- 
scoustic  velocity,  msss  of  the  charged  psrtlclec,  and  electron  therm*:  velocity 


3 .  fact  tat  ion  of  field-aligned  irregularlt fee 

Setting  y  -  0  In  (1>  glvee  the  threshold  field  ( | EJ 1  -  |fc  e,|): 

-  (3/0)CJ(kjek*)k]2k"2w2(w2-nJr2(2*(.1/mt)k3vJ(;“2)f|k5  2uj*k2) ‘MvV^*  )2 


•  e  el  J  t '  “e  e  e  1 

♦(ll/6)kj.|/nj)3,2(.t/»I>(.t/Qt)2kjv3r2.(.1/.t)k3v[3/n3X|E1  |2.'|El(tJ-l)  ■Sll) 
These  expressions  can  be  much  simplified  analytically  for  the  excitation  of 
short-  end  large-scale  Irregularities:  (1)  for  (l '21(0^*  lk2v2/{i2»>l ,  t.e., 
*3*10  m  If  the  following  ionospheric  parameters  arc  used:  u^/m^-2  9**104, 
vt~I0* (m/sec) .  Q  /2s~1.4  MHt.  The  threshold  field  and  the  growth  rate  can  be 
expressed  approximately  as  |  187l”2(v/m)-21 (mv/m)  and  T-l . 2*|J Ej2 ' |Ej  |^2 

«O.09)>l-1.3}l3^>4.4(aec*1)  for  Xj  •  3  m.  The  threshold  field  required  for 
exciting  the  short-scale  Irregularities  le  quite  small  compared  with  the  Lang¬ 
muir  field  (say,  0.7(v/m))j  (2)  for  (ll/t)(mt/m#)k2v|Q“2««l,  l.e.,  ij»»30  m 

Ik«>.  |t,  |th-6.2»10"2»,(»/.>'31(»*/»)  .ml  ,.J.l«10'2(ll*1.05,103i]2(|El|2/ 

Our  proposed  mechanism  pradlcts  that  in  the  case  of  exciting  large-scale 
irregularlt lea.  It**  threshold  fields  are  required  for  modes  with  smaller  seel# 
eltee.  Moreover ,  mode#  with  smaller  scale  alsee  have  larger  growth  rates.  This 
may  axplalm  th*  dattctlon  of  dlacrst*  ret  her  than  continuous  spectra  in  the 
Duncan  and  tabaka'*  axparinanca.  It  ahould  be  stressed  that  the  instability  of 
concern  can  generate  both  the  abort-  and  larga-acala  Irregular  It  lee  This 
Instability  may  contribute  additlvely  to  the  excitation  of  Ionospheric  Irregu¬ 
lar  it  lee  with  other  mechanism*  dlscueeed  In  th#  literature.  Since  either  0-  or 
X-mode  boater  wave*  can  eleo  generate  lerge-ecele  lrragularltlae  via  the  eelt- 
focuelng  instability,  ut  the  instability  discussed  In  this  psper  can  only  be 
excited  by  O-mods  heater  waves,  we  expect  that  more  Intense  large-scale  Irregu¬ 
larities  can  be  produced  by  0-mode  than  by  X-nod*  heater  waves.  This  predic¬ 
tion  cm  be  teoted  by,  for  example,  radio  Star  or  satellite  scintillation 
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The  conditions  for  the  upper  hybrid  modes  to  be  excited  parametrically  by  an  intense  electromagnetic  wave  near  its  cut- 
oil  have  been  examined.  They  show  that  upper  hybrid  side-bands  can  easily  be  produced  in  the  ionosphere  by  an  obliquely 
incident  wave  of  ordinary  but  not  extraordinary  polarization. 


Lee  1 1 1  has  shown  that  upper  hybrid  and  diffusion 
modes  can  be  excited  parametrically  by  an  extraordi¬ 
nary  pump  wave.  The  criterion  of  the  instability  is 
pccj0  >  il]  +  »>g,  where  w(),  S2C,  and  pc  arc  pump 
wave  frequency,  electron  gyro  frequency,  and  elec¬ 
tron-ion  collision  frequency,  respectively.  However, 
this  criterion  cannot  be  satisfied  in  cither  laboratory 
plasma  or  space  plasma.  We  show  in  this  letter  that 
upshifted  and  downshifted  upper  hybrid  modes  can 
be  excited  concurrently  with  a  purely  growing  mode 
by  an  obliquely  incident  ordinary  wave  transmitted 
into  the  ionosphere  from  the  ground.  The  established 
criterion  prohibits  the  excitation  of  these  plasma 
modes  by  an  extraordinary  wave. 

The  process  under  consideration  is  a  parametric  in¬ 
stability  which  occurs  near  the  cutoff  of  electromag¬ 
netic  waves  in  magnetized  plasmas.  A  monochronic 
dipole  pump  field,  K  =  *p  cxp(  iui((f),  can  be  rea¬ 
sonably  assumed,  let  the  parametrically  excited  pure¬ 
ly  growing  (i.c.,  /.cro  frequency)  mode  take  the  form 
of  exp(7f  -  i/br),  where  a  positive  value  of  y  rep¬ 
resents  the  growth  rate  of  the  instability.  The  wave 
vector  ( k )  of  this  mode  is  taken  to  be  along  the  x 
axis  of  a  rectangular  system  of  coordinates.  The  z  axis 
coincides  with  the  imposed  dc  magnetic  field  ( Bq ). 
Then,  the  excited  upshifted  and  downshifted  upper 


hybrid  modes  can  be  represented  by  exp[-i(ta  +  wf)] 
and  cxp[i(ikx  -  tor)]  respectively,  where  to  =  <o0  +  fy. 

The  mode-coupling  equation  for  exciting  the  upper 
hybrid  waves  by  scattering  the  pump  waves  off  the 
electrostatic  density  perturbations  can  be  obtained 
from  electron  continuity  equation  and  momentum 
equation.  It  is  given  by 


where  is  the  electrostatic  potential  of  the  upper 
hybrid  modes,  5ns  is  the  electron  density  perturba¬ 
tion  due  to  the  excited  zero-frequency  mode,  ut  is 
the  electron  thermal  velocity,  is  the  electron 
plasma  frequency,  and  =  xb/bx  +  yb/by.  In  ob¬ 
taining  (1 ),  the  approximation  of  uniform  medium 
has  been  made.  Substituting  <j> ,  =  0  j  exp[i(Aar  -  u>/)] , 
<^2  =  0  2  exP  I  K**  +  u>r)]  ■  Sfis  =  8ns  exp(yr  i kx). 
and  E p  =  *p  exp(  i u>0t)  into  equation  ( 1 )  leads  to 

$l  =A(6n*/n0),  =  zl(5ns/n0) ,  (2) 


where 
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=  <^pe  (^0  +  ~  *Vpy 

*  iw0O  +  l'e(2w0  "  wl*)’ 

where  a  =  u*k  +  fig  +  v*  ~  wj  and  to^  =  oupc 
+  3/fc2u2;  the  subscripts  l  and  2  designate  the  down¬ 
shifted  and  upshifted  upper  hybrid  modes  respec¬ 
tively.  y  <  w0,  vt  has  been  assumed  in  obtaining  (2). 

The  other  equation  relating  bns  to  <pt,  derived 
from  electron  and  ion  momentum  equations,  elec¬ 
tron  heat  equation,  and  quasi-neutrality  approxima¬ 
tion,  can  be  shown  to  be 

B(6njn0)  =  c$\  -  c*$2  ,  (3) 

where 

B  =  Cyfteft  i/ve  +  k2c2)(y  +  2vtmjm,  +  *2ufo/ft2) 

+  h*2cl , 

C  =  l(kielmi)[2i>eu0l(uo  -flc)l(V  -  iup>,fte/w0) 
-  'v(*2/*d)  * 

where  upjt  and  vpy  are  the  penpendicular  components 
of  the  induced  electron  velocity  in  the  presence  of 
pump  field;me,mj,cs,  ftp  and  2 rr/fcD  are  electron 
mass,  ion  mass,  ion-acoustic  velocity,  ion  gyrofrequen- 
cy,  and  Debye  length,  respectively.  In  the  case  of  our 
interest,  co0  >  fte  and  k2/kQ  <  1  are  assumed.  Sub¬ 
stituting  (2)  into  (3)  yields 

B  =  CA*  +  C*A  .  (4) 

For  the  oblique  incidence  of  ordinary  pump 
waves,  the  perpendicular  components  of  the  pump 
field  are  approximately  related  as  ep>,  =  -ie^ .  The 
induced  electron  velocity  is,  then,  given  by  vpx 
=  -ieefxl[mc(coQ  -  ftc)]  and  vpy  =  eepxl\mc( w0 
-  fte)J .  (4)  can  be,  therefore,  written  as 

B^a0kepjt/mc|2{:ia(wo  +  fte)  ^wgl/CwQ  ftt.) 

-(^/Ac^Kuofwo  +  fteW^e^  +  o2)-'  (5) 

after  the  substitution  of  these  expressions  into  it, 
where 

°0  =  3(OTc/™i)I*2ye/wO<wO  -  nc)l  • 

The  threshold  power  for  exciting  the  instability  is 
determined  by  setting  y  =  0  in  (5),  namely. 


leepx/mclfh  -  3u2(mc/fn,  +  k2u2/2il 2) 

X  ltO0(w0  ftc)/u>pu](o2  +^2U)q) 

X  {2[a(coQ  +  ftc)  -  f2wu|/(co()  ftc) 
(*2/*£)u,0(a;0  +  ftc)}  1  •  (6) 

Since  the  RHS  of  (6)  must  be  positive,  it  requires  that 

o>  \  (k2 Iky  )coq(coq  -  ftc)  +•  +  ftc) .  (7) 

Inequality  (7)  is  the  condition  for  the  instability  of 
concern  to  be  excited  by  an  ordinary  pump  wave  in 
magnetized  plasmas  (c.g.,  the  ionosphere).  If  the  pump 
wave  frequency  is  large  compared  with  the  electron 
gyrofrcqucncy ,  (7)  can  be  simplified  to  be 

o>i’l  +  \(k1lkl)ul  .  (8) 

The  threshold  power  has  a  minimum  value  at  a  =  o0 
=  0  +  (02  +t'2o^)1/2,  where  0=^2  +  (<r2/*2))(w^/ 2). 
Solving  (5)  gives  the  growth  rate  of  the  instability.  It  is 

y  =  -avK  +  ize  la 2  +  (2k2v2vebl^l2) 

X(kpJ2/|epx|2h  Dl,/2.  (9) 

where 

fl  =  (mc/mi)  +  ^(*2uf/ft2) , 

b  ~  (mjmj)  +  5  (At2u2 /ft2 ) . 

The  typical  ionospheric  parameters  in  the  F  region 
arc:  r>c  «  1  kHz,  ft0/ 2tr  «=  1 .4  MHz,  1000  K, 

mjtriy «  3.4  X  lO"5,  cj,K./2n  «6  MHz,  2tr /kD  *  3 
X  10--*  m.  Let  us  investigate  the  possibility  of  exciting 
upper  hybrid  modes  with  the  wavelength  of,  say,  3  m 
in  the  ionosphere  by  an  HF  (say,  6  MHz)  ordinary 
pump  wave.  It  is  seen  that  o0  ~  and  the  condi¬ 
tion  for  the  instability  |eq.  (7)  or  (8)|  can  be  well 
satisfied.  The  required  minimum  field  intensity  can 
be  estimated  from  (6).  It  is  about  3.2  (mV/m).  The 
corresponding  growth  rate  is  1 .8  s  1  if 
kn* ■  l2/kpr Inijn  *  's  assumed.  This  threshold 

field  intensity  is  far  below  that  (a  few  hundreds  of 
mV/m)  employed  in  the  ionospheric  heating  experi¬ 
ments  performed  at,  for  example,  Boulder  and  Puerto 
Rico.  Therefore,  upper  hybrid  modes  can  be  easily 
produced  by  ordinary  pump  waves  in  the  ionosphere. 
The  comparison  between  the  predictions  and  the 
relevant  observations  in  the  experiments  will  be  pub¬ 
lished  elsewhere. 
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In  the  ease  of  the  cxtruordinuiy  pump  wavc,t- 
=  i e p v .  The  extraordinary  mode  tends  to  propagate 
along  the  geomagnetic  field  near  its  cutoff  [2}  and 
almost  becomes  a  right-hand  circularly  polarized 
wave.  The  induced  electron  velocity  can  be  shown  to 
be  t5px  =  -ieepjt/|wL,(to0  +  f2t)]  and  vpy 
=  -eepv/{/we{to()  +  C2c)J .  The  threshold  power  for  ex¬ 
citing  the  instability  by  an  extraordinary  pump  wave 
can  be  similarly  derived  as 

\eepxlmc\ih  =  §  [o2(w0  +  ^C)2I^\ 

X  (mjm{  +  *^/2i2|!){*'e<w0  +  fl2)2  +  co(2o2} 

X  ((iuq(u(|  S2c)  ti2(w0  +  S22) 

■  <*2/*d)(««  -  a5»‘_1  •  0°) 

F-'ot  the  instability  to  occur,  il  icquircs  that 
o  >  +  i22)/tOf)(u)Q  f2c) 

+  (Al“/Ap)(wq  +  S22  )(wq  +  S2C )/ cjq  .  (II) 

Near  the  cutoff  of  the  extraordinary  waves  in  the 
ionosphere, 

<Jh=fic/2+i<nc2+4t o2c)'/2. 

Note  that  a  =  w2*  +  S22  +  k2  Wy.  Therefore, 
o  =  2i22  +  *2u2  +n2  (Oc/2)(n2  +4w2c),/2  . 

wltich  can  never  be  positive  in  the  ionosphere.  Thus, 
we  may  conclude  that  upper  hybrid  modes  cannot  be 
excited  in  the  ionosphere  by  an  HF  pump  wave  of 
extraordinary  polarization  transmitted  from  the 
ground.  However,  the  accessibility  of  the  extraordi¬ 
nary  pump  wave  to  the  upper  hybrid  resonance  zone 
is  possible,  if  the  wave  is  incident  upon  the  plasmas 
from  cither  the  higher  density  side  or  the  higher 
magnetic  field  side.  The  typical  example  for  the 
former  case  is  the  incidence  of  the  wave  from  the 
topside  of  the  ionosphere.  The  latter  case  is  achiev¬ 
able  in  the  laboratory  (e.g.,  let.  [3 ] ),  the  inhomogene¬ 
ous  magnetic  field  can  be  generated  in  the  magnetic 
fusion  devices,  e.g.,  Tokamaks. 

In  addition  to  the  criterion  for  the  parametric  ex¬ 
citation  of  upper  hybrid  modes  in  the  ionosphere, 
there  arc  some  new  approaches  in  our  work,  which 


lead  to  different  results  in  comparison  with  the  previ¬ 
ous  relevant  studies  (e.g..  refs.  [4-  8]).  Since  the  ex¬ 
cited  low  frequency  modes  are  purely  growing  modes, 
four  wave  interaction  (i.e.,  pimp  wave,  up-shifted 
and  down-shifted  excited  high-frequency  inodes  puie- 
ly  growing  mode)  should  be  analyzed  as  described  in 
this  paper.  Except  Diinant’s  work  [4] .only  three- 
wave  interaction  was  considered  before.  Moreover,  a 
new  nonlinear  effect,  viz.,  a  nonoscillatory  beating 
current  is  included  in  our  analyses.  This  beating  cui- 
rent  at  zero  frequency  is  driven  by  pump  wave  in  the 
oscillatory  electron  density  perturbations  of  the  ex¬ 
cited  high-frequency  waves.  It  tends  to  compensate 
part  of  the  well-known  nonlinear  effects,  i.e.,  the 
ponderomotivc  force  and  the  differential  Ohmic  heat¬ 
ing.  This  new  nonlinear  effect  not  considerd  in  the 
previous  work  leads  to  a  marked  difference  between 
our  results  and,  e.g.,  Dimant’s  [4j  in  the  threshold 
fields  and  the  growth  rates  of  the  excited  instabilities. 
Finally  we  should  point  out  that  the  present  work 
is  restricted  to  a  uniform  medium.  The  restriction  of 
our  results  as  applied  to  the  ionospheric  heating  ex¬ 
periments  will  be  discussed  elsewhere. 

Wc  thank  the  two  referees'  comments  which  draw 
our  attention  to  the  interesting  work  in  the  previous 
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by  NSF  under  grant  numbei  ATM-8008296  at 
Polytechnic  Institute  of  New  York  and  by  AFGL 
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The  Role  of  Parametric  Decay  Instabilities 
in  Generating  Ionospheric  Irregularities 
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We  show  that  purely  growing  instabilities  driven  hv  ihe  saturation  spectrum  of  parametric  decay 
instabilities  can  produce  a  broad  spectrum  of  ionospheric  it  regularities.  I  he  threshold  field  |/-;,h|  of  the 
instabilities  decreases  with  the  scale  lengths  X  of  the  ionospheric  irregularities  as  \H,h\  *  X  'in  the 
small-scale  range  (<15  m)  but  becomes  an  increasing  function  as  *  X  in  the  large-scale  range.  The 
minimum  threshold  field  (-2  mV/m)  occurs  at  x  21)  m  Ihe  growth  rate  y  of  the  instabilities 
maximizes  at  X  —  6  m  and  decreases  drastically  us  y  r  x  with  scale  lengths  larger  than  a  few 
kilometers.  The  excitation  of  kilometer-scale  irregularities  is  strictly  restricted  by  Ihe  instabilities 
themselves  and  by  the  spatial  inhomogeneity  of  the  medium  These  results  are  drawn  from  the 
analyses  offour-wave  interaction.  Ion-neutral  collisions  impose  no  net  effect  on  the  instabilities  w  hen 
the  excited  ionospheric  irregularities  have  a  field-aligned  nature. 


1.  Introduction 

Parametric  decay  instabilities  are  expected  to  be  excited 
by  O  mode  heater  waves  in  ionospheric  heating  experiments 
[Perkins  and  Kaw,  1971].  A  saturation  spectrum  of  Lang¬ 
muir  waves  within  a  small  cone  around  the  geomagnetic  held 
can  be  produced  by  these  instabilities  via  cascading  process¬ 
es  and  nonlinear  ion  Landau  damping  |e.g.,  Fejer  and  Kuo. 
1973 a\  Perkins  et  al..  1974].  Although  experimental  spectra 
of  such  Langmuir  waves  are  still  not  available,  it  has  been 
generally  believed  that  plasma  line  enhancement  observed  at 
Arecibo  [Carlson  et  al..  1972]  is  the  evidence  of  the  excita¬ 
tion  of  parametric  decay  instabilities.  Based  on  the  theoreti¬ 
cal  spectra  of  parametric  decay  instabilities,  plasma  line 
enhancement  at  Arecibo  can  be  successfully  reproduced 
[Fejer  and  Kuo.  19736;  Perkins  et  al..  1974], 

The  relevant  roles  of  parametric  decay  instabilities  in 
causing  other  unexpected  modification  effects  due  to  iono¬ 
spheric  heating  have  been  actively  explored,  such  as  6300- A 
airglow  enhancement  [Fejer  and  Graham.  1974;  Weinsiock. 
1974;  Nicholson.  1977],  the  generation  of  short-scale  field- 
aligned  ionospheric  irregularities  [Perkins,  1974;  Lee  and 
Fejer,  1978],  the  possible  formation  of  a  soliton  and  its 
collapse  [  Weatherall  el  al..  1982],  etc.  However,  parametric 
decay  instabilities  do  not  seem  to  be  the  only  cause  leading 
to  the  generation  of  short-scale  irregularities.  Their  domi¬ 
nant  role  is  questioned  in  particular  by  the  discovery  of  the 
overshoot  phenomenon  of  plasma  line  intensity  at  Arecibo 
[S/iowen  and  Kim,  1978],  Other  competitive  mechanisms 
have  been  suggested  with  continuing  effort  [e.g..  Das  and 
Fejer,  1979;  Inhesterel  al..  1981  and  references  therein;  Lee 
and  Kuo,  1982].  The  excitation  of  ionospheric  irregularities 
by  powerful  HF  radio  waves  is  probably  a  gross  effect 
contributed  additively  by  different  mechanisms. 

The  purpose  of  this  paper  is  to  investigate  the  generation 
of  ionospheric  irregularities  by  the  saturated  spectrum  of 
parametric  decay  instabilities  via  purely  growing  instabilities 
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( Perkins .  I974|.  In  his  study  of  three-wave  interaction, 
Perkins  1 1974]  has  already  commented  on  the  importance  of 
four-wave  interaction  for  exciting  purely  growing  modes.  In 
this  paper  we  analyze  a  process  of  four-wave  interaction 
including  beating  current  in  addition  to  differential  ohmic 
heating  and  nonlinear  l.orentz  forces  as  the  nonlinear  effects 
for  purely  growing  modes.  The  mechanism  for  the  genera¬ 
tion  of  field-aligned  ionospheric  irregularities  is  the  filamen- 
tafion  instability  of  Langmuir  waves.  The  beating  current  at 
zero  frequency  is  the  nonoscillatory  part  of  the  current 
driven  by  the  total  high-frequency  fields  on  the  total  high- 
frequency  density  perturbations.  It  appears  in  the  electron 
continuity  equation.  The  differential  ohmic  heating  force  is 
the  hydrodynamic  expansion  force  of  electron  gas  caused  by 
the  wave-induced  heating.  It  appears  as  an  additional  pres¬ 
sure  gradient  term  in  the  electron  momentum  equation.  I  he 
nonlinear  Lorentz  force  is  derived  from  the  convective  term 
of  the  electron  momentum  equation  and  is  reduced  to  the 
ponderomotive  force  in  the  unmagnetized  plasma  case.  It 
will  be  seen  that  the  differential  ohmic  heating  force  predom¬ 
inates  over  the  ponderomotive  force  in  the  consideration  of 
present  interest,  in  contrast  to  other  cases,  e  g.,  type  III 
solar  radio  bursts  [Bardwell  and  Goldman.  1976:  Weatherall 
el  al..  1981],  where  ponderomotive  force  is  the  dominant 
nonlinear  effect. 

The  plan  of  this  paper  is  as  follows.  The  excitation  of 
parametric  decay  instabilities  by  ()  mode  heater  waves  is 
briefly  described  in  section  2.  In  section  3.  purely  growing 
instabilities  driven  by  the  saturation  spectrum  of  parametric 
decay  instabilities  are  analyzed  in  spatially  uniform  plasmas. 
The  uniform  medium  theory  is  justified  by  the  result  that 
large-scale  (kilometers)  ionospheric  irregularities  cannot  be 
excited  favorably  by  the  purely  growing  instabilities.  The 
influence  of  spatial  inhomogeneity  on  the  generation  of 
ionospheric  irregularities  is  discussed  in  section  4.  and 
conclusions  are  finally  drawn. 

2.  Parametric  Decay  Insi  ami  ities 

An  O  mode  heater  wave  can  decay  into  one  l.angmiiu 
wave  and  one  ion  acoustic  wave  near  its  reflection  height  in 
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the  ionosphere.  After  (he  excited  Langmuir  wave  grows  and 
acquires  a  sufficiently  large  intensity,  it  can  decay  into  a 
daughter  Langmuir  wave  and  an  ion  acoustic  wave.  A 
cascade  to  lower-frequency  Langmuir  waves  proceeds  and 
finally  results  in  a  saturation  spectrum  of  Langmuir  waves 
[Pej-.-r  unJ  Kuo.  I973n;  Perkins  et  til. .  1974)  whose  distribu¬ 
tion  in  wave  vector  space  is  confined  in  a  narrow  cone 
around  the  axis  taken  along  the  geomagnetic  field. 

These  parametrically  excited  Langmuir  waves  are  now 
considered  as  pump  waves  to  excite  ionospheric  irregular¬ 
ities  (purely  growing  modes)  and  wave  vector  up-shifted  and 
down-shifted  daughter  Langmuir  waves.  The  process  under 
consideration  can  be  represented  by  the  following  matching 
relations: 


an  -  0  a>i  -  a>2  —  04i 

k.  =  k,  +  k,  k|  =  k}  -t  kn 

where  0.  1 ,  2,  and  3  are  used  to  represent  respectively  the 
down-shifted  daughter  Langmuir  wave,  the  Langmuir  pump 
wave,  the  up-shifted  daughter  Langmuir  wave,  and  the 
purely  growing  mode  (ionospheric  irregularities). 

To  begin  with  our  analyses,  the  Langmuir  pump  waves 
(wi.  k: )  are  represented  approximately  by  the  following 
dispersion  relation: 


(2) 


where  <u,,.  A;).  and  1),.  are  electron  plasma  frequency,  elec¬ 
tron  Debye  length,  and  electron  gyrofrequency  respectively. 
A  Cartesian  system  of  coordinates  is  chosen  with  its  z  axis 
along  the  geomagnetic  field,  and  the  wave  vectors  lie  on  the 
x-z  plane.  Since  no  damping  term  is  included  in  (2),  unatten¬ 
uated  Langmuir  pump  waves  have  been  assumed,  this  is 
because  the  Langmuir  pump  waves  with  saturated  intensities 
are  maintained  by  CW  O  mode  heater  waves.  For  simplicity, 
the  Langmuir  pump  waves  are  assumed  to  be  no!  depleted. 
In  fact,  the  process  discussed  in  section  3  can  modify  the 
saturation  spectrum  of  Langmuir  waves  determined  by,  for 
example,  cascading  processes  or  nonlinear  ion  Landau 
damping  and  c  <r  x  viewed  as  a  saturation  mechanism  of  the 
parametric  dac->  instability.  Therefore  the  effect  of  pump 
depletion  should,  .n  general,  be  considered  in  the  following 
analyses.  However,  this  is  beyond  the  scope  of  the  present 
work. 


3  Pur i-:i  v  Growing  Instabilities 
Excite i>  by  Langmuir  Pump  Waves 

The  coupled  mode  equations  for  daughter  Langmuir 
waves  (on.  k<i;  or-,  k>)  in  terms  of  Langmuir  pump  waves  (cui, 
ki)  and  ionospheric  irregularities  (m.  kd  can  be  derived 
from  the  electron  continuity  equation  and  electron  momen¬ 
tum  equation  as 


V 


n_  /  a 

Hi  \dr 


oW 


,  H 
ilr-  :- 
et 


fku_ 

+  I  —  I  H.  'oi/rt,' 

k,  i 


t  v,  j  ( k  i  •  k„) 


4>I 


(4.1 


where  4>i,  <t>;,  and  <J>o  are  electrical  potentials  of  the  Lang¬ 
muir  pump  wave,  up-shifted  daughter  Langmuir  wave,  and 
down-shifted  daughter  Langmuir  wave;  n0.  fin.  and  ly  are 
unperturbed  electron  density,  induced  nonoscillatory  elec¬ 
tron  density  perturbations,  and  electron-ion  collision  fre¬ 
quency  respectively:  the  superscript  asterisk  is  used  to 
indicate  complex  conjugate;  and  <uk/  =  o>p;(l  +  3 k,2lkp).j  = 
0,  2,  where  is  the  Debye  wave  number.  In  obtaining  (3) 
and  (4),  uniform  background  density  n0  has  been  assumed. 

The  other  coupled  mode  equation  can  be  obtained  from 
the  following  fluid  equations  for  electrons  and  ions: 
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—  n,.  +  V  ■ 
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n,V. 


=  0 


m\ST+ v'’vlv' 
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— - <  E  -  mll,.V,  x  i 
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(5) 


(b) 


H  ,  2 

-  t,.  +  ir«v  •  V,.  -  —  V  •  <7?  /  +  R , e  V ,  )7, 

61  ine 

-  2k  ~<Te  -  T»]  T  V'<">  (7) 

M 

H 

~  n,  +  V  •  h.V,  =  0  (8) 

Ht 

m(  -  +  v,  •  7  ]V,  =  -  —  +  fE  +  MSiy,  x  C 

\Ht  j  n, 

•t  vtm  —  (V,  -  Vj)  -  vinM\,  (9) 

n, 

where 

K  '  -  Mn{,Tdmi’r) 

P  '  -  jliv'il ,  ):(«(i7V/«f,  ) 

and  where  m  and  M  are  the  masses  of  electrons  and  ions;  n. 
V.  T.  11.  and  P  are  their  density,  velocity,  temperature, 
gyrofrequency,  and  thermal  pressure:  7»  is  (he  plasma 
temperature  (measured  by  electron  volts)  in  thermal  equilib¬ 
rium;  is  (he  ion-neutral  collision  frequency;  and  the  angle 
brackets  in  (7)  mean  the  time  average  taken  over  a  few  high- 
frequency  periods. 

According  to  the  matching  relations  given  by  ( I ).  waves 
(),'  ‘I,’  ‘2.’  and  \V  may  be  designated  respectively  by 

fip„  exp  (i(k„  •  r  mO]  bp i  exp  |i(ki  •  r  -  ..<n| 

bf>2  exp  |i(k;  •  r  <u/)j  fipxcxptyi  *  tki  •  r) 
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Fig.  1.  Configuration  of  four- wave  interaction.  Langmuir  pump 
wave  (k, >  propagating  along  the  geomagnetic  field  excites  a  down 
shifted  daughter  Langmuir  wave  (k,,)  and  a  purely  growing  mode 
(k3).  An  up-shifted  daughter  Langmuir  wave  <k2)  is  also  excited  as 
the  beat  product  of  the  Langmuir  pump  wave  (k,)  and  the  purely 
growing  mode  (k,). 

where  <u  =  w,  +  iy  and  a  positive  y  represents  the  growth 
rate  of  the  instabilities.  Equations  (5H9)  in  linearized  form 
can  be  written  as 


/  8ne\ 

yj  —  I  +  ik3  •  6V,.  +  iky  •  J  = 
\n0  ) 


=  ~{Tltlnu)ik\fm,  »  cK, 

-t  x  :  +  iy/u(8V,  -  6V,)  (9't 

where  P,„  =  y  +  t>m.  The  nonlinear  L.orentz  force  term  is  not 
included  in  (9'),  because  ions  are  too  massive  to  respond 
quickly  in  the  high-frequency  fields.  Eliminating  Ft.  6V,. 
6V,,  and  ST,,  from  (5'M9').  we  obtain  the  following  coupled 
mode  equation  for  the  purely  growing  mode: 

-IwJVi  +  — f  k'-  ^ 


1(5 


J  =  n0''(n,.Vf0*  +  riH>*Vr,  +  n,.,*\r2  n,.2\,,m) 

is  the  superposition  of  nonoscillatory  beating  currents  driven 
by  the  Langmuir  pump  wave  in  the  oscillatory  electron 
density  perturbations  of  the  excited  high-frequency  waves 
and  vice  versa;  n„  is  the  unperturbed  plasma  density. 


/F  =  -ik,STr  -  /k,r„ 


-  eEj  -  mil,6V,  x  j  -  mvr(S\r  -  6V,)  (6") 

where 

F  =  m(k,  •  Vh>*V„  -  ko  •  V„Vrt* 

+  k2-VH*Vf:  -  k.-V.yV,,*) 

is  the  superposition  of  nonlinear  Lorentz  force-  acting  on  an 
electron  which  oscillates  under  the  influence  w"  the  Lang¬ 
muir  pump  wave  and  experiences  a  force  exerted  by  the 
excited  high-frequency  waves  and  vice  versa.  This  force 
reduces  to  the  ponderomotive  force  in  the  unmagnetized 
plasma  case  and  includes  an  additional  term.  -(V,  x  V  x 
V,>,  in  the  magnetized  plasma  case.  E3  is  the  low-frequency 
electric  field.  The  inertial  term  has  been  neglected  in  (6'j 
since  y  <c  is  assumed  and  can  be  justified. 

8T,  =  <j <M  -  iTokj '  5V,.) 


v  ‘  ?  I’.im/M)  •  “  a  .'ft  -  •  *,,>  'I  '7| 

Ul„ 


where  //  =  2K.m(V,.|  •  V,.0*  +  ■  V,.,*)  is  the  differential 

ohmic  heating  term  arising  from  electron  energy  dissipation 
in  the  high-frequency  wave  fields. 

ySn,  +  in,Jk\  ’  6V,  -  0  (ft') 


bky,  11,  ky%x 

M  '  v,  M 


/  n,.u,  a,.-  , 
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+•  l  cik\x" - +  I  — —  II  h‘  ■+  — —  k\~k\-~  1  (k  •  J)  (10) 

*V  \  3  \  »y*  ''  / 


a  -  y  +  2 Mm/M)  +  (3/2/i„)(A;,..-«.  +  ku:R .") 

h  -  ky1  +  Aj,2(a,.n,/i',T,„) 

and  Fx.  F,.  and  F.  are  the  three  components  of  F. 

The  dispersion  relation  can  be  derived  from  (3).  (4).  and 
(10).  Instead  of  giving  its  complicated  general  expression,  we 
analyze  the  dispersion  relation  for  a  special  case,  wherein 
the  Langmuir  pump  wave  propagates  along  (he  geomagnetic 
field.  This  actually  represents  the  most  favorable  situation, 
because  this  Langmuir  wave  has  the  largest  intensity  in  the 
saturation  spectrum  of  parametric  decay  instabilities.  The 
study  of  this  case  is  therefore  sufficient  to  highlight  the 
problems.  Since  the  up-shifted  (k;)  and  down-shifted  (k,,) 
daughter  Langmuir  waves  are  in  symmetry  with  respect  to 
the  Langmuir  pump  wave  (see  (3)  and  (4)).  the  purely 
growing  mode  (kd  is  exactly  perpendicular  to  the  geomag¬ 
netic  field  in  this  particular  case.  Illustrated  in  Figure  I  is  the 
configuration  of  this  four-wave  interaction  process  in  wave 
vector  space: 

k i  —  k\Z  k3  -  Li 
k;  =  M  +  kx:  k„  -  -*vf  4  k 

Hence  (10)  reduces  to 
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a0  -  y  +  Ir.tmlM )  t 

where  r,.  is  the  electron  gyroradius.  With  the  aid  of  (2).  (3) 
and  (4)  become 

** = /  ~kxW{w  *  w) _  2<ji  v 
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<oi‘  --  «>,,-( I  +  3Af /A,r)  u,*  =  7Vm 

and  y  **.  k  has  been  assumed. 

In  the  presence  of  high-frequency  electrostatic  wave 
fields,  the  electron  density  perturbations  are  given  by 

Sn,  -  -A,'<l>/4ire  j  =  0.  1,2 

and  the  electron  velocity  perturbations  can  be  represented 
approximately  by 

<0,  11,  ,  .  A,,  1 

v,  =  -{elm)  x  — — ~v  A„  +  iy  — - Try  AM  +  z  <t>; 
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/V/  oil  (w,“  "  12,  ) 

Substituting  <  14)— <  17)  into  (II)  and  with  the  aid  of  (12)  and 
(13).  we  obtain  the  dispersion  relation 
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On  the  right-hand  side  of  ( I8«)  the  three  terms  labeled  by  h  . 
H,  and  J  correspond  respectively  to  the  three  nonlinear 
effects,  namely,  nonlinear  Lorentz  forces,  differential  ohmic 
heating,  and  beating  currents.  The  first  term  inside  the 
brackets  labeled  by  F  comes  from  the  ponderomotive  force, 
and  the  second  term  is  derived  from  -<V,,  x  V  x  V,.).  which 
reduces  to  zero  in  the  unmagnetized  plasma.  It  is  seen  that 
the  nonlinear  Lorentz  force  and  differential  ohmic  heating 
force  add  together  for  a*,  >  12,  and  are  partially  counterbal¬ 
anced  by  beating  currents  with  the  net  result  given  in  (18/;). 
Note  that  a»  =  y  +  lvr{mlM )  +  vekylrtr,  where  ( m/M )  -  3.4 
x  10  5  and  r,  -  2  cm  in  the  F  region.  Therefore  the 
contribution  of  the  ponderomotive  force  in  generating  field- 
aligned  ionospheric  irregularities  is  smaller  than  that  of 
differential  ohmic  heating  by  a  factor  of  Ayr,.2  for  irregular¬ 
ities  with  scale  lengths  (A,)  less  than  irr,.(2M!m)x  2  -  15  m 
and  by  a  factor  of  (m/M)  otherwise. 

If  the  following  ionospheric  parameters  are  used  for  the  F 
region  heating  experiments:  m/M( O' )  ~  3.4  x  10  \  iy  -  I 
kHz,  12,./2tt  ~  1.2  MHz.  v,  -  1.3  x  10' m/s.  co,/2jr  ~  6  MHz. 
K „  ~  0.5  Hz.  2w/A,  -  0.4  m.  and  £,  =  )A ,ct>,!  -  0.7  V/m, 
equation  H8/>)  can  be  simplified  as 


f  -y  — Ai2u,2  -  -  [2(m/M)(e/m)2|<l>,!2A|4Ax4(A|2  +  A,2)  '| 

M 
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Fig.  2.  Threshold  field  intensities  |£lh|  as  a  function  of  scale 
lengths  of  the  ionospheric  irregularities  excited  by  purely  growing 
instabilities. 
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It  is  interesting  to  note  that  since  :»  I ,  vin  is  finally 

canceled  in  (19).  The  reason  can  be  understood  from  (10) 
that  P,„  can  be  retained  in  the  expression  unless  the  purely 
growing  modes  are  exactly  field-aligned  (i.e. .  k\;  -  0).  In 
other  words,  i>,„  has  no  net  effect  on  field-aligned  purely 
growing  modes.  It  can  be  shown  that  the  retardation  on  the 
polarization  drift  of  field-aligned  modes  caused  by  damp¬ 
ing  is  compensated  by  an  £n  x  fl  drift  due  to  v,„  drag  force  in 
the  direction  perpendicular  to  the  plane  formed  by  the 
geomagnetic  field  and  the  polarization  drift. 

The  threshold  field  |£,h|  =  is  determined  from  (19) 
by  setting  y  -  0,  i.e., 
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Solving  (19)  for  y  gives  the  growth  rate 


y  =  ~*v(  r.  +  if  *3  VI  +  «v' 


M 
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+  v*iV 


(20) 


(21) 


The  results  of  (20)  and  (21)  are  plotted  in  Figures  2  and  3 
respectively.  Also  plotted  in  Figure  3  is  the  growth  time  of 


the  instabilities  defined  to  be  the  time  interval  for  the 
amplitude  of  excited  modes  to  exceed  their  thermal  level  by 
five  e  folds  of  magnitude. 

The  threshold  field  decreases  with  the  scale  lengths  of 
irregularities  in  the  small-scale  range  but  becomes  an  in¬ 
creasing  function  in  the  large-scale  range.  This  interesting 
behavior  of  threshold  field  versus  scale  lengths  can  be  seen 
from  (20)  where 

|£J  ~  (V3  (IHmU'HvrkMx)  *  M  2 

\  «  jrr,(2M//N)' 2  -  15  m 

j£"lh|  —  1 .2(mlM)v'v^mle)(vriu{k\lksfl,.)  y  A i  K  -*>  15  m 

While  |£(h|  has  a  minimum  value  of  1.89  mV/m  around  Xi  -■ 
20  m,  the  growth  rate  (time)  has  a  maximum  value  of  5f>.()  s  1 
(minimum  value  of  0.09  s)  near  \i  =  6m.  Although  meter- 
scale  irregularities  have  the  fastest  growth  rates,  larger- scale 
(tens  and  hundreds  of  meters)  irregularities  can  still  be 
excited  by  the  present  instabilities  in  reasonable  time  (a  few 
to  tens  of  seconds).  However,  the  thresholds  and  growth 
time  of  kilometer-scale  irregularities  increase  drastically. 
Consequently,  kilometer-scale  irregularities  cannot  be  excit¬ 
ed  by  the  present  instabilities  as  favorably  as  by  the  thermal 
self-focusing  instabilities  \ Perkins  and  Valeo,  1974,  Cronin 
and  Fejvr,  1974;  Berner  et  «/.,  1975). 

4.  Discussions  and  Conclusions 

Along  the  line  of  Perkins'  [  1974)  approach,  we  consider 
the  saturation  spectrum  of  parametric  decay  instabilities  as 
the  source  of  heating-induced  ionospheric  irregularities. 
Instead  of  three-wave  coupling,  we  analyze  four-wave  inter¬ 
action  whereby  the  Langmuir  pump  wave  can  excite  two 
daughter  Langmuir  waves  (wave  vector  up-shifted  and 
down-shifted  sidebands)  and  ionospheric  irregularities  (pure¬ 
ly  growing  modes).  This  is,  essentially,  the  filamentaiion 
instability  of  Langmuir  waves  in  a  magnetized  plasma. 
Moreover,  nonoscillatory  beating  currents  as  well  as  differ- 


Fig.  3.  Growth  rates  y  and  growth  times  7  as  a  function  of  the 
scale  lengths  of  the  excited  ionospheric  irregularities.  C  urves  A  and 
B  represent  respectively  the  growth  rates  and  growth  times  of  purely 
growing  instabilities.  Growth  time  is  defined  to  be  the  time  interval 
for  the  amplitude  of  excited  modes  to  exceed  their  thermal  level  bv 
five  e  folds  of  magnitude. 
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ential  ohmic  heating  and  nonlinear  Lorentz  force  are  the 
nonlinear  effects  in  the  formation  of  field-aligned  ionospheric 
irregularities.  Because  the  excited  low-frequency  modes  are 
purely  growing  modes,  the  beating  current  always  provides  a 
stabilization  effect  to  counterbalance  partially  the  differential 
ohmic  heating  and  nonlinear  Lorentz  forces.  The  principal 
results  of  our  analyses  are  summarized  as  follows.  (I)  A 
broad  spectrum  of  ionospheric  irregularities  can  be  produced 
with  scale  lengths  ranging  from  meters  to  hundreds  of  meters 
under  reasonable  experimental  conditions.  (2)  The  analyses 
of  four-wave  interaction  reflect  new  features  of  purely 
growing  instabilities  not  discovered  before  in  the  analyses  of 
three-wave  interaction.  The  threshold  field  |£,h|  of  the  insta¬ 
bilities  decreases  with  the  scale  lengths  X  of  the  ionospheric 
irregularities  as  |£(hj  *  X  2  in  the  small-scale  range  (<  1 5  m) 
but  becomes  an  increasing  function  as  |£,h[  *  X  in  the  large- 
scale  range.  The  threshold  field  has  a  minimum  value  of  ~2 
mV/m  at  X  ~  20  m.  The  growth  rate  y  of  the  instabilities 
maximizes  at  X  —  6  m  and  decreases  drastically  as  y  *  X  J 
with  scale  lengths  larger  than  a  few  kilometers. 

There  is  one  important  point  to  be  made  regarding  the 
maximum  scale  lengths  of  ionospheric  irregularities  excited 
by  electrostatic  pump  waves.  As  shown  in  Figures  2  and  3, 
the  instabilities  do  not  favor  the  excitation  of  large-scale 
(kilometers)  irregularities.  On  the  other  hand,  the  ionospher¬ 
ic  inhomogeneity  also  imposes  a  limitation  as  explained 
below.  To  generate  large-scale  irregularities  as  the  beat 
product  of  two  Langmuir  waves,  it  requires  that  the  coher¬ 
ence  of  these  two  short-wavelength  Langmuir  waves  be 
maintained  over  a  distance  of  many  wavelengths  in  the 
transverse  direction.  Under  quiet  and  spatially  uniform  (with 
typical  scale  height  of  a  few  tens  of  kilometers)  ionospheric 
conditions,  the  maximum  scale  lengths  of  the  excited  irregu¬ 
larities  are  determined  by  the  instabilities  themselves.  Other¬ 
wise,  they  are  primarily  determined  by  the  spatial  inhomo¬ 
geneity  of  the  medium. 

It  should  be  pointed  out  that  although  growth  rates  are 
determined  by  ionospheric  and  experimental  parameters,  the 
location  of  their  maximum  (around  Xi  -  A  m  in  Figure  3) 
does  not  change  significantly  with  different  parameters  (e.g., 
varying  j£||:  0.3-1. 2  V/m).  The  result  that  meter-scale  irreg¬ 
ularities  have  the  fastest  growth  rates  generally  agrees  with 
the  experiments.  However,  whether  the  growth  rates  of 
irregularities  maximize  at  the  meter-scale  lengths  as  predict¬ 
ed  by  this  theory  cannot  be  thoroughly  checked  from  the 
available  backscatter  radar  measurements  [e.g.,  Ruo  and 
Thome,  1974;  Minkoff,  I974|  and  awaits  future  experimental 
verification.  The  apparent  experimental  difficulty  is  that,  for 
instance,  irregularities  with  a  scale  length  of  10  m  can  be 
picked  up  by  a  backscatter  radar  of  15  MHz.  but  HF  radar 
echoes  cannot  be  simply  attributed  to  the  scattering  from 
irregularities,  because  refraction  is  important  too.  Lower- 
frcquency  HF  probing  radar  signals  cannot  even  reach  the 
heated  ionospheric  F  region. 

The  use  of  unattenuated  Langmuir  pump  waves  is  not  as 
unrealistic  as  it  appears  to  be.  As  mentioned  before,  they  arc 
waves  with  saturated  intensities  maintained  by  CW  heater 
waves.  As  soon  as  the  heater  waves  are  turned  off.  they  die 
out  within  milliseconds.  The  induced  ionospheric  irregular¬ 
ities  do  not  decay  quickly,  however.  The  reason  is  that  the 
field-aligned  nature  of  these  irregularities  prevents  ihcm 
from  being  removed  promptly  hy  the  diffusion  damping. 
Meter-scale  irregularities  can  persist  for  seconds,  and  large- 


scale  (hundreds  of  meters)  irregularities  can  even  survive  for 
hours. 

Although  our  attention  has  been  addressed  primarily  to 
the  generation  of  field-aligned  ionospheric  irregularities,  the 
filamentation  of  Langmuir  pump  waves  is  also  an  important 
result  in  our  work.  Schmidt  (1975)  shows  that  the  breakup  of 
solitons  is  just  a  continuation  of  the  filamentation  instability. 
A  comparison  between  the  present  work  and  the  recent  work 
on  a  soliton  and  its  collapse  [Weatherall  el  al.,  1982;  Hafizi 
et  al.,  19821  cannot  be  made  directly,  however.  This  is 
because  only  ponderomotive  force  is  included  to  analyze  the 
nonlinear  evolution  of  Langmuir  waves  via  the  oscillating 
two-stream  instability  \Weatherall  et  al.,  1982)  or  via  the 
scattering  of  Langmuir  waves  off  ion  acoustic  quasi-modes 
[Hafizi  et  al..  !982|.  In  our  future  work,  we  plan  to  take  a 
broadband  Langmuir  pump  in  the  analysis  of  the  modula- 
tional  instability.  Then,  the  nonlinear  evolution  of  this 
instability  can  be  studied  and  compared  with  the  results  of 
the  present  linear  stability  analysis.  As  mentioned  before, 
cascading  processes  and  nonlinear  ion  Landau  damping 
[ Fejer  and  Kuo,  1973a;  Perkins  et  al..  1974)  have  been 
suggested  as  the  saturation  mechanisms  of  the  parametric 
decay  instability.  The  process  discussed  in  our  paper  can 
also  be  considered  as  a  saturation  mechanism  of  the  para¬ 
metric  decay  instability.  It  occurs  on  a  much  slower  time¬ 
scale  than  the  cascading  processes  or  the  nonlinear  ion 
Landau  damping.  It,  then,  contributes  to  depleting  and 
broadening  the  saturation  spectrum  of  Langmuir  waves 
determined  by.  for  example,  the  cascading  processes.  There¬ 
fore  the  depletion  and  broadening  effects  of  the  present 
process  on  the  Langmuir  pump  waves  should  be  taken  into 
account  in  our  future  analyses.  If  pump  depletion  is  the 
saturation  mechanism  of  the  instability  proposed  here,  the 
intensity  of  the  purely  growing  modes  can  then  be  deter¬ 
mined. 

In  conclusion,  purely  growing  instabilities  driven  by  the 
saturation  spectrum  of  parametric  decay  instabilities  can 
excite  not  only  meter-scale  but  also  tens  and  even  hundreds 
of  meter-scale  irregularities  within  seconds  (see  Figure  3) 
under  optimum  (i.e.,  quiet  and  spatially  uniform)  ionospher¬ 
ic  conditions.  The  excitation  of  kilometer-scale  irregularities 
is  strictly  restricted  by  the  instabilities  themselves  (i.e..  their 
spatial  coherency)  and  by  the  spatial  inhomogeneity  of  the 
medium.  It  is  true  that  parametric  decay  instabilities  encoun¬ 
ter  competitive  processes  as  indicated  probably  by  the 
overshoot  phenomenon  of  plasma  line  intensity  [.S'/nm  en  and 
Kim,  1978).  However,  since  the  threshold  fields  of  purely 
growing  instabilities  are  a  few  millivolts  per  meter  in  the  tens 
of  meters  scale  range  (see  Figure  2),  parametric  decay 
instabilities  may  still  play  a  significant  role  in  generating 
intermediate-scale  ionospheric  irregularities  in  the  presence 
of  other  competitive  processes. 
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Faraday  Polarization  Fluctuations  of  Transionospheric  Propagation 

M.  C.  Lit 

Regis  College  Research  Center,  Weston,  Massac  husetts  021V) 

Rapid,  intense  Faraday  polarization  fluctuations  of  satellite  HF  and  VHF  beacon  signals  have  been 
shown  to  he  caused  by  wave  scattering  from  density  irregularities.  Except  in  the  case  of  HF  signals  (e  g..  20 
MHz),  irregularities  with  Gaussian  spectra  cannot  affect  the  polarization  of  satellite  signals  as  effectively  as 
irregularities  with  power  law  spectra.  Ii  is  adequate  to  use  the  theory  of  single  scattering  for  describing  the 
Faraduy  polarization  fluctuations  even  in  cases  involving  large  density  fluctuations,  such  as  those  of 
41-MHz  signals  near  the  magnetic  equator  and  those  of  136-MHz  signals  near  the  equatorial  anomaly 
region. 


Inikodiktiun 

Rapid,  intense  fluctuations  in  Faraday  rotation  angles  have 
been  observed  at  middle  and  high  latitudes  in  low  orbit  satellite 
beacon  signals  of  20,  40,  and  54  MHz  (Kent,  1959;  Par- 
thasaralhy  and  Reid ,  1959;  Yeh  and  Swenson.  1959;  McClure, 
1964;  Roger,  1965],  Similar  F'araday  polarization  fluctuations 
IFPF)  of  136-MFIz  signals  from  geostationary  satellite  have 
also  been  seen  at  locations  like  Ascension  Island,  Calcutta,  and 
Delhi,  situated  near  the  equatorial  anomaly  region  [ Klobuchar 
and  Aarons,  1980;  Das  Gupta  and  Muilra,  1980;  Lee  et  al„  1982 
and  references  therein].  In  the  observations  at  Ascension 
Island.  Ijce  el  ul.  [1982]  noted  the  coexistence  of  FPF  with 
strong  l.  band  amplitude  scintillations  as  displayed  in  Figure  I. 
These  two  phenomena  appeared  together  during  the  post¬ 
sunset  hours  and  disappeared  simultaneously  near  or  after 
midnight. 

F  PI  of  satellite  signals  have  been  interpreted  as  the  de¬ 
polarization  of  linearly  polarized  waves  by  wave  scattering 
from  ionospheric  irregularities  [Yeh  and  Liu.  1967;  Lee  et  al.. 
1982].  Yeh  and  Lin  ( 1967]  show  that  small  density  fluctuations 
will)  (i.itissian  spectra  are  sufficient  to  cause  FPF  of  2()-MHz. 
satellite  signals.  In  the  case  of  136-MHz  I  PF.  Ia-c  el  al.  [1982] 
show  that  small-scale  (  <200  ml  density  fluctuations  with  power 
law  other  than  Gaussian  spectra  are  responsible  for  this  effect. 
High  ambient  plasma  densities  have  been  shown  to  be  required 
to  cause  the  F  PF  of  1 36- MHz  signals. 

I  he  purposes  of  this  paper  are  to  investigate  the  different 
nature  of  ionospheric  irregularities  that  give  rise  to  the  FPF  of 
satellite  I  IF  and  VHF  beacon  signals  and  to  discuss  the  validity 
of  weak  scattering  theory  and  the  effect  of  multiple  scattering  in 
the  studv  of  F  PF  hi  the  ionosphere. 

I  AKAI1AY  Pill  AKI/ AVION  F I  IXI 1 IA  III  INS  I FPI) 

A  linearly  polarized  satellite  beacon  signal  can  be  decom¬ 
posed  into  an  ordinary  and  an  extraordinary  mode  for  its 
transionospheric  propagation  If  a  slab  of  uniform  density  ii- 
regiilanlies  is  assumed  to  exist  in  the  I  region,  each  mode  can 
interact  with  dcnsitv  irregularities  in  a  different  manner  as 
described  below  Assuming  Ilia)  tile  irregularity  scale  sizes  of 
interest  are  much  greater  Ilian  the  wavelength  of  the  signal,  one 
can  neglect  the  refraction  caused  by  spatial  inhomogeneily  and 
the  reflection  from  density  irregularities  For  simplicity,  one- 
may  only  consider  single  scattering  process,  which  introduces 
different  complex  shifts  into  the  original  ordinary  and  extra¬ 
ordinary  modes,  i.e.,  E„  t(z„)  -  F, ,  exp  (  i<f> , ).  where 
F.„ ,  |F„  )  refers  to  the  received  wave  held  of  the  ordinary 

Copyright  «'■  1982  by  the  American  Geophysical  Union 

Paper  number  I A 1606.  751 

0148 -0227/82/001  A- I606S01. 00 

The  U.S.  Qovarnmant  It  authoring  to  roproduct  and  tall  thlt  report. 

Parmlttlon  lor  further  reproduction  by  others  mutt  t>«  obtained  from 
•bo  copyright  ownor. 


(extraordinary)  mode  in  the  presence  of  density  irregularities, 
K,  ,(E, . )  is  the  wave  field  of  the  ordinary  (extraordinary)  mode 
in  the  absence  of  density  irregularities,  and  tf>tf<j>_)  the  phase 
shift  caused  by  wave  scattering  from  density  irregularities.  The 
phase  shifts  are  given  by 


<t>t  =  -tz0r,  |  J'r 


AN(r)  exp  [  - 


ikt(r  +  R  -  2„)] 
Rr 


where  rr(  -  2  82  x  10  15  m)  is  the  classical  electron  radius, 
A Nfr)  is  the  local  density  fluctuations,  k,(k_)  is  the  wave 
number  of  the  ordinary  (extraordinary)  mode,  V  is  the  total 
volume  of  density  irregularities,  and  as  illustrated  in  Figure  2, 
z„  is  the  distance  between  the  satellite  and  the  ground-based 
receiver,  r  that  between  the  satellite  and  the  scatterer  V,,  and  R 
that  between  the  scatterer  and  the  receiver,  respectively.  For 
quasi-longitudinal  propagation  in  the  ionosphere.  kt  =  lt[l 
4  X  Y  cos  11/21 1  -  X)]  where  X.  Y.  and  0  are  the  squared  ratio 
of  the  ambient  plasma  frequency  (/w)  to  the  wave  frequency  (/), 
the  ratio  of  electron  cyclotron  frequency  (12)  to  the  wave  fre¬ 
quency.  and  the  propagation  angle  of  the  signal  with  respect  to 
Ihe  geomagnetic  held;  k  f2nf.e\  (I  -  Xl1  2  is  the  wave 
number  for  propagating  in  an  isotropic  medium  (i.e..  V  -  0). 

An  induced  Faraday  rotation  angle  |Af2)  due  to  wave  scatter¬ 
ing  by  density  irregularities  is  just  the  half  of  the  differential 
phase  shift  between  the  ordinary  and  the  extraordinary  modes. 


All  =  \ftft ,  -  4> .  I  =  -  :„re 

f  .,  A/V(r|  exp  [-i  kfr  +  R-  :„)] 


sin  [rk(r  +  R  -  r0)] 


where  i:  =  X  Y  cos  0/2(1  -  A"). 

An  induced  wave  field  that  is  perpendicular  to  the  original 
one  can  be  visualized  to  be  associated  with  the  induced  Fara¬ 
day  rotation  angle.  The  ratio  of  the  induced  to  the  original 
wave  fields  is  approximately  equal  to  the  induced  Faraday 
rotation  angle,  if  the  differential  phase  shift  between  two  char¬ 
acteristics  modes  is  less  than  one  radian. 

I  he  physical  quantity  of  interest  is  the  variance  of  the  in¬ 
duced  Faraday  rotation  angle,  i.e.. 

otfi-vvUUy 

Ji  Ji  "i^GG 


cosffctr,  4  R, 


kfr 2  4  R  -  -  -ol] 


sin|f*(r,  f  R,  -  ;,i)|  sin  |ut(r,  4  R. 
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Fig.  I.  Sample  data  of  the  geostationary  satellite  signals  received  at  Ascension  Island  on  January  2,  1980  The  amplitude 
scintillations  of  the  136-MHz  signals  transmitted  from  SIRIO  and  those  of  the  /.  band  (I  541  GHz)  transmitted  from 
MARISAT-(1)  are  displayed  in  channels  (A)  and  (B).  respectively.  Shown  in  channels  (C)  and  (D)  are  the  Faraday  rotation 
measurements  of  136-MHz  signals  in  a  ramp  formal  with  one  channel  displaced  by  90  relative  to  the  other.  Rapid,  intense 
fluctuations  make  the  last  two  channels  almost  undistinguishablc  [Lee  cl  at.,  1982  ] 


where  the  angle  brackets  indicate  that  the  spatial  average  is  where  D  is  the  thickness  of  irregularities  layer.  *  the  zenith 
taken  over  the  volume  V ;  <AN(r,)AN*(r2)>  is  the  auto-  angle:  F'lr)  =  1 2C(2z  —  z0)  -F  4/fzz0  —  —  (C/2)21 1 } 

correlation  function  of  ionospheric  irregularities.  1=0,  ±e.  It  is  clear  in  (2)  that  if  F[±r.)  ±  h{0),  <|Afl|2> 

vanishes  and  no  FPF  can  be  observed.  Therefore,  FPF  requires 


Ionospheric  Irregularities 

Power  law  spectra  and  Gaussian  spectra  with  a  scale  size  of  1 
km  have  been  seen  to  be  associated  respectively  with  the  early 
and  the  late  developed  phases  of  irregularities  in  the  observ¬ 
ation  of  satellite  scintillations  [Bosu  et  ai,  1980u]  To  simplify 
the  problem  mathematically,  statistically  isotropic  density  ir¬ 
regularities  are  considered  and  assumed  to  be  generated  by 
spatially  homogeneous  processes,  namely,  the  correlation  func¬ 
tion  can  be  expressed  as  <AN(r1)AN*(r2)>  =  (AN)2  expf-(|r, 

-  where  (AN)2  is  the  mean  square  fluctuations  of 
the  ambient  plasma  density;  m  -  I  for  the  power  law  type  and 
m  =  2  for  the  Gaussian  type  of  ionospheric  irregularities 

Making  the  following  coordinate  transformation  in  (1)  x  - 
Jfx,  +  x2),  y  -  ify,  +  y2),  z  =  fcz,  +  z2),  {  =  x,  -  x2,  n  =  y, 

-  >2.  C  =  zi  -  *1,  and  integrating  ( I )  over  x  and  v  from  -  x. 
to  -F  oo,  one  obtains 


<  I 


Afi  |2)  =  jt(AN)2r,2  j*  dz  jjji/t  di\  dt 


expt-(CJ  +n1  +  i2)mlil "1 

•  (2F(0)  sin  |F(0M^  +  rj2)| 

-  He)  sin  \Hf  M^2  *  r?2)) 

•  /•<  >)  sin  !/•'<  -rMf’  *  n'H) 


S<0,0,0) 


f-  ijc  2  (ieometrv  of  iransionosphenc  propagation  and  wave  scatter¬ 
ing  hv  density  irregularities 


C) 
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that 

4e[zz0  -  z2  -  <//2)2]  2  2/<2z  -  z„)  (?) 

in  F(  ± z),  Cs  have  replaced  (*s  in  (3)  since  £  is  of  the  order  of  I. 

On  the  above  condition  it  can  be  shown  that 

<|An|l>  =  222r.2<(A  JV/N)j>,NjD 

•  sec  zCezJZo  -  zJ-l/(2nz0)]',J  (4) 

for  irregularities  with  power  law  spectra  [Lee  et  a/.,  1982], 
where  a  =  2n/k,  and  z„  is  the  distance  between  the  satellite  and 
the  mean  altitude  of  density  irregularities,  and  that 

<|An|2>  =  *-3Wr,2N2<(AN/N)2>0 

•  [ezjUo  -  =*l/co]20  sec  x  (S) 

for  irregularities  with  Gaussian  spectra  [  Yeh  and  Liu ,  1967]. 
Equations  (4)  and  (S)  can  be  used  to  estimate  the  required 
density  fluctuations  (<(AAf/Af)2>,/2)  for  causing  different  de¬ 
grees  of  FPF  (<|Afl|2>1/2)  under  the  given  ionospheric  condi¬ 
tions. 

Equaling  the  right-hand  sides  of  (4)  and  (S)  leads  to 

<(Aiy/iV)2)ul/2 
«A  N/N)2V'2 

=  2«i'2/j,V#  1,2[z0  'zjz„  -  zjtlc  (cos  fl)]  3,4 

0  ~/*2//2)J,VJ  (6) 

which  gives  the  ratio  of  Gaussian  to  power  law  types  of  density 
fluctuations  in  causing  the  same  level  of  FPF ;  where  /*,  O,  c. 
and  /  arc  the  ambient  plasma  frequency,  electron  cyclotron 
frequency,  speed  of  light  in  vacuum,  and  the  signal  frequency, 
respectively 

If  the  scale  sizes  of  irregularities  are  assumed  to  be  much  less 
than  the  slab  thickness  of  irregularities,  one  may  neglect  (f/2)2 
in  (3)  and  obtain 

2(/*,/>2in,:/*i  ' 

cos  (l  :J:„  :j  «;  H2zm  -  :„)  (7) 

on  the  substitution  of  i:  -  V  V  cos  0,  2(1  -  V).  X's  have  been 
replaced  in  (3)  by  r„'s.  which  represent  the  distance  between  the 
satellite  and  the  mean  altitude  of  density  irregularities.  Taking 
the  equality  in  |7).  one  can  determine  the  outer  scale  sizes  (/)  of 
irregularities  involved  in  causing  the  FPF  of  satellite  beacon 
signals  under  the  given  ionospheric  conditions 

FPF  of  Gfostahonaky  Satfj.mif  Hfacon  Sk.nai.s 
As  mentioned  before.  FPF  of  geostationary  satellite  beacons 
at  136  MH/  have  only  been  reported  from  ground-stations 
such  as  Ascension  Island.  Calcutta.  Delhi,  etc  .  near  the  cqua 
torial  anomaly  region  during  the  current  solar  maximum 
period  In  the  observations  made  at  Ascension  Island  (31  S 
dip),  very  high  /,',/  '  <16-2(1  MHz)  were  indicated  from  the  TI  C 
measurements  (  I  A  Klobuchar.  private  communication.  198 1 1 
and  from  the  airborne-ionosonde  measurements  of  ionization 
density  profiles  (f  J  Weber  and  J  <i  Moore,  private  com¬ 
munication.  198 1 1  While  large  FPF  of  low  orbit  satellite 
40-MHz  signals  have  been  seen  at  middle-latitude  stations,  e  g  . 
Cambridge,  J  ngland  |  Kent.  1459).  Faraday  rotation  measure¬ 
ments  of  geostationary  satellite  (ATS  6|  41-MHz  signals  made 
at  Huancayo  (1.3  N  dip),  Peru  in  1975  did  not  show  rapid, 
intense  fluctuations  (J.  A  Klobuchar.  private  communications. 


Fig  3a.  ((AN/V)2)1  2  versus  <|Afi|2>‘  2  for  power  law  type  of 

irregularities  to  cause  the  FPF  of  geostationary  satellite  41-MFIz  bea¬ 
cons  at  Huancayo  (curves  (A)  and  (B))  and  the  FPF  of  136-MHz 
beacons  at  Ascension  Island  (curve  (C))  and  at  Huancayo  (curve  (D)). 
The  propagation  angles  (0)  and  the  zenith  angles  (jr|  of  136-MHz 
signals  at  Ascension  Island  and  at  Huancayo  are  52°,  10  and  85  ,  70". 
respectively.  0  and  /  of  41-MHz  signals  at  Huancayo  are  82"  and  70  . 
The  ambient  plasma  densities  N  (or  frequencies/*)  at  Ascension  Island 
and  at  Huancayo  are  taken  to  be  4.05  x  10' 2  m  3  (or  18  MHz)  and 
1.25  x  10' 1  w  1  (or  10  MHz),  respectively.  But  in  the  solar  minimum 
case  at  Huancayo  (curve  |B)),  S  (or/*)  is  assumed  to  be  6.13  x  10" 
m  ’(or 7 MHz). 


1981 ).  The  fj-'i  at  Huancayo  were  generally  less  than  10  MHz  in 
this  solar  quiet  year. 

l.ct  us  assume  that  ionospheric  irregularities  have  a  thickness 
(l>)  of  200  km  at  a  mean  altitude  of  600  km(r0  -  zj.  Figure  3a 
shows  the  theoretically  calculated  <(A N!N)i')r'  1  versus  <  I  Afl 
|  ■’>'  ■’  (equation  (4))  for  power  law  type  of  irregularities  to  cause 
the  41-MFJz  FPF  at  Huancayo  and  the  136-MHz  FPF  at 
Ascension  Island  and  at  Huancayo;  <  |  AI2|2>'  2  =  0.1  radians 
corresponds  to  appreciable  FPF  and  <  |  AO|2)'  2  =  0.5  radians 
to  large  FPF.  The  outer  scale  sizes  of  irregularities  /  (equation 
(7))  involved  are  displayed  in  Figure  3b.  V  versus/*  (equation 
(6)1  arc  plotted  in  Figure  3c.  The  product  of  F  and 
<(A N/b/)2>J  2  gives  the  required  Gaussian  type  of  density  fluc¬ 
tuations  <(A/V/jV)2>,;  1  2  for  causing  the  same  level  of  FPF. 

The  136-MHz  FPF'  at  Ascension  Island  turns  out  to  be 
caused  by  small-scale  (<  200  ml  irregularities  with  about  I'L 
density  depletion  It  is,  however,  impossible  for  Gaussian  type 
of  irregularities  to  cause  the  136-MHz  FPF.  since  absurdly 


0  12  3  4 

tO  10  10  10  10 

I  (  m) 


I  ig  l/i  I  he  outer  stale  sizes  Id  of  irregularities  versus  ambienl 
plasma  Iicqueneies  I  /,l  in  causing  I  he  41-MHz  FPF  at  Huancayo 
(curve  it’ll  and  the  136-MHz  FPF  at  Huancayo  (curve  (A))  and  ai 
Ascension  I  aland  (curve  i  All 
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Fig.  3r.  T  versus  /„.  The  product  of  T  and  (lAN/N)1)  ful  gives  the 
required  Gaussian  type  of  density  fluctuations  <(AA//A/),)C1'2  for  caus¬ 
ing  the  same  level  of  41-MHz  FPF  at  Huancayo  (curve  (A))  and  for 
causing  that  of  136-MHz  FPF  at  Huancayo  (curve  (C))  and  at  Ascen¬ 
sion  Island  (curve  (B)). 

large  density  depletions  (e.g.,  <(AA/A/)2>1,2  =  834",,  for  <  |  Aft 
I2)1'2  =0.1)  are  required.  The  136-MHz  FPF  arc,  therefore, 
not  expected  during  the  late  phase  of  ionospheric  irregularities. 
Since  both  the  L  band  amplitude  scintilation  and  the  136-MH/ 
FPF  are  caused  by  small-scale  (  <  1  km)  density  irregularities, 
the  coexistence  of  these  two  phenomena  observed  at  Ascension 
Island  can  be  so  explained.  In  contrast,  the  observation  of 
136-MHz  FPF  at  Huancayo  is  not  expected.  This  is  because 
Huancayo  is  located  near  the  magnetic  equator,  where  the  f„F. 
is  relatively  low  and  the  propagation  angles  0  of  geostationary 
beacons  are  nearly  perpendicular  to  the  geomagnetic  field, 
therefore  as  shown  in  Figure  3.  meter  scale  irregularities  with 
unrealistically  large  fluctuations  (1.4%  for  <  |  AO  | 2 > 1  2  =0.1) 
are  required. 

The  theory  predicts  that  the  41-MHz  FPF  at  Huancayo,  if 
any.  are  caused  by  irregularities  with  an  outer  scale  of  several 
hundreds  of  meters.  This  leads  to  the  expectation  of  the  coexist¬ 
ence  of  41-MHz  FPF'  and  the  /.  band  amplitude  scintillation. 


F"ig.  4a.  Ionospheric  irregularities  at  middle-  or  high-latitudes  are 
assumed  to  have  a  thickness  of  50  km  at  a  mean  altitude  of  400  km. 
Curves  (A),  (B),|D),  and  (F)  correspond  to  the  FPF  of  low  orbit  satellite 
20.  40,  54.  and  136-MHz  signals,  respectively,  caused  by  power  law 
type  of  irregularities  within  plasma  density  background  of  1.25  x  10' 2 
m  5  (or/,  =  10  MHz).  Both  the  propagation  angles  (0)  and  the  zenith 
angles  (y)  are  taken  to  be  zeros  for  simplicity  But  irregularities  of  200 
km  thick  within  high  plasma  density  environment  of  4.05  x  10*  *  m  ’ 
(or/,  -  18  MHz.)  have  been  assumed  in  obtaining  curve  (Cl  for  the 
FPF  of  low  orbit  satellite  136-MHz  beacons  at  Ascension  Island. 
Curve  (F.)  shows  the  condition  for  the  F  PF  of  20  MHz  signals  caused 
by  irregularities  with  Gaussian  rather  than  power  law  spectra. 


Fig.  46.  Curves  (A),  (B),  (C),  and  (I))  show  the  outer-scale  sizes  of 
irregularities  involved  in  causing  the  FPF  of  136.  54.  40,  and  20-MHz 
signals,  respectively,  within  different  plasma  environments  ( /,). 


No  large  FPF  of  A'l  S  6  41-MHz  signals  observed  at  Huancayo 
in  1975  probably  can  be  attributed  to  low  /„F,  in  the  year  of 
low  sunspot  numbers  and  small  A N  near  the  magnetic  equator. 
Indeed,  the  1.  band  amplitude  scintillations,  which  arc  pri¬ 
marily  affected  by  A/V.  show  small  fluctuations  under  low  sun¬ 
spot  conditions  observed  at  Huancayo  fffu.su  el  id.,  19806],  It  is 
possible  to  sec  large  I  PF  of  geostationary  41-MHz  beacons 
at  Huancayo  during  the  current  solar  maximum  period.  How¬ 
ever.  if  A  A'  and  A  in  the  anomaly  region  are  greater  than  those 
at  the  magnetic  equator,  as  luroio  el  ul  f  1981]  have  discussed, 
the  F  PF  of  geostationary  41-MHz  beacons  recorded  at  Huan¬ 
cayo  will  not  be  as  large  as  those  of  geostationary  136-MHz 
beacons  measured  at  Ascension  Island 

FPF  ni  Low  Orbit  Saiiii.ih  Beacon  Sk.nai.n 

Since  the  FPF  of  geostationary  satellite  136-MHz  beacons 
have  been  observed  at  Ascension  Island,  it  is  natural  to  specu¬ 
late  that  the  FPF'  of  low  orbit  satellite  136-MHz  beacons  could 
also  be  seen  at  the  same  location.  As  shown  in  Figures  4a.  46. 
and  4c,  the  ionospheric  conditions  for  causing  the  latter  are  not 
very  different  from  those  for  the  former.  It  also  may  be  ex¬ 
plained  by  small-scale  irregularities  with  power-law  spectra. 

Kilometer-scale  irregularities  with  low  density  fluctuations 
are  sufficient  to  give  rise  to  the  FPF  of  low  orbit  satellite 
beacons  at  20,  40.  and  54  MHz  as  illustrated  in  F'igures  4a  r 
Fxccpt  in  the  case  of  20-MHz  s.gnals,  Gaussian  type  of  irregu¬ 
larities  cannot  effectively  cause  the  FPF  of  satellite  signals 

Imagining  that  four  polarimetcrs  at  Ascension  Island  are 
turned  on  to  record  low  orbit  satellite  beacons  of  20. 40.  54,  and 
136-MHz  separately,  one  expects  to  see  the  FPF  of  20-MHz 
signals  for  most  of  the  time  Once,  all  the  polarimetcrs  have 


F'lg  4,  L  versus  /,  Curves  (A I,  (B).  and  ((')  foi  20.  40.  and  S4.MH7 
signals  Hi  middle  or  high  latitudes  Curve  ID)  for  136-MHz  signals  -it 
Ascension  Island. 
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measured  the  FPF  of  20-  to  1.16- MM/  signals,  the  1 36- Mil/ 
FPF  fades  out  first  with  the  I.  hand  amplitude  scintillation;  the 
54-MH/.  FPF  and  then,  the40-MH/  FPF  die  out  later  on;  even 
in  the  late  phase  of  irregularities  characterized  hy  Gaussian 
spectra,  the  20-MH/.  FPF  can  still  persist  for  quite  a  while. 

Discussion 

Single  wave  scattering  from  ionospheric  irregularities  has 
been  used  to  study  the  FPF  of  satellite  HF  to  VHF'  beacon 
signals  in  this  paper.  It  shows  (hat  weak  kilometer-scale  density 
fluctuations  are  sufficient  to  cause  the  FPF  of  20-  to  54-MH/ 
signals  at  locations  not  near  the  magnetic  equator.  Its  use  is 
questionable  for  its  extension  to  cases  wherein  relatively  strong 
density  fluctuations  are  involved,  such  as  the  F  PF'  of  41-MHz 
near  the  magnetic  equator  and  that  of  1 .16- MHz  near  the 
anomaly  region.  The  validity  of  weak  (single)  scattering  theory 
as  applied  to  the  cases  concerned  is.  therefore,  briefly  discussed 
as  follows 

l  et  W'lOl  represent  the  power  of  satellite  beacon  signals 
scattered  per  unit  solid  angic.  per  unit  incident  power  density, 
and  per  unit  volume  by  irregularities;  d>  is  the  angle  between 
the  direction  of  incidence  and  the  direction  of  scattering.  If  the 
wavelength  is  much  shorter  than  the  irregularity  scale  sizes,  the 
scattering  has  a  sharply  directed  character,  namely,  most  of  the 
scattered  power  M(<J>|  is  restricted  within  a  small  solid  angle  <1*. 

W|«l>l  -  X/rr,.2/\(A/V|2>|  I  f  4n2/24 >2/a2]  2  (X) 

for  scattering  by  irregularities  with  power  law  spectra  |  Hooker 
and  Gordon.  I450|,  anil 

W  (4*1  -=  r. '  -7,  -’/  '  .  ( A,V|2>  exp  (  ir2(2<l*2  7.2)  (4) 

for  scattering  by  irregularities  with  Gaussian  spectra  f Fryer, 
|441] 

The  power  n  removed,  per  unit  length,  from  a  beam  of  unit 
power  density,  incident  at  a  zenith  angle  y  on  the  irregularity 
layer  (see  Figure  21  is  given  by 

<t  -  </<J>2n<l>W'(<t>)  sec  /  (10) 

After  the  incident  wave  traverses  the  layer  of  irregularities,  the 
power  density  of  the  unscatiercd  wave  becomes  /’  =  l‘„  exp 
(  -fH  n  dr)  ---  P„  exp  (-1)  where  /’„  is  the  incident  power 
density,  /)  the  layer  thickness  of  irregularities,  and  E  /(,’« r 

Weak  scattering  requires  that  the  unscattcred  power  density 
should  not  be  very  much  less  than  the  incident  power  density, 
namely.  £  has  to  be  small  in  comparison  with  I . 

E,  -  2/r.  2/2N2<(AN//V)2>0  sec  /  (II) 

from  (to  and  (10)  corresponding  to  power  law  type  of  irregu¬ 
larities;  and 

£,,  -  *'  2/r,2A2N2<(A/V//V>2>0  sec  /  (12) 

from  (4)  and  (10)  corresponding  to  Gaussian  type  of  irregu¬ 
larities.  -  0.1  if  one  substitutes  into  (II)  those  relevant  iono¬ 
spheric  parameters  required  for  causing  appreciable  FPF  (i.e., 
/|AQ|2>"2  —  0.1)  of  41-MHz  signals  at  Huancayo  as  de¬ 
scribed  in  Figure  !  Zr  -  0  27  corresponds  to  the  ionospheric 
conditions  lor  causing  appreciable  F  PF  of  1.16-MHz  signals  at 
Ascension  Island.  Weak  scattering  is,  therefore,  adequate  for 
the  description  of  FPF  even  in  the  presence  of  large  density 
fluctuations.  However,  since  E,  can  be  larger  than  1  in  cases 
involving  large  FPF  (e  g..  <  |  ADI2)1'2  >  0.5),  it  is  interesting  to 


estimate  the  cflecl  of  multiple  scattering  on  the  current  analy  ses 
based  on  weak  (single)  scattering  theory. 

Multi-scattered  waves  have  larger  angles  of  deviation  4>  than 
the  singly  scattered  wave  f Fejer.  145.1],  Consequently,  multiple 
scattering  can  result  in  larger  phase  shifts  of  ordinary  and 
extraordinary  modes.  Since  the  induced  Faraday  rotation  angle 
is  proportional  to  the  differential  phase  shift  between  the  two 
modes,  multiple  scattering  is  essentially  equivalent  to  the  effect 
of  increasing  the  density  fluctuations  <(A N/N)1)'  2  for  causing 
large  F  PF  via  single  scattering.  In  other  words,  (  |  Al)  j2> 12  x 
<(AN/A/)2>li2  is  not  true  for  large  FPF.  <(A N/N)2)1  2,  there¬ 
fore,  have  been  overestimated  for  large  FPF.  say,  <j 2 
-  0.5,  on  the  basis  of  single  scattering  theory.  The  over-esti¬ 
mation  is.  however,  not  serious,  since  there  is  just  a  factor  of 
5  between  the  <( AJV/.V )2 > 1  2  for  the  appreciable  FPF  (i.e., 
<  |Afi|2>‘  2  -  0.1)  and  that  for  large  FPF  (i.e.,  ( |Afl|2)'  2  = 
0.5)  No  qualitative  or  significant  quantitative  changes  arc  ex¬ 
pected  in  the  results  presented  in  this  paper  after  multiple 
scattering  is  taken  into  account. 
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The  nighttime  polarization  fluctuations  of  linearly  polarized  136  MHz  satellite  signals  received  at 
Ascension  Island,  located  near  the  southern  crest  of  the  equatorial  anomaly,  have  been  shown  to  be  the 
manifestation  of  depolarization  effect  due  to  the  diffractive  scattering  by  small-scale  (<200  m)  density 
irregularities  with  power  law  spectra.  The  theory  can  explain  its  coexistence  with  L  band  scintillation. 
The  absence  of  this  phenomenon  at  equatorial  locations  off  the  anomaly  crests  is  attributed  to  two  factors: 
(I)  the  ambient  plasma  densities  are  relatively  low,  and  (2)  the  propagation  angles  of  satellite  signals 
are  more  nearly  perpendicular  to  the  geomagnetic  field 


I  INTRODUCTION 

Intense  and  fast  fluctuations  in  Faraday  rotation 
angles  have  been  observed  from  low-orbit  satellite 
beacon  signals  transmitting  on  20  MHz  by  Parthas- 
arathy  and  Reid  [1959],  Yeh  and  Swenson  [1959], 
and  Roger  [1965],  on  40  MHz  by  Kent  [1959],  and 
on  54  MHz  by  McClure  [1964],  Kaushika  and  de 
Mendonca  [1974]  observed  that  the  Faraday  rotation 
angle  of  the  137-MHz  VHF  signals  from  a  geosta¬ 
tionary  satellite  sometimes  exhibited  fluctuations  of 
relatively  slower  periods,  ranging  from  a  few  seconds 
to  about  an  hour,  when  the  local  ionosonde  showed 
spread  F  echoes  at  the  low-altitude  station  of  Sao  Jose 
dos  Campos  (24°S  dip),  Brazil.  Recently,  Klobuchar 
and  Aarons  [1980]  and  Das  Gupta  and  Maitra  [1980] 
have  reported  intense  and  fast  fluctuations  of  the  Far¬ 
aday  rotation  angles,  in  association  with  strong  and 
fast  amplitude  scintillation,  respectively,  from  Ascen¬ 
sion  Island  (31°S  dip)  and  Calcutta  (32°N  dip),  both 
situated  near  the  crests  of  the  equatorial  anomaly. 
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Similar  observations  of  fast  fluctuations  of  Faraday 
rotation  angles  associated  with  strong,  fast  amplitude 
scintillation  have  been  observed  as  far  away  from  the 
equator  as  Delhi  (44°N  dip),  India  (Y.  V.  Somaya- 
julu,  private  communication.  1981 ),  but  they  have  not 
been  seen  at  Arequipa  (7.9°S  dip),  Peru,  located  near 
the  geomagnetic  equator,  nor  have  they  been  reported 
by  Yeh  et  al.  [  1979a,  b\,  in  observations  taken  at 
Natal  (7.9°S  dip),  Brazil. 

In  our  observations  from  Ascension  Island,  we 
have  noted  the  coexistence  of  Faraday  polarization 
fluctuations  with  strong  L  band  amplitude  scintilla¬ 
tion.  They  appear  together  during  the  postsunset 
hours  and  disappear  at  the  same  time  near  or  after 
midnight.  Saturated  136  MHz  amplitude  scintillations 
with  a  fast  fading  rate  are  always  seen  during  the  ob¬ 
servations  of  Faraday  polarization  fluctuations.  A 
sample  record  is  shown  in  Figure  1 .  In  an  experiment 
at  Ascension  Island  in  December  1979,  we  separately 
recorded  the  amplitude  of  the  left-hand  circular  and 
right-hand  circular  components  of  136-MHz  signals 
from  the  SIRIO  geostationary  satellite  observed  at  80° 
elevation  transmitting  nominally  linearly  polarized 
waves.  We  found  that,  during  the  times  of  strong,  fast 
amplitude  fading  when  polarization  fluctuations  oc¬ 
cur,  there  was  a  loss  of  correlation  of  the  opposite 
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Fig  I.  Sample  data  of  the  geostationary  satellite  signals  received  ai  Ascension  Island  on  January  2.  1980. 
Channels  (A)  and  (B)  show  the  amplitude  scintillation  of  the  136-MHz  signals  transmitted  from  SIRIO  and  that 
of  the  L  band  (1.541  GHz)  transmitted  from  MARISAT  I.  respectively.  Shown  in  channels  (C)  and  (D)  are  the 
Faraday  rotation  measurements  of  136-MHz  signals  in  a  ramp  format  with  one  channel  displaced  by  90°  relative 
to  the  other.  Intense  and  fast  polarization  fluctuations  make  the  last  two  channels  almost  undistinguishable. 


sense  circularly  polarized  wave  components  (i.e.,  or¬ 
dinary  and  extraordinary  modes),  indicating  that,  dur¬ 
ing  those  times,  the  concept  of  simple  Faraday  rota¬ 
tion  was  not  valid  at  our  observation  frequency.  It  is 
the  purpose  of  this  paper  to  explain  the  reasons  for 
those  observations  and  to  suggest  why  polarization 
fluctuations  of  VHF  radio  waves  have  not  been  ob¬ 
served  very  near  the  magnetic  equator. 

Ionospheric  density  irregularities  are  responsible 
for  the  irregular  fluctuations  in  the  Faraday  rotation 
angles  of  satellite  signals.  We  believe  that  the  polar¬ 
ization  fluctuations  as  shown  in  Figure  I  indicate  the 
depolarization  of  linearly  polarized  transionospheric 
satellite  signals.  Two  mechanisms  have  been  sug¬ 
gested  for  interpreting  the  polarization  fluctuations  of 
HF  satellite  signals,  specifically,  at  the  frequency  of 
20  MHz.  One  concerned  the  effect  of  refraction  im¬ 
posed  by  large  spatial  inhomogeneities  that  give  rise 
to  large  path  separation  of  the  ordinary  and  extraor¬ 
dinary  modes,  and,  consequently,  these  two  charac¬ 
teristic  components  of  a  linearly  polarized  signal  can 


scintillate  independently  [Roger,  1965).  The  other 
one  proposed  by  Yeh  and  Liu  |1967]  is  the  process 
of  wave  scattering  by  density  irregularities.  Induced 
fluctuations  in  Faraday  rotation  angles  due  to  this 
scattering  may  cause  depolarization  effect  on  the  lin¬ 
early  polarized  satellite  signals. 

In  Roger’s  mechanism,  it  is  the  ray  path  separation 
of  the  two  modes  but  not  the  density  irregularities  that 
is  primarily  investigated  for  the  polarization  fluctua¬ 
tions  of  satellite  signals.  The  magnitude  of  path  sep¬ 
aration,  which  is  large  enough  to  decorrelate  the  two 
modes,  is  estimated  roughly  as  A r  =  //((A4>)J)I/J, 
where  /  is  the  scale  size  of  density  irregularities,  and 
((A«t»)2)l/J  is  the  rms  phase  variation  of  the  signals. 
Roger  assumed  that  /  =  1  km,  <(A<|>)2>1/:  =  2  rad  for 
20-MHz  signals  and  obtained  A r  =  0.5  km.  The  cur¬ 
rent  in  situ  data,  however,  indicate  a  power  law  type 
of  irregularity  power  spectrum  encompassing  scale 
sizes  from  meters  to  tens  of  kilometers.  Moreover, 
Roger’s  mechanism  cannot  be  used  in  the  case  of  con¬ 
cern  here  to  explain  the  close  relationship  of  polari- 
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zation  fluctuations  with  L  band  scintillations  and  oth¬ 
ers. 

By  contrast,  no  significant  path  separation  of  or¬ 
dinary  and  extraordinary  modes  is  necessary  in  Yeh 
and  Liu’s  mechanism.  The  wave  scattering  by  density 
irregularities  is  considered  to  be  the  principal  process 
responsible  for  the  polarization  fluctuations  of  satel¬ 
lite  signals.  Nevertheless,  Yeh  and  Liu’s  theory  can¬ 
not  be  extended  to  the  case  of  VHF  signals  con¬ 
cerned,  as  explained  later.  As  a  matter  of  fact,  density 
irregularities  with  Gaussian  spectra,  as  assumed  in 
Yeh  and  Liu’s  theory,  tum  out  to  be  ineffective  in 
causing  the  depolarization  of  VHF  satellite  signals. 

Along  the  line  of  Yeh  and  Liu’s  approach,  we  show 
that  the  polarization  fluctuations  of  linearly  polarized 
136-MHz  satellite  signals  is  the  manifestation  of  de¬ 
polarization  effect  due  to  the  diffractive  scattering  by 
small-scale  (<200  m)  density  irregularities  with  power 
law  spectra.  The  absence  of  this  phenomenon  at  equa¬ 
torial  locations  off  the  anomaly  crests  is  attributed  to 
two  factors:  (1)  the  ambient  plasma  densities  are  rel¬ 
atively  low,  and  (2)  the  propagation  angles  of  satellite 
signals  are  more  nearly  perpendicular  to  the  geomag¬ 
netic  field.  It  can  be  shown  that  these  two  factors  lead 
to  the  requirement  of  rather  unrealistically  large  frac¬ 
tional  fluctuations  of  ionization  density  for  the  de¬ 
polarization  of  VHF  signals  to  be  seen  at  other  equa¬ 
torial  locations. 

This  paper  is  organized  as  follows.  The  theory  of 
single  wave  scattering  by  density  irregularities  is  pre¬ 
sented  in  section  2.  In  section  3,  the  variances  of 
phase  shift,  field  amplitude  ratio,  the  Faraday  rotation 
fluctuations,  and  the  induced  field  component  of 
transionospheric  signals  are  formulated  and  evaluated 
for  density  irregularities  of  both  power  law  and  Gaus¬ 
sian  types.  Interpretation  of  observations  and  discus¬ 
sion  are  given  in  section  4. 

2.  TRANSIONOSPHERIC  WAVE  SCATTERING 
BY  DENSITY  IRREGULARITIES 

For  simplicity,  density  irregularities  are  modeled 
as  a  uniform  mean  square  fractional  fluctuations  of 
ionization  density  <(AN/N)2)  within  a  slab  of  ioniza¬ 
tion  in  the  F  region  of  mean  density  N.  A  linearly 
polarized  wave  transmitted  from  a  geostationary  sat¬ 
ellite  5  is  received  on  the  ground  at  zenith  angle  x- 

The  center  of  the  coordinate  system  is  chosen  to  be 
located  at  the  geostationary  satellite  5  as  illustrated 
in  Figure  2.  In  Figure  2,  Z0  represents  the  distance 
between  the  satellite  and  the  receiver,  a  represents 


S  (0,0,0) 


A. 


Fig.  2.  Geometry  of  the  transionospheric  propagation  con¬ 
cerned.  The  center  of  the  coordinate  system  is  located  at  the  geo¬ 
stationary  satellite  S.  =  35.000  km.  a  =  34.300  km.  h  =  500 
km.  O  =  200  km.  9  =  52°.  and  x  =  10°. 

that  between  the  satellite  and  the  topside  of  density 
irregularities,  b  represents  that  between  the  receiver 
and  the  bottomside  of  density  irregularities,  D  rep¬ 
resents  the  thickness  of  the  slab,  and  0  represents  the 
propagation  angle.  The  unscattered  ray  path  within 
the  density  irregularities  is  D  sec  x  In  the  case  of 
concern  here,  Z„  =  35.000  km,  b  =  500  km,  D  = 
200  km,  a  =  34,300  km,  0  =  52°,  and  x  =  10°  are 
assumed. 

For  quasi-longitudinal  propagation  in  the  iono¬ 
sphere,  the  refractive  index  of  electromagnetic  waves 
can  be  approximated  as 

"I=*  1  ~/n<I  ±«/.>  (1) 

at  frequencies  large  compared  with  the  electron  gy- 
rofrequency,  where  fl  - /*//*,  i.e.,  the  squared  ratio 
of  the  ambient  plasma  frequency  to  the  wave  fre¬ 
quency,  fl  =  fl//,  i.e..  the  ratio  of  the  electron  gy- 
rofrequency  to  the  wave  frequency,  and  fl,  =  fl  cos 
0.  The  plus  and  minus  signs  in  (1 )  refer,  respectively, 
to  the  extraordinary  and  ordinary  waves  that  are  the 
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location.  That  is, 

E„(r)  =  (2£/r)l*  cos  (tkr)  +  y  sin  (tkr) J 


S 


R 


Fig.  3.  Geometry  of  the  wave  scattering  by  density  irregularities. 

two  characteristic  modes  of  a  linearly  polarized  wave 
in  the  ionosphere. 

Since  the  refractive  index  n  is  defined  as  kc/ta, 
expression  ( 1 )  can  be  written  as 

k„  =  t(l  +  e)  *,  *  *(1  -  e) 

as  the  wave  numbers  of  the  ordinary  and  extraor¬ 
dinary  waves,  respectively,  where  f.  defined  as 
-  f\)  has  a  small  value  (~  1(T4)  for  the 
VHF  propagation  in  the  ionospheric  F  region,  and  k 
=  (to/ c)(  1  ~  fl)'n  is  the  wave  number  in  an  isotropic 
medium  (i.e. .  e  =  0). 

In  the  far-field  approximation,  the  wave  fields  of 
these  two  modes  received  at  a  distance  of  r  from  the 
source  are  given,  respectively,  by 

Eto(r)  =  (E/r)(X  +  i$)  exp  |iU  -  kjr)\  (2) 

and 

Ejr)  =  (E/r)(x  -  iy )  exp  |i(u>/  -  k,r)]  (3) 

The  summation  of  (2)  and  (3)  yields  the  field  of  the 
linearly  polarized  spherical  wave  received  at  the  same 


•  exp  (i(u>t  -  Jkr)l  (4) 

showing  that  the  polarization  vector  of  the  linear¬ 
ly  polarized  wave  is  rotated  by  an  angle  of  tkr  - 
(1/2 ){k0  -  kx)r  known  as  the  Faraday  rotation. 

In  the  presence  of  density  irregularities,  namely 
Pn  =  (p»)  +  4/n  (where  (pN)  corresponds  to  the  mean 
ionization  density  and  A/*  corresponds  to  the  density 
irregularities),  an  additional  change  of  Faraday  rota¬ 
tion  angle  is  introduced  by  wave  scattering  from  den¬ 
sity  irregularities.  Since  satellite  signals  are  received 
at  a  high  elevation  angle  and  the  wavelength  is  as¬ 
sumed  to  be  much  shorter  than  the  scale  sizes  of  the 
most  intense  density  irregularities  concerned,  wave 
reflection  and  the  refraction  due  to  horizontal  strati¬ 
fication  may  be  neglected.  Only  the  forward  scatter¬ 
ing  is  considered  here. 

Let  us  first  consider  the  scattering  of  ordinary 
waves  by  density  irregularities.  The  incident  electric 
field  at  scatterer  V ,  (see  Figure  3)  is  represented  by 
Elo  =  (E„/r)  exp  [i(u>l  -  k0r)\.  For  single  forward 
scattering  with  small  angles  of  deviation,  the  resultant 
field  (E,)  of  the  scattered  waves  received  at  RiO,  O, 
Zq)  is  approximately  given  by 

dsr  A/>)  exp  [-!*,,(#•  +  R  -  Zq)I  (5) 

Rr 

where  V  is  the  total  volume  of  density  irregularities. 
The  total  field  of  ordinary  waves  received  at  R(0.  O, 
Zg)  is  composed  of  the  unscattered  and  the  scattered 
parts,  i.e., 

|E„(Zo)/^]  exp  ~  koZ0)] 

f .  Z0to2  [  A fUr)  exp  [-(*>  +  R  -  Z0)]| 

]/'  Sr  ) 

(6) 

If  the  scattered  waves  are  weak  compared  with  the 
unscattered  waves,  (6)  may  be  also  expressed  as 

E^(Zo)  »  Eto(Zo)  exp  ( (7) 

This  means  that  weak  wave  scattering  by  density  ir¬ 
regularities  introduces  a  complex  phase  shift  into  the 
original  wave.  Expanding  the  exponential  function  in 


50 


DEPOLARIZATION  OE  VHE  SATELLITE.  SIGNALS 


403 


(7)  to  the  first  order  of  <j>0  and  equating  it  to  (6)  yields  Faraday  rotation 

.  Z0  w2  f  A/J(r)  exp  ( -iku(r  +  R  —  Z0)|  <i>„  -  d>,  Z„  u»;  f  ,,  A /;(r> 

4>„  =  -1“ - 2  or - - -  (8  All  = - - -  =  ~  - - ■,  </  r  — — 

4n  c  Jv  Rr  2  4ir  <-  J  Rr 


The  phase  shift  (<J>^)  of  extraordinary  waves  due  to 
weak  scattering  by  density  irregularities  can  be  sim¬ 
ilarly  obtained  as 

.  .  to2  f  ,  A f2N(r)  exp  [  -ikjr  +  R  -  ZJ\ 

- Tr -  !9) 

The  field  of  the  linearly  polarized  waves  received  at 
R(0,  O,  Z0)  is,  therefore,  given  by 

E;  =  EJZ0)  exp  ( -  %)  +  EjZa)  exp  ( -  /'<)>,)  (10) 

which  reduces  to  (4)  in  the  absence  of  density  irreg¬ 
ularities. 

A  new  xy  coordinate  plane  is  so  chosen  that  the 
new  x  axis  coincides  with  the  polarization  vector  of 
the  unscattered  linearly  polarized  wave  at  R(0,  O, 
Z0),  i.e.,  (4)  is  expressed  as 

E„(Z0)  =  (2 E/Z0)  exp  [i(oir  -  *Z0)]f  (11) 

In  this  new  system  of  coordinates,  expression  (10) 
assumes  the  following  form 

E;(Zo)  =  (2£/Zq)  exp  [»(«/  -  kZ0) |  exp  1  — 1(4>„  +  4>,)/2] 

a  comparison  between  (11)  and  (12)  shows  that  den¬ 
sity  irregularities  cause  phase  shift,  amplitude  vari¬ 
ation,  and  Faraday  rotation  of  the  linearly  polarized 
waves.  An  induced  field  component  along  the  y  axis 
can  be  visualized  to  be  associated  with  the  induced 
Faraday  rotation  and  to  result  in  the  depolarization  of 
the  linearly  polarized  waves.  They  are  summarized 
as  follows. 


Phase  shift 


A<|>  =  Re 


Z^o,2  f  ,  A/2(r) 

4 


■  exp  |  -  iklr  +  R  -  Z„)|  sin  (» k(r  ■»  R  -  /.,,)] 
Change  of  polarization  vector 

\  2  I  4tt  c'  Rr 


■  exp  | -ikir  +  R  -  Z„)|  sin  | F*(r  +  R  -  Z„)j  ( 16) 

for  |<j>„  —  <t>,!  <  I  rad.  As  noted  in  (8)  and  (9).  d> ,  and 
<!>,  are  in  general  complex,  the  real  parts  of  which 
correspond  to  phase  shifts,  and  the  imaginary  parts 
of  which  indicate  changes  in  logarithmic  amplitudes 
of  wave  fields.  Therefore,  the  induced  Faraday  ro¬ 
tation  angles  (Ail)  and  the  polarization  vector  change 
(£/£,)  have  complex  values  as  expressed  in  (15)  and 
(16).  respectively.  Physically,  it  means  that  All  and 
( EjE ,)  result  from  the  differential  phase  shift  (i.e.. 
Re  (d>„  -  tf>,))  and  the  differential  changes  in  loga¬ 
rithmic  amplitudes  of  ordinary  and  extraordinary 
modes.  The  spatial  averages  of  A4>.  1n(E*/£R).  All. 
and  (£,/£,)  are  zeroes;  their  variances,  which  are  of 
main  physical  interest,  are  not,  however. 


3.  VARIANCES  OF  A<t>.  I n  (h'l/K,),  Alt,  AND  (A.,/6  ) 

Let  us  first  evaluate  the  induced  Faraday  rotation 
( 1 5)  or  the  change  of  polarization  vector  ( i  6)  for  |d>„ 
-  <j>x|  <  1  rad,  i.e.. 


/  =  - 


Zp  [ 
4tt  c2J¥ 


d'r 


A/y(r) 

Rr 


■  exp  | -ik(r  +  R  -  Z„)] 


•  sin  [e/t(r  +  R  -  Z,,)] 


(17) 


In  the  system  of  coordinates  (see  Figure  3)  used  to 
evaluate  (17),  r  =  (x2  +  y‘  +  z2)'  2  and  R  =  |(Z„  - 
z)2  +  .t2  +  y2]l/2.  Since  Z„,  Z„  -  z  (v2  +  v2)1'2  is 
assumed,  therefore,  r  -  z  +  (,v2  +  y2)/2z  and  R  - 
(Z„  -  z)  +  U2  +  yJ)/2 (Z„  -  z)  lead  to 

r  +  R  —  Za—  (,r2  +  y2)ZJ2z(Z„  -  z)  (18) 


•  cos  (eMr  +  R  -  Z„) \  sin  (Wr  +  R  -  Z0)]  (13) 


Amplitude  variation 


•  cos  (e*(r  +  R  -  Z0)]  cos  |*(r  +  R  -  Z0)\  (14) 


Equation  (18)  is  used  to  substitute  for  (r  +  R  -  Z0) 
in  (17),  and  Rr  in  (17)  is  simply  replaced  by  z(Z0  - 
z).  These  approximations  are  based  on  the  assumption 
that  the  outer  scale  of  density  irregularities  responsi¬ 
ble  for  the  depolarization  effect  is  small  compared 
with  the  distance  between  the  altitude  of  density  ir¬ 
regularities  and  the  receiver.  This  assumption  also 
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assures  forward  scattering  with  small  angles  of  de¬ 
viation. 

The  variance  of  Aft  or  {EJE,)  over  space  can  be, 
therefore,  approximated  as 


f  d\ 

Jv 

[*<*?  + 

y])Z« 

•  cos 

2z  ,(Z(J 

•  sin 

ek(x2  + 

y,  )Z.t) 

„  2z,(Zp 

-  Z|) 

(19) 


where  ((A/£(r,))( fN(r2))*)  is  the  autocorrelation  func¬ 
tion  of  density  irregularities.  To  carry  out  the  inte¬ 
gration  in  ( 19).  it  is  convenient  to  make  the  following 
coordinate  transformation: 


i  =  h  =  v i  -  ,v2  C  =  -  z2 

X  =  I/2U,  +  ,r.)  Y  =  1/2  (v,  +  va)  (20) 

Z  *  1/2(2,  +  2j) 


Integrating  (19)  over  X  and  Y  from  -*  to  +*  leads 
to 


,2  ru*l>MXn  rl>%CLK 


16n<  1 


di 


oj: 


^((A/i(r,))(A/v(r.))*> 


•  {2G(0)  sin  |O'(0)(4"  +  rf)1  -  G(e)  sin  |G(e)(42  +  t|2*I 
-  G(-t)  >in  1 1 ,(  —  f )( +  rpll  (21) 

where 

GO)  -  Z,Jc{2ij2'/  -  Z„)  +  4/lZZ,,  -  Z2  -  (5/2):|}‘ 1 


where 


I  ~  0,  in 

The  form  of  the  autocorrelation  function  has  to  be 
determined  for  the  further  evaluation  of  (21).  Two 
types  of  power  spectra  are  assumed  for  density  irreg¬ 
ularities:  (I)  power  law  type  and  (2)  Gaussian  type, 
which  have  been  seen  to  be  associated  respectively 
with  the  early  and  the  late  phases  of  density  irregu¬ 
larities  in  the  observation  of  satellite  scintillation 
\Basu  el  al.,  1980J. 

For  simplicity,  density  irregularities  are  assumed 
to  be  generated  by  spatially  homogeneous  processes. 
The  autocorrelation  functions,  therefore,  depend  only 


on  the  coordinate  differences  £  =  x,  -  x2,  i)  =  y, 
-  y2,  and  £  =  z,  -  z2.  Moreover,  only  statistically 
isotropic  density  irregularities  are  considered.  In  this 
case  the  autocorrelation  functions  are  only  functions 
of  (£2  +  V  +  iY1. 


3.1.  Density  irregularities  with  power  law  spectra 

For  isotropic  density  irregularities  of  power  law 
type,  the  autocorrelation  function  is  taken  to  be 

«A/«(r,))(A/y(rj))*>  =  <(A ftf) 


■  exp  [-/'(l-2  +  V  +  £2>l/2l 

where  ((A fl)2)  corresponds  to  the  mean  square  fluc¬ 
tuations  of  ionization  density  and  l  is  the  correlation 
length.  The  integration  in  (21)  over  (j  and  tj  ran  be 
conveniently  evaluated  with  the  transformation  of  (£. 
•q)  into  the  polar  coordinates  (r,  4>),  namely. 

[  di  [  dcyG  exp  |-/~'(ij2  +  r,2  +  CV/2] 

J  ~x  J  -  * 

■  sin  [G((;2  +  t)2)1  =  it  I  dt 

-'0 

•  exp  l-/-'(|G|)l/2(r  +  |C|t2)1/Jl  sin  (r)  (22) 

Although  the  general  analytic  solution  of  (22)  can¬ 
not  be  obtained,  an  approximate  form  that  is  appro¬ 
priate  for  the  case  of  kl  »  1  is  shown  in  the  appendix 
to  be 


f 

J 0 


dt  exp  l  — '(|Gj)_  ,/2(r  + 


IGICY2)  sin  (t) 


=-  exp  (-/-'ltl)(l  -  r'\G\-'n) 


(23) 


Since  G/2  ~  kl  »  1,  the  second  term  in  parenthe¬ 
ses  of  (23)  is  less  than  the  first  one  by  an  order  of 
(kl)1'2.  However,  if  the  second  term  in  (23)  is  ne¬ 
glected,  the  integration  in  (21)  over  4  and  t)  cancels 
and  thus  (|/|2>  =  0.  This  indicates  that  the  depolari¬ 
zation  of  satellite  signals,  if  any.  is  a  second  order 
ionospheric  effect. 

Using  the  identity  of  integration  given  by  (23).  (21 ) 
yields 


<|/|J>  - 


*:<(A/2)2) 
16(1  -  </2))2 


l 


d*D%ecx 


dZ 


o, 

J -Dux* 


d{l  '12(2 Z  -  Z0)/ZJt\ul 


■  (114  4  e2(ZZ0  -  Z2  -  <{/2)2)/(2 Z  - Z nlll1'2 
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4-  ||(  -  e2(ZZ„  -  Z2  -  (C/2)J)/(2 Z  -  Z„)lll/i 

-  2KI,/J}exp(-/-,KI)  (24) 

It  is  quite  clear  in  (24)  that  (|/|2)  =  0  if  the  terms 
involving  e  in  the  integrand  are  much  less  than  £.  We 
assume  that  e2[ZZ0  -  Z2  -  (4/2)2]/(ZZ  -  Z„)  ~  2nd 
is  comparable  to  (  s  /,  where  d  is  the  distance  be¬ 
tween  the  receiver  and  the  mean  altitude  of  density 
irregularities.  Then 

/  —  ltd  =  2d(fN/f)2(il/f){  I  -/2//V  cosfl  (25) 


can  be  used  to  determine  the  outer  scales  of  irregu¬ 
larities  that  are  responsible  for  the  1 36-MHz  Faraday 
polarization  fluctuations. 

Neglecting  £  in  the  braces  of  (24)  and  replacing  the 
limits  of  integration  ±D  sec  x  by  ±*,  since  D  »  /. 
we  can  easily  carry  out  the  integration  of  (24)  over 
£  and  obtain 


<l'l2>  = 


*2<(A/2)2) 

2(1 


(26) 


Because  a  D,  the  integral  in  (26)  can  be  simply 
represented  by  the  product  of  the  interval  of  integra¬ 
tion  and  the  integrand  wherein  all  Z’s  are  replaced  by 
Zm,  where  Zm  is  the  distance  between  the  satellite  and 
the  mean  altitude  of  density  irregularities.  We  finally 
obtain  the  following  result 


<|/|2>  = 


k({\N/N)2)(f2„)2 
2(1  -  (f2N))2 


( tkd)u2D  sec  x 


During  the  current  solar  maximum  period,  total 
electron  contents  (TEC’s)  as  large  as  100  units  (1  unit 
=  1016  m~2)  have  been  frequently  measured  at  As¬ 
cension  Island.  This  indicates  a  large  F,  maximum 
plasma  density  Nm  (=  100  units/250  km  ~  4. 1  x  l  o' 2 
m“\  where  250  km  is  the  assumed  slab  thickness  of 
the  ionosphere),  namely,  a  large  ambient  plasma  fre¬ 
quency  ( —  1 8  MHz).  Ionograms  obtained  at  Ascen¬ 
sion  Island  from  airborne  measurements  also  indicate 
large  f0F2  between  16  MHz  and  20  MHz  (E.  Weber 
and  J.  Moore,  private  communication,  1981). 

If/V  =  18  MHz,  Qr  =  1.5  MHz,  and  d  =  600  km 
are  used  in  (25),  the  polarization  fluctuations  of  1 36- 
MHz  satellite  signals  received  at  Ascension  Island 
(8  =  52°)  turns  out  to  be  caused  by  density  irregu¬ 
larities  with  an  ‘outer  scale  length’  of  <200  m.  It  is 
interesting  to  note  that  recent  rocket  |  Rina  el  al., 
1981|  and  satellite  in  situ  ( Basu  el  al.,  1 98 1 1  mea¬ 
surements  at  F  region  heights  show  a  break  in  spectral 
slope  generally  around  several  hundred  meters. 


Therefore,  the  outer  scale  length  so  termed  here  may 
refer  to  that  portion  of  irregularity  power  spectrum. 

At  Ascension  Island,  the  Faraday  rotation  mea¬ 
surements  were  made  at  an  elevation  angle  of  80°, 
namely,  x  =  10°  in  (27).  If  the  thickness  of  density 
irregularities  is  taken  to  be  200  km  and  the  ambient 
plasma  frequency  is  taken  to  be  18  MHz,  (27)  shows 
that  0.33%  rms  fractional  fluctuations  of  ionization 
density  can  cause  appreciable  polarization  fluctua¬ 
tions  of  1 36-MHz  signals,  i.c.,  (|/|2>l/2  =  0.1,  and  that 
1.65%  ionization  density  fluctuations  may  result  in 
large  depolarization  effect,  i.e.,  (|/|2)l/;  =  0.5. 

Similarly,  ((A4>)2)  and  (\n(Ev/E, )2)  can  be  derived 
from  ( 13)  and  ( 14),  respectively,  as 

<(A6)2>  =  4 r;/V:«A(V//V)2)X2/D  see  x/(  1  -  </;»  (28) 

and 

<1  !«(£;/£„) |2>  =  ('1<(r^>)yi  o*  x(M)'/2  (29) 

where  /  represents  the  scale  sizes  of  density  irregu¬ 
larities  that  cause  the  phase  fluctuation  of  VHF  sat¬ 
ellite  signals,  r,(~2.82  x  10  l5m)  is  the  classical 
electron  radius,  and  \  the  wavelength  of  VHF  signals. 


3.2.  Density  irregularities  with  Gaussian  spectra 
<(A/2(r,))(4/>2))*>  -  «A/2)>expl-/  2<£2  +  V  t  (’t| 


is  the  autocorrelation  function  used  to  represent  the 
density  irregularities  with  Gaussian  spectra.  In  situ 
measurements  indicate  that  Gaussian  spectra  are  as¬ 
sociated  with  the  late  developed  phase  of  irregularities 
with  a  correlation  length  of  I  km  \Basu  et  at.,  I980|. 
Therefore,  I  is  taken  to  be  1  km  in  (30).  Substituting 
(30)  into  (21)  and  integrating  with  respect  to  f  and 
t),  we  obtain 

«>-!»/ 

'a  *  * 

■cxp(-/  2{"){2F(0)  -  Fir)  -  F(-n) f  (31) 

where 


where  t  =  0,  ±e. 
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Since  2£(2Z  —  Zj/Z^kl'  is  of  the  order  of  2/kl  ~ 
7.0  x  10  4  and  4 eZ(Z0  -  Z)/Zkl2  of  the  order  of 
4 ed/kl2  -  5.0  x  10  \  FiO)  and  F(±e)  have  the  gen¬ 
eral  expression  of  (1  +  8:)  1  with  8  1.  which  can 

be  well  approximated  as  (1  -  82).  The  expression  in 
the  braces  of  (31)  is,  consequently,  simplified  to  be 
2[4eZ(Z0  -  Z)/Zkl2\2.  Straightforward  integration  of 
(3 1 )  over  £  and  Z  leads  to 

<|/|2>  =  4~V/2<(AN//V)j><1  -  fl)hfl)*(Cl)2d2l  'D 


•  sec  x  cos2  6  (33) 

Substituting  into  (33)  the  same  parameters  as  used  in 
the  case  of  power  law  type  irregularities  shows  that 
unrealistically  large  density  fluctuations  are  required 
for  causing  the  polarization  fluctuations  of  136-MHz 
signals,  namely  {(hN/N)2)'12  =  843%  for  <|/|2>,/2  = 
0.1  and  ((A N/Nf)'n  =  4215%  for  <|/|2>  =  0.5.  A  rea¬ 
sonable  ionospheric  density  depletion  (((kN/N)2)'12) 
with  Gaussian  spectra  is  about  a  few  percentage 
points,  say  5%.  According  to  (33),  the  ambient 
plasma  frequency  must  be  as  high  as  65  MHz  so  that 
the  Gaussian  type  of  irregularities  is  able  to  cause  the 
appreciable  depolarization  (i.e.,  <|/|2>"2  =  0.1)  of 
136-MHz  signals.  Such  a  high  ionospheric  plasma 
frequency  is,  however,  unrealistic  under  the  normal 
conditions. 

The  derivation  of  <( A<J>)?>  and  <1  n(E'R/E„)2)  can  be 
easily  done  in  the  case  of  Gaussian  type  irregularities. 
They  are 


«A4>)  >  = 

and 


,v  _  irl/2/X2r2yV2((A)V/)V)2) 


4(1  -</i))2 


{7(0)  +  7(1)} 


(34) 


(IN  £/»/£*)!>  = 


2.  _  vU2lK2r]N2{UkN/NY) 


4(1  -  </.v»J 


{7(0)  -7(1)} 


(35) 


where 


7<v) 


-t 


U4/(Z„  -  /)/Z„kl2 12} 


where 


a  -  0.  I 

Equations  (34)  and  (35)  can  be  reduced,  respectively, 
from  the  (26)  and  (27)  of  Yeh  and  Liu  1 1967),  wherein 
7(e)  reduces  to  7(0),  since  e  —  10  4  for  the  136-MHz 
wave  propagation. 


4  DISCUSSION  AND  CONCLUSION 

In  this  paper,  wave  scattering  by  density  irregular¬ 
ities  with  power  law  spectra  has  been  proposed  as  the 
process  for  cat:  ing  the  intense  and  fast  fluctuations 
of  the  SIRIO  136-MHz  signals  received  at  Ascension 
Island.  Based  on  the  above  analysis,  the  following 
things  can  be  understood:  (1)  the  nature  of  density 
irregularities  that  result  in  the  polarization  fluctua¬ 
tions  of  satellite  signals,  (2)  why  the  polarization  fluc¬ 
tuations  of  136-MHz  signals  has  a  close  correlation 
with  the  occurrence  of  the  L  band  scintillation,  (3) 
why  the  polarization  fluctuations  of  136-MHz  signals 
has  only  been  seen  at  locations  near  the  equatorial 
anomaly  crests,  and  (4)  the  frequency  limit  beyond 
which  density  irregularities  will  not  affect  signifi¬ 
cantly  the  polarization  vector  of  satellite  signals. 

We  have  already  shown  that  the  polarization  fluc¬ 
tuations  of  136-MHz  satellite  signals  result  from  the 
wave  scattering  by  density  irregularities  whose  power 
spectra  exhibit  power  law  rather  than  Gaussian  dis¬ 
tribution.  In  situ  measurements  [ Basu  el  al.,  1980] 
show  that,  while  density  irregularities  have  relatively 
flat  power  spectra  obeying  a  power  law  slope  in  their 
early  phase,  density  irregularities  in  their  late  phase 
are  seen  to  have  a  sharp  slope  essentially  character¬ 
istic  of  a  Gaussian  distribution.  This  indicates  that, 
during  the  change  of  power  spectra  from  power  law 
to  Gaussian  types,  relatively  short-scale  (say,  less 
than  I  km)  density  irregularities  decay.  As  shown  in 
(25),  the  polarization  fluctuations  of  136-MHz  satel¬ 
lite  signals  turns  out  to  be  caused  by  those  short-scale 
density  irregularities  with  an  outer  scale  length  less 
than  about  200  m. 

Density  irregularities  with  scale  lengths  of  a  few 
hundreds  of  meters  are  responsible  for  the  L  band 
scintillation.  Therefore,  the  L  band  scintillation  is  not 
expected  during  the  late  phase  of  density  irregularities 
characterized  with  Gaussian  spectra,  indicating  the 
decaying  of  short-scale  density  irregularities.  Obser¬ 
vations  of  the  257-MHz  scintillation  at  Ascension  Is¬ 
land  in  the  presence  and  absence  of  the  L  band  scin¬ 
tillation,  indeed,  indicate  the  characteristic  increasing 
of  the  correlation  lengths  of  the  257-MHz  signals 
[Basu  el  al.,  1981 J.  That  both  the  L  band  scintillation 
and  the  polarization  fluctuations  of  the  136-MHz  sig¬ 
nals  are  caused  by  relatively  short-scale  density  irreg¬ 
ularities  is  likely  to  account  for  the  coexistence  of 
these  two  phenomena. 

The  ambient  plasma  density  (TV)  and  the  propaga¬ 
tion  angle  (0)  are  the  most  location  dependent  param- 
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Fig.  4.  Fractional  fluctuations  ({(AA//A/)2)1")  of  the  plasma  density  N  versus  the  outer  scale  lengths  /  of  density 
irregularities  or  the  ambient  plasma  density  N,  as  required  to  cause  the  polarization  fluctuations  of  1 46  MHz 
satellite  signals  at  Ascension  Island.  Scale  (A)  for  {|fi,),/‘  =  0.1,  scale  (B)  for  (|/|>)1'1  =  1.0. 


eters.  As  mentioned  before,  the  ambient  plasma  den¬ 
sity  measured  at  locations  near  the  Appleton  anomaly 
crests  may  be  larger  than  those  at  other  equatorial  lo¬ 
cations  by  a  factor  of  2-4  during  the  solar  maximum 
period  [Rastogi,  1966;  Aarons  el  at.,  1981],  Let  us 
look  at  the  possibility  of  observing  the  1 36-MHz  po¬ 
larization  fluctuations  near  the  magnetic  equator, 
e.g..  Natal,  Brazil  (geographic  coordinate  5.8S°S, 
35.23°W,  9.6°S  dip). 

The  propagation  angle  at  Natal  is  about  80°  (cf.  52° 
at  Ascension  Island).  If  the  ambient  plasma  frequency 
at  Natal  is  taken  to  be  10  MHz  (cf.  18  MHz  at  As¬ 
cension  Island),/*  cos  8  obtained  at  Natal  is  less  than 
that  at  Ascension  Island  by  a  factor  of  11.5.  There¬ 
fore,  according  to  /  «  /*  cos  8  in  (25),  the  density 
irregularities  causing  the  polarization  fluctuations  of 
136  MHz  signals  at  Natal,  if  any,  have  an  outer  scale 
shorter  than  20  m  (cf.  200  m  at  Ascension  Island). 
On  the  other  hand,  since  the  ratio  of  /*( cos  0)172  at 
Natal  to  that  at  Ascension  Island  is  35.6,  (|/|2)  * 


((AA//A/)2)/*(cos  8)l/2  in  (27)  requires  that  ((AN/A/)2) 
at  Natal  be  larger  than  that  at  Ascension  Island  by 
35.6  times  to  cause  the  same  level  of  1 36-MHz  po¬ 
larization  fluctuations,  i.e.,  ((AN/A/)2)172  =  10%  for 
(|/)2)172  =  0.5.  This  example  shows  clearly  that  low 
ambient  plasma  densities  and  large  propagation  an¬ 
gles  at  the  magnetic  equator  make  it  impossible  to 
observe  the  polarization  fluctuations  of  1 36-MHz  sat¬ 
ellite  signals  at  equatorial  locations  except  near  the 
anomaly  crests. 

In  (25),  /  */*  *  N  indicates  that,  in  an  environment 
of  high  ambient  plasma  density,  density  irregularities 
with  longer  scale  lengths  can  also  be  involved  in  pro¬ 
ducing  the  polarization  fluctuations  of  satellite  sig¬ 
nals.  Eliminating  fN  from  (25)  and  (27)  yields 

((AN/N)2)172  =  C(\r\2)'nl  1/4  (36) 

where  C  =  4(£l/f)dk~  ,/4(0  sec  x)  1/2  cos  8  -  16.5 
for  the  assumed  ionospheric  condition  at  Ascension 
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Fig.  5.  Plot  of /(f)  =  exp  |-/  '(iG|)  '"(/  +  iGlt1)'71}  sin  (f)  for 
Gl1  »  I.  The  unhatched  area  represents  approximately  the  inte¬ 
gration  of  /(f)  over  f  from  0  to 

Island:  d  =  600  km,  D  =  200  km,  fl  =  1.5  MHz, 
/=  136  MHz,  6  =  52°,  and  \  -  10°.  It  is  clear  from 
the  plotting  of  (36)  in  Figure  4  that,  in  the  presence 
of  density  irregularities,  high  ambient  plasma  density 
provides  the  favorable  environment  for  the  VHF  po¬ 
larization  fluctuations  to  occur.  This  may  explain  why 
the  polarization  fluctuations  of  136-MHz  signals  has 
been  frequently  observed  at  Ascension  island  during 
the  current  solar  maximum  period. 

For  illustration,  the  frequency  limit  beyond  which 
density  irregularities  cannot  cause  significant  polari¬ 
zation  fluctuations  of  satellite  signals  is  loosely  de¬ 
fined  as  that  at  which  (|/|2>,/2  =  0.1  requires  1%  of 
density  fluctuation  under  the  aforesaid  ionospheric 
condition.  This  frequency  is  determined  as  235  MHz 
from  (27),  and  the  outer  scale  length  of  density  ir¬ 
regularities  is  determined  as  28  m  from  (25). 

It  should  be  pointed  out  that  the  theory  of  singly 
wave  scattering  discussed  in  this  paper  cannot  explain 
the  observed  saturated  amplitude  scintillations  of  1 36- 
MHz  signals.  It  is  well  known  that  the  saturation  of 
VHF  amplitude  scintillation  is  attributed  to  multiple 
scattering.  Therefore,  deductions  of  the  present  the¬ 
ory  for  the  amplitude  scintillation  may  not  be  true. 
The  validity  of  single  scatter  approximation  as  applied 
to  the  problem  of  polarization  fluctuations  can  be  jus¬ 
tified  \Lee,  198 1 J.  Multiple  scattering  is  expected  to 
give  rise  to  the  saturation  of  polarization  fluctuations 
in  the  presence  of  large  density  irregularities.  How¬ 
ever.  the  effect  of  multiple  scattering  does  not  impose 
significant  quantitative  changes  in  the  results  derived 
from  the  theory  of  single  wave  scattering  Simple 
analyses  presented  in  this  paper  turn  out  to  be  ade¬ 
quate  for  the  understanding  of  the  136-MHz  polari¬ 
zation  fluctuations  and  associated  phenomena. 


APPENDIX  AN  APPROXIMATE  SOLUTION 
OF  dt  exp  l -/  '(|(i|)  '  :(t  +  \G  {’)'  j  mu  (/) 

For  Gl 2  ~  kl,  which  is  much  greater  than  one 
as  assumed,  the  integrand  in  (22),  i.e.,  exp 
l-r'(\G\r'n(t  +  |Gj£2)'/2}  sin  (f)  can  oscillate  for 
many  cycles  before  ii  dies  out.  As  illustrated  in  Figure 
5,  only  the  unhatched  area  does  not  vanish  and  rep¬ 
resent  the  approximate  integration  of  (22). 

Integration  by  parts  twice  in  (22)  yields 


r exp  l-r'(|C|)-,/2(f  +  |G|£2),/2]  sin  (r) 


=  exp(-r'|{|) 


-f  dt  sin  (/) 


■  exp  [-r'(|G|)-|/2((  +  |G|£2),/2]  (Al) 

Since  exp  [-r'(|G|rl/2tr  +  |G|£2)I/2]  vanishes  effec¬ 
tively  at  distances  of  the  order  of  the  correla¬ 
tion  length,  the  double  differentiation  of  exp 
l-/"'(|G|)~,/2(r  +  |G|C2),/21  with  respect  to  t  may  be 
approximated  by 

max  {exp  |-/  '(|G{)*  !'2(f  +  |G|£V'2]}/G2/4 

=  exp(-r'|t|)/G2/4  (A2) 

Sin  (t)  in  ( A 1 )  may  be  represented  approximately  by 
the  parabolic  function  (t)'u.  The  integral  in  (Al)  can 
therefore  be  approximated  as 

|  dt  sin  (r)-^pexp  |  -/~'(iC|)~  'n(t  +  |G|£J)I/2] 

-  |G|/2  •  (|G|/2),/2  •  expt-r'Itb/G2/4 

=  exp  (~r'\C)/(\G\)<nl  (A3) 

Substituting  (A3)  into  (Al)  leads  to  the  approxi¬ 
mate  form  of  (22)  for  Gl2  »  1,  namely, 


Jo 


sin  (r)  exp  |-r'(|G|)',/2[r  +  |C|(2]' 


=»  exp  (-r'lSDll  -r'(\G\)->n]  (A4) 

Indeed,  that  the  approximate  form  in  (A4)  carries  the 
right  functional  relationship  with  G  and  /  is  confirmed 
by  the  exact  integration  in  (22)  for  £  =  0. 
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VHF  Faraday  polarization  fluctuations 
and  strong  /.-band 
amplitude  scintillations 
near  Appleton  anomaly  crests 
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The  Faraday  rotation  angle  of  a  linearly  polarized  transionos- 
pheric  very  high  frequency  (VHF)  signal  near  the  crests  of  the 
Appleton  anomaly,  around  the  recent  solar  maximum  period, 
often  exhibits  fast  and  intense  fluctuations,  when  the  amplitude 
of  the  signal  shows  saturated  and  very  fast  scintillations' These 
fluctuations  usually  occur  during  the  post-sunset  hours  and 
decay  around  midnight.  The  occurrence  pattern  has  been  found 
to  be  similar  to  that  of  microwave  amplitude  scintillations  at 
Z-band.  We  present  here  certain  observed  characteristics  of 
the  fast  polarization  fluctuations  and  report  the  results  of  an 
experiment,  which  shows  that  there  is  a  loss  of  correlation 
between  fadings  on  the  two  circularly  polarized  component 
modes  of  the  signal  during  the  periods  of  VHF  polarization 
fluctuations.  Hence,  the  usual  concept  of  Faraday  rotation  is 
no  longer  valid  during  these  limes.  The  above  features  are 
discussed  in  terms  of  our  present  knowledge  of  small  scale 
irregularity  structures  in  the  equatorial  ionosphere. 

Equatorial  ionospheric  irregularities  have  been  extensively 
investigaied  during  the  pasl  decade  by  techniques  such  as  rocket 
and  satellite  in  stiu  measurements,  radar,  airglow  and  propaga¬ 
tion  effects  on  iransionosphcric  signals'.  It  has  been  shown  that 
the  plumes  of  equatorial  irregularities  associated  with  deple¬ 
tions  in  background  electron  content  can  he  studied  in  a  very 
simple  and  inexpensive  manner  from  simultaneous  measure¬ 
ments  of  amplitude  scintillations  and  Faraday  rotation  of  a 
VHF  signal  propagating  through  the  ionosphere4 While  the 
occurrence  of  depletions  of  the  integrated  electron  content  is 
commonly  observed  with  Faraday  rotation  records  obtained 
near  the  magnetic  equator"  ’,  a  new  feature,  fast  and  intense 
fluctuations  of  the  Faraday  rotation  angle,  has  recently  been 
reported  from  stations  situated  near  the  crests  of  the  Appleton 
anomaly' :,  a  region  characterized  by  severe  amplitude  scintilla¬ 
tions  during  solar  maximum  years"  '.  Polarization  fluctuations 
of  136-MHz  signals  have  been  observed  al  Ascension  Island 
<31°S  dip),  which  is  siluated  close  to  the  southern  crest  of  the 
equatorial  anomaly  in  Ihc  African  zone,  and  at  Calcutta  (32°N 


dip)  and  Ahmedahad1  1  i34°N  dipl  near  the  northern  crest  of 
ihc  anomaly  in  the  Indian  sector.  These  have  sometimes  been 
observed  even  al  Delhi  |4()°N  dip)  i  V.  V.  Somayajulu,  personal 
communication)  near  the  northern  edge  of  the  anomaly,  but 
no  such  phenomenon  has  been  observed  at  Arequipa  (9°Sdip), 
Peru  or  at  Natal  tl()®S  dip),  Brazil  which  are  close  to  the 
magnetic  equator.  Kaushika  and  deMcndonca"  possibly  obser¬ 
ved  the  same  type  of  fluctuations  in  the  polarization  records  of 
137-MHz  signals  recorded  at  Sao  Paulo  (26°  S  dip)  which  is 
also  situated  near  the  crest  of  the  anomaly,  but  they  had  no 
amplitude  scintillation  measurements.  In  their  attempt  to 
associate  the  fluctuations  with  atmospheric  gravity  waves,  they 
were  concerned  with  slower  fluctuations  of  period  2  min  and 
longer. 

Figure  I  shows  a  typical  polarimeter  record  obtained  during 
post-sunset  hours  al  Ascension  Island  (7.9°  S,  14.4°  W,  31°  S 
dip).  The  channel  at  the  top  is  the  amplitude  of  the  136-MHz 
transmission  from  the  geostationary  satellite  SIRIO  (IS°W) 
and  the  two  bottom  channels  are  the  Faraday  polarization  angle 
of  the  signal  with  one  shifted  from  the  other  by  90°.  The 
polarization  fluctuations  are  often  too  intense  to  render  the 
Faraday  data  useful.  The  station  also  recorded  signals  from  the 
satellite  Marisat  (also  at  15°  W)  at  257  MHz  and  1.54  GHz  for 
amplitude  scintillation  measurements.  The  sub-ionospheric 
points  of  the  two  satellites  arc  essentially  identical.  The  third 
channel  from  the  top  shows  amplitude  scintillations  of  the 
1.54-GHz  signal.  Polarization  fluctuations  are  found  when  the 
amplitude  of  the  VHF  signal  shows  saturated  and  very  fast 
scintillations,  which  generally  occur  during  the  pre-midnight 
hours.  We  have  found  during  more  than  2  yr  of  observations 
around  the  recent  solar  maximum  period  that  the  polarization 
fluctuations  coexist  with  intense  /.-band  scintillations  at 
1 .54  GHz.  During  the  post-midnight  period  polarization  fluctu¬ 
ations  decay,  even  though  VHF  amplitude  scintillation  remains 
saturated,  hut  with  slower  fading  rates. 

The  change  in  fading  rates  of  VHF  amplitude  scintillations, 
with  and  without  polarization  fluctuations,  is  clearly  evident  in 
the  power  spectra  of  136-MHz  scintillation  measurements  as 
presented  in  Fig.  2 a  and  6  respectively.  The  corresponding 
power  spectra  obtained  at  257  MHz  and  1.54  GHz  are  also 
shown.  The  approximate  corner  frequencies  of  the  power 
spectra  are  indicated  by  arrows.  The  corner  frequency  of  a 
power  spectrum  gives  a  measure  of  the  fading  rale,  a  higher 
corner  frequency  corresponding  to  a  faster  fading  rate.  In  the 
present  case,  a  lowering  of  the  corner  frequency  is  obvious  in 
the  late  phase,  as  shown  in  Fig.  26.  when  no  polarization 
fluctuations  are  observed. 

Power  spectrum  analysis  of  scintillation  measurements  is 
generally  used  to  obtain  information  on  the  characteristics  of 
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the  ionospheric  irregularities  responsible  for  scintillations. 
However,  the  interpretations  of  the  power  spectrum  in  condi¬ 
tions  of  weak  and  strong  scintillations  involve  different  scatter¬ 
ing  mechanisms  by  the  ionospheric  irregularities.  According  to 
the  theory  of  radiowave  propagation,  for  weak  scintillation 
conditions,  the  amplitude  fading  rate  of  a  transionospheric 
signal  from  a  geostationary  satellite  is  proportional  to  i'll,, 
where  v  is  the  drift  velocity  of  the  irregularities  and  I,  the 
Fresnel  dimension*'.  Because  lf  oe>/A  (where  A  is  the  radiowave 
length),  the  fading  rate  should  increase  with  observing 
frequency  in  weak  scatter  conditions.  In  contrast,  in  the  ease 
of  severe  scintillations,  the  wave  propagation  has  to  be  con¬ 
sidered  in  the  framework  of  multiple  scatter  theory  and  the 
recorded  scintillation  data  pertain  to  the  ground  diffraction 
pattern,  which  cannot  be  related  to  the  irregularity  pattern  in 
the  ionosphere  in  a  simple  manner1’.  The  fading  rate,  in  severe 
scintillation  conditions,  is  inversely  related  to  the  autocorrela¬ 
tion  distance  of  the  ground  pattern.  The  autocorrelation  interval 
is  determined  by  the  drift  velocity  as  well  as  by  the  strength  of 
scattering.  As  the  strength  of  scattering  (in  other  words,  scin¬ 
tillation  index)  increases,  the  autocorrelation  distance 
decreases",  and  an  increase  in  the  fading  rate  is  observed". 
The  scintillation  index  is  usually  inversely  related  to  the  observ¬ 
ing  frequency.  Thus,  in  the  case  of  severe  scintillations,  with 
an  increase  of  the  observing  frequency  the  autocorrelation 
distance  should  increase  and  hence,  a  slowing  down  of  the 
amplitude  fading  rate  occurs 

In  the  equatorial  region  amplitude  scintillations  observed  at 
136  MHz  and  250  MHz  arc  usually  severe  isaturaled  or  near 
saturated)  for  most  of  the  time  the  phenomena  arc  observed. 
Near  the  crests  of  the  Appleton  anomaly,  severe  scintillations 
in  excess  of  20  dB  are  quite  frequent  even  at  /  -band  (such  as 
1.54-GHz  observations  at  Ascension  Islandl  during  the  pre¬ 
midnight  hours,  particularly  around  the  recent  solar  maximum'. 
The  power  spectra  depicted  in  Fig.  2  arc  representative  of 
scintillations  observed  at  the  corresponding  frequencies  in  the 
developed  and  late  phases  of  equatorial  night-time  ionospheric 
irregularities  respectively. 

During  an  intensive  campaign  of  measurements  in  December 
1979,  in  order  to  investigate  the  nature  of  fast  polarization 
fluctuations,  the  left  and  right  circular  components  of  the  136- 
MHz  signal  were  recorded  separately  in  addition  to  the  routine 
observations.  The  data  were  recorded  on  a  magnetic  tape  in 
analog  form,  which  was  digitized  at  an  IX  samples  s  '  rate  for 
later  analysis.  The  cross-correlation  between  the  amplitude 
fadings  of  the  two  component  modes  was  computed  over  3-min 
intervals.  The  scintillation  index  S4  irefs  15,  lb),  which  is  a 
measure  of  the  r.m.s.  amplitude  fluctuations  about  the  mean 
level,  is  also  computed  for  each  3  min  interval  at  different 
observing  frequencies.  Figure  3  shows  the  variation  of  the 
cross-correlation  coefficient  between  the  amplitude  fadings  of 
the  two  circularly  polarized  component  waves  at  1 36  MHz  with 
changes  in  the  scintillation  index  S4  observed  at  1.54  GHz. 

The  correlation  coefficient  decreases  as  St  increases  and 
attains  low  values  (<0.2)  when  ,V4  is  0.5  or  more.  The  intervals 
with  St  >0.5  at  1.54  GHz  generally  exhibit  fast  polarization 
fluctuations  at  136  MHz.  In  other  words,  the  two  circularly 
polarized  components  undergo  quasi-independent  fluctuations 
when  strong  /.-band  and  polarization  fluctuations  are  obtained. 
If  the  two  modes  become  quasi-independent,  the  usual  concept 
of  Faraday  rotation  breaks  down.  The  close  correlation  of 
occurrence  of  1 36-MHz  polarization  fluctuations  and  the  I  54- 
GHz  scintillation,  clearly  seen  in  Fig.  3,  implies  that  either 
irregularities  of  nearly  the  same  scale  size  are  responsible  for 
the  two  phenomena  or  that  the  irregularities  causing  the  two 
phenomena  coexist  in  the  ionosphere. 

l.ee  ft  of.17  have  shown  that  ionospheric  power  law  density 
irregularities  with  an  outer  scale  size  in  the  range  50-2lH)m 
are  responsible  for  the  1 36-MHz.  polarization  fluctuations. 
Faraday  polarization  fluctuations  result  from  differential  phase 
shift  and  the  differential  change  in  the  logarithmic  amplitude 
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D*.  2  Power  spectra  of  amplitude  scintillations  at  different 
frequencies  for  a  V  min  sample  centred  on  the  lime  indicated,  rt. 

2 1  h  46  min  30  s  I  it;  Faraday  polarization  channels  show  intense 
and  rapid  fluctuations  during  this  period  h,  22  h  -i  nun  Alls  u  i : 
no  fast  polarization  fluctuation  was  observed  during  this  period 
I  be  approximate  corner  frequencies  are  given  bv  arrows. 

of  the  two  circularly  polarized  modes  introduced  by  wave  scat¬ 
tering  from  these  short-scalc  density  irregularities.  The  theory 
shows  that  a  deviation  of  -  I  -2%  in  a  high  ambient  ionization 
density,  say,  not  less  Ihnn  3.0  x  Id' 'in  corresponding  to  a 
plasma  frequency  of  ovci  IS  Mil/,  is  required  to  provide  a 
favourable  environment  for  the  I  36-MHz  fluctuations  to  occur. 
When  the  r.m.s.  differential  phase  departure  exceeds  — (3.5  rad. 
the  two  circularly  polarized  modes  are  expected  to  become 
practically  uncorrclatcd  Scintillations  at  1 .54  GHz  are  most 
effectively  produced  by  irregularities  with  scale  sizes  of  a  few 
hundred  metres  ('300  m  in  ihc  present  case),  which  arc  of  the 
same  order  as  the  Fresnel  zone  dimension.  Both  the  1 36-MHz 
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r*  1  Variation  of  cross-correlation  coefficients  of  I  he  fluctu¬ 
ations  of  the  two  opposite  circularly  polarized  components  of  the 
I.Vi-MII?  signal  with  the  amplitude  scintillation  hides  .V4  at 
I  54 fill/.  The  intervals  ,V4  -0.5  correspond  to  intense  and  fast 
I'aradav  polari/alion  lluetuations. 

polarization  fluctuations  and  the  1.54-Gllz  scintillations 
involve  density  irregularities  of  a  few  hundred  metres  scale  size. 
If  (he  irregularity  power  spectrum  follows  a  power  law,  with 
diminishing  spectral  density  at  shorter  scale  lengths,  and  if  the 
overall  strength  of  the  irregularities  is  high  enough  to  have 
spectral  density  in  the  50-200  m  scale  size  range  sufficient  to 
cause  polarization  fluctuations  at  1 36  MHz,  it  is  implied  that 
irregularities  at  a  ^300  m  scale  size  range  will  also  have  enough 
spectral  power  to  produce  appreciable  amplitude  scintillations 
at  1 .54  GHz.  Thus,  the  presence  of  fast  polarization  fluctuations 
at  136  MHz  may  be  taken  as  an  indication  of  the  possible 
scintillations  at  /.-hand  and  lower  frequencies.  The  reverse  may 
not,  however,  always  he  true. 


From  multi-technique  observations  Basu  et  al have  repor¬ 
ted  that,  while  during  the  pre-midnight  period  the  power  spectra 
of  equatorial  irregularities  follow  a  power  law,  during  post¬ 
midnight  hours  the  overall  strength  of  the  irregularities  is 
eroded  and  there  is  a  transition  to  a  quasi-gaussian  type  of 
spectrum  indicating  the  decay  of  shorter-scale  (typically  <  1  km) 
irregularities.  As  the  1.54-GHz  scintillation  became  weak,  with 
the  decay  of  short-scale  density  irregularities,  the  two  circularly 
polarized  modes  were  seen  to  have  high  cross-correlation 
coefficients,  reflecting  no,  or  very  small,  polarization  fluctu¬ 
ations. 

Except  in  the  HF  band  (for  example,  20  MHz),  the  polariz¬ 
ation  fluctuations  of  satellite  signals  cannot  be  effectively  caused 
hy  irregularities  with  quasi-gaussian  spectra19.  Near  the 
Appleton  anomaly  crests,  in  solar  maximum  conditions,  the 
ambient  ionization  density  (n )  is  very  high  and  the  high  level 
persists  during  the  evening  hours.  As  a  result,  the  absolute 
value  of  the  irregularity  intensity  (An)  is  also  very  high  in  the 
anomaly  region".  In  contrast,  in  regions  near  the  magnetic 
equator,  the  relatively  small  n  and  An  and  the  relatively  large 
angles  of  propagation  with  respect  to  the  geomagnetic  held 
make  it  impossible  to  give  rise  to  significant  polarization  fluctu¬ 
ations  of  136-MHz  signals'7. 
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